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ABSTRACT 
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A mathematical model for use in real-time piloted simula- 
tion of a 1985-era tilt rotor passenger aircraft is presented. 
The model comprises the basic six degrees-of- freedom equations 
of motion, and a large angle of attack representation of the 
airframe and rotor aerodynamics, together with equations and 
functions used to model turbine engine performance, aircraft 
control system and stability augmentation system. 

A complete derivation of the primary equations is given 
together with a description of the modeling techniques used. 
Data for the model is included in an Appendix. 
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SUMMARY 
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This report documents the equations , functions , control 
systems diagrams, and data required in the real-time simulation 
of a 1985-era tilt rotor passenger aircraft. The simulatio*' 
mathematical model was intended for use on the NASA-Ames 7 , 

Simulator for Advanced Aircraft to study the handling qual < 
of large tilt rotor aircraft. The model could also be ust ' :o 
research on advanced terminal area control. :rs. 


The mathematical model consists of the rigid body equations 
for motion of the aircraft in roll, pitch, and yaw about a 
moving center of gravity. The equations differ somewhat from 
the classical equations because of the necessity of accounting 
for the motion of the tilting rotors and nacelles. The math 
model is " full- force" , that is, the representation of the aero- 
dynamics of tlie rotors and airframe is suitable for the large 
angles of attack encountered in VSTOL flight and can represent 
pure rearwards and sidewards flight from hover. The aero- 
dynamics of the airframe and the interference between components 
was estimated from a combination of theory and experimental 
data. The forces and moments acting on the large 5 6- foot dia- 
meter hingeless rotors were obtained from a regression analysis 
of test data on a smaller rotor of similar construction and 
properties. The control system models p:' lot controls, a thrust 
management system and a stedDility augmentation system. 
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LIST OF SYMBOLS 
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Symbol 

Definition 

Units 

A 

Rotor disc area (per rotor) 

ft2 

AR 

Aspect ratio, b^/s 

ND 

Ale 

Lateral cyclic angle in rotor wind 
axes 

deg 

Aic 

Lateral cyclic angle in swashplate 
axes 

deg 


Lateral cyclic amgle in swashplate 
axes resolved through swashplate 
phase angle 

deg 

a 

Speed of sound or acceleration 

ft/sec or 
ft/sec2 

a 

Acceleration 

ft/sec2 

(ag/a) 

Ratio of lift-curve slope in ground 
effect to lift-curve slope out of 
ground effect 

ND 

ao^a32 

Coefficients in wing lift and drag 
equations 

— 

Bq 

Percent brake pedal deflection 

ND 

B.L. 

Aircraft butt line 

in. 

Blc 

Longitudinal cyclic angle in rotor 
wind axes 

deg 

Bic 

Longitudinal cyclic angle in swash- 
plate axes 

deg 

O »* 

®lc 

Longitudinal cyclic angle in swash- 
plate axes resolved through swash- 
plate phase angle 

deg 

b 

Span of lifting surface (wing, tail, 
etc) 

ft 

c 

Chord 

ft 

Cd 

Drag coefficient, ^ 

qS 

ND 
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Symbol 

Definition 

Units 

CDo 

Drag coefficient at zero lift 

ND 

ACd 

Drag coefficient increment 

ND 

^DS 

Drag coefficient referred to rotor 
slipstream dynamic pressure/ D/qgS 

ND 

Cl 

Lift coefficient/ L/qS 

ND 

CLo 

Average lift coefficient 

ND 

ACl 

Lift coefficient increment 

ND 

CLs 

Lift coefficient referred to rotor 
slipstream dynamic pressure/ L/qsS 

ND 

CLq 

Li ft- curve slope 

1/rad 

CL{ 

Lift increment due to flap deflection 

1/deg 


Rolling moment coefficient, '^/q bS 

ND 


Rolling moment coefficient referred 
to rotor slipstream dynamic pressure, 
i/qgbS 

ND 

Cm 

Pitching moment coefficient, M/qSc 

ND 

CMo 

Wing pitching moment coefficient as a 
function of flap deflection; pitching 
moment coefficient of fuselage or 
nacelles at zero angle of attack 

ND 

A Cm 

Pitching moment coefficient increment 

ND 

CMs 

Pitching moment coefficient referred 
to rotor slipstream dynamic pressure, 
M/qgSc 



Change in wing/body pitching moment 
coefficient as a function of flaperon 
deflection 

ND 

Cn 

Yawing moment coefficient, N/qSb 

ND 

CNq 

Yawing moment coefficient of fuselage 
or nacelles at zero angle of attack 

ND 
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Symbol Definition Units 

CNs Yawing moment coefficient referred to ND 

rotor slipstream dynamic pressure/ 

N/qgSb 

Cnf Rotor normal force coefficient/ ND 

NF/pTTfi2R4 

%Fo Rotor normal force coefficient with ND 

zero cyclic pitch 

Cp Rotor power coefficient/ 550RHP ND 

Cpq Rotor power coefficient with zero ND 

cyclic pitch 

CpM Rotor hub pitching moment coefficient/ ND 

PM/pTTfi2R5 

Cpf 4 Rotor hiib pitching moment coefficient ND 

° with zero cyclic pitch 

CsF Rotor side force coefficient, ND 

SF/pTTfi2R4 

^SFn Rotor side force coefficient with zero ND 

cyclic pitch 

C<j Rotor thrust coefficient/ T/piTfi2R4 jjq 

CTo Rotor thrust coefficient with zero ND 

cyclic pitch 

Cxg Rotor thrust coefficient referred to ND 

rotor slipstream dyna»ric pressure, 

T/qgA 

Cy Side force coefficient, Y/qS ND 

Cyj 4 Rotor yawing moment coefficient, ND 

P TTfi^RS 

CyMq Rotor yawing moment coefficient with ND 

zero cyclic pitch 

CYjj Lift-curve slope of vertical tail 1/rad 

Cq Coefficient of equation that defines ND 

pitching moment coefficient as a 
. function of flap deflection 
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Symbol 

Definition 
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Units 

Cl 

Coefficient of equation that defines 
pitching moment coefficient as a 
function of flap deflection 

1/rad 

C2 

Coefficient of equation that defines 
pitching moment coefficient as a 
function of flap deflection 

l/rad2 

D 

Rotor diameter 

ft 

(D/T) 

Aircraft download-to-thrust ratio 

ND 

DNF1-.5 

Coefficients in the equation for the 
change in normal force coefficient 
with lateral cyclic angle 

1/deg 

DpMi^e 

Coefricients in the equation for the 
change in hub pitching moment coeffi- 
cient with lateral cyclic angle 

1/deg 

Dsfi^s 

Coefficients in the equation for the 
change in side force coefficient with 
lateral cyclic angle 

1/deg 

“STn 

Damping coefficients of the landing 
gear oleo struts 

Ib/ft/sec 


Coefficients in the equation for the 
change in hub yawing moment coeffi- 
cient with lateral cyclic angle 

1/deg 

dCjgp/dAic 

Change in normal force coefficient 
with lateral cyclic angle 

1/deg 

<^Cnf/<3Biq 

Change in normal force coefficient 
with longitudinal cyclic angle 

1/deg 

dCpM/dAic 

Change in hub pitching moment coeffi- 
cient with lateral cyclic angle 

1/deg 

dCpM/dBic 

Change in hub pitching moment coeffi- 
cient with longitudinal cyclic angle 

1/deg 

dCpfi/ dQ 

Change in hub pitching moment coeffi- 
cient with pitch rate 

1/rad/sec 

dCsF/dAic 

Change in side force coefficient with 
lateral cyclic angle 

1/deg 
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Symbol 

Definition 

Units 

dCsp/dBic 

Change in side force coefficient with 
longitudinal cyclic angle 

1/deg 

dCyM/dAic 

Change in hub yawing moment coeffi- 
cient with lateral cyclic angle 

1/deg 

dCYM/dBic 

Change in hub yawing moment coeffi- 
cient with longitudinal cyclic angle 

1/deg 

dCyfi/ dR 

Change in hub yawing moment coeffi- 
cient with yaw rate 

1/rad/sec 

dCM/dCL 

Change in wing pitching moment with 
lift coefficient 

ND 

da/d8 

Change in fuselage sidewash angle 
with sideslip angle 

ND 

El 

Product of modulus of elasticity and 
moment of inertia 

Ib-in^ 

EIq 

Product of modulus of elasticity auid 
moment of inertia at wing root 

Ib-in^ 

EnFi^5 

Coefficients in the equation for the 
change in normal force coefficient 
with longitudinal cyclic angle 

1/deg 


Coefficients in the equation for the 
change in hub pitching moment coeffi- 
cient with longitudinal cyclic angle 

1/deg 

^SFi-,5 

Coefficients in the equation for the 
change in side force coefficient with 
longitudinal cyclic angle 

1/deg 


Coefficients in the equation for the 
change in hub yawing moment coeffi- 
cient with longitudinal cyclic angle 

1/deg 

Ert# EvT 

Oswald efficiency of horizontal or 
vertical tail 

ND 

F 

Generalized force or force on nacelle 

lb 

FPR 

Lateral-directional SAS function 

— 

FRl 

Lateral-directional SAS function 

— 

F(^ 

Lateral-directional SAS function 
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Symbol 

Definition 
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Units 

* 


F^l 

Lateral-directional SAS function 

j 

t 


F(j>2 

Lateral-directional SAS function 


- 


Fi|;l 

Lateral-directional SAS function 

— 



F<|)2 

Lateral-directional SAS function 

— 



Fa 

Aerodynamic force on nacelle 

lb 

-r 


Fgzn 

Landing gear oleo strut vertical force 

lb 



Fsn 

Landing gear oleo strut lateral force 

lb 



Fx 

Longitudinal generalized force 

lb 

* 



Lateral generalized force 

lb 



Fz 

Vertical generalized force 

lb 

\ 


^yn 

Landing gear oleo strut longitudinal 
force 

lb 



25 

4-1 

Multiplier on rotor normal force 

ND 

* 




Multiplier on rotor power 

ND 



fpM 

Multiplier on rotor hub pitching 
moment 

ND 



fQ 

Multiplier on rotor torque 

ND 



fsF 

Multiplier on rotor side force 

ND 




Multiplier on rotor thrust 

ND 

• 


fyM 

Multiplier on rotor hub yawing moment 

ND 



G 

Generalized moment 

ft-lb 



GEF 

Ground effect factor 

ND 



Ggi 

Governor gain 

deg/sec/rad/ 

sec 



Gg2 

Governor gain 

deg/sec/rad/ 

sec 



Gg3 

Governor gain 

deg/sec/deg 

> 
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Symbol 

Definition 

Units 

L. 

Gp 

Lateral directional SAS gain 

in/rad/sec 

i 

Gprl 

Lareral directional SAS gain 

in/rad/sec 


Gp6s 

Lateral directional SAS gain 

in/in 


Gq 

Longitudinal SAS gain 

deg/rad/sec 

\ 

i 


Lateral directional SAS gain 

in/rad/sec 


Gr2 

Lateral directional SAS gain 

in/rad/sec 

< 

4*-' 


Lateral directional SAS gain 

in/rad/sec 


=6P 

Lateral directional SAS gain 

in/rad 

<1 « 

^er 

Lateral directional SAS gain 

in/rad 

\ 

Gg6r 

Lateral directional SAS gain 

in/in 


®6B1 

Longitudinal SAS gain 

deg/in 



Longitudinal SAS gain 

deg/in 


®6TH 

Governor throttle gain 

deg/in 


Ge 

Longitudinal SAS gain 

deg/rad/sec 



Lateral directional SAS gain 

in/rad/sec 



Lateral directional SAS gain 

in/in 


G,j;6 x: 

Lateral directional SAS gain 

in/in 


g 

Gravitational constant 

f t/sec^ 


H 

Height 

ft 

' 

HP 

Horsepower 

— 

L . 

*^W'FUEL 

Horizontal distance between wing 
mass element center of gravity and 
fuel center of gravity 

ft 

1 

t 

1 

H ’ 

W'NF 

Horizontal distance between wing 
mass element center of gravity and 
fixed nacelle center of gravity 

ft 

1 


Horizontal distance between wing 
mass element center of gravity and 
fixed nacelle center of gravity 

ft 
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Symbol 


Definition 


h ■ 


Height or angular momentum 



hp 


hp 


hf 

hw 


Angular momentum of nacelle c±>out 
aircraft center of gravity 

Distance from wing pivot plane to 
fuselage mass element center of 
gravity 

Height of pivot above wing chord line 
or angular momentum of nacelle about 
the pivot 

Landing gear oleo strut deflection 
during ground contact 

Distance from wing pivot plane to 
wing mass element center of gravity 


ho 


Angular momentum of an element of 
mass about its own center of gravity 


hj^ Wing vertical bending deflection 

h/D Rotor hub height to rotor diameter 

ratio 


he 




^XXq 



Distance from aircraft center of 
gravity to bottom of right main gear 
following a positive pitch rotation 

Distance from aircraft center of 
gravity to bottom of right main gear 
following a positive roll 

Mass moment of inertia 

Vehicle mass roll moment of inertia 
about center of gravity 

Mass roll moment of inertia of air- 
craft components about their own 
center of gravity 

Mass roll moment of inertia of fuse- 
lage mass element about its center 
of gravity 
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Units 

ft or Ib-ft- 
sec 

i 

Ib-ft-sec 

ft 

ft 

ft 

ft 

Ib-f t-sec 

ft 

ND 

ft 

ft 

slug-f t^ 
siug-ft- 

slug-ft‘ 

slug-ft- 
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Symbol 

( 

■XX 




^XX 


Lyy 


■•yyo 


I (F) 

Ayy 


(W) 

l-yy 


I ' 

yy 


■X2 




t(F) 

"X2 


.(W) 

■X2 


I ' 
X2 


■22 


■zz, 


Definition 


Mass roll moment of inertia of wing 
mass element about its center of 
gravity 

Mass roll moment of inertia of the 
tilting portion of each nacelle about 
its center of gravity 

Vehicle mass pitch moment of inertia 
about center of gravity 

Mass pitch moment of inertia of air- 
craft components about their centers 
of gravity 

Mass pitch moment of inertia of fuse- 
lage mass element about its center of 
gravity 

Mass pitch moment of inertia of wing 
mass element about its center of 
gravity 

Mass pitch moment of inertia of the 
tilting portion of each nacelle about 
its center of gravity 

Vehicle mass product of inertia sdsout 
center of gravity 

Mass product of inertia of aircraft 
con^onents about their own centers of 
gravity 

Mass product of inertia of fuselage 
mass element about its center of 
gravity 

Mass product of inertia of wing mass 
element about its center of gravity 

Mass product of inertia of the tilting 
portion of each nacelle about its 
center of gravity 

Vehicle mass yaw moment of inertia 
about center of gravity 

Mass yaw moment of inertia of aircraft 
components about their own centers of 
gravity 


Units 

slug-ft^ 

slug-ft^ 

slug-ft^ 

slug-ft^ 

slug-ft^ 

slug-ft^ 

slug-ft^ 

slug-ft^ 

slug-ft^ 

slug-ft^ 

slug-ft^ 

slug-ft^ 

slug-f t^ 
slug-ft^ 
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Syirbol Definition Units 


(F) 

^zz 


Mass yaw moment of inertia of fuselage 
mass element about its center of 
gravity 

slug-ft^ 

l(W) 

zz 


Mass yaw moment of inertia of wing 
mass element about its center of 
gravity 

slug-ft^ 

^zz 


Mass yaw moment of inertia of the 
tilting portion of each nacelle about 
its center of gravity 

slug-ft^ 

X 


Incidence angle 

deg or rad 



Unit vector in i direction 


*^xx 


Dummy inertia, lzz“^yy 

slug-ft^ 

Jyy 


Diiinmy inertia, Ixx'^zz 

slug-ft^ 

^zz 


Dummy inertia, lyy“Ixx 

slug-ft^ 

A 

i 


Unit vector in j direction 

— 



Wing slipstream correction factor 

ND 

KPI ^ 
T 

KP4 

T 

Coefficients of curve fit equation 
for wing download as a function of 
rotor height/diameter ratio 

ND 

Kmi 

T 

Km4 

T 

Coefficients of curve fit equation 
for wing pitching moment as a func- 
tion of rotor height /diameter ratio 

ND 

K.; 


Multiplier on slipstream rolling 
moment coefficient 

ND 

Kn 


Miltiplier on slipstream yawing 
moment coefficient 

ND 

Kst 

n 


Landing gear spring constants 

Ib/ft 

%1-^KwiO 

Coefficients for wing bending 
equations 

— 



Multiplier on longitudinal cyclic 
pitch available from longitudinal 
stick 

in/in 
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Symbol 

Definition 

Units 

K6e 

Ratio between longitudinal stick 
motion and elevator deflection 

deg/ in 

K6r 

Multiplier on longitudinal cyclic 
pitch available from pedal displace- 
ment 

in/in 

^'^RUD 

Ratio between pedal and rudder 
deflection 

deg/in 


Multiplier on longitudinal cyclic 
pitch and differential collective 
available from lateral stick 

in/in 

^6 ' s 

Lateral cyclic pitch/degree of 
differential collective pitch 

deg/ deg 

Ke 

Wing stiffness in torsion 

ft-lb/rad 

Kq 

Coefficient of fuselage drag coeffi- 
cient equation to account for drag 
due to sideslip 

l/rad^ 

Kl 

Coefficient in fuselage drag coeffi- 
cient equation 

l/rad^ 

K2 

Coefficient in fuselage drag coeffi- 
cient equation 

1/rad 


Coefficient in fuselage lift coeffi- 
cient equation 

1/rad 

K4 

Coefficient in fuselage lift coeffi- 
cient equation 

1/rad^ 

K5 

Coefficient in fuselage pitching 
moment coefficient equation 

1/rad 

H 

Coefficient in fuselage pitching 
moment coefficient equation 

l/rad2 

K? 

Coefficient in fuselage side force 
coefficient equation 

1/rad 

KS 

Coefficient in fuselage side force 
coefficient equation 

1/rad 

K9 

Coefficient in fuselage yawing moment 
coefficient equation 

1/rad 
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Symbol Definition Units 


^10 

Coefficient in fuselage yawing moment 
coefficient equation 

l/rad2 

K20 

Wing/body interference effects on C . 

>. 3 

1/rad 

K21 

Wing planform effects on C q 

C P 

1/rad 

K22 

Wing planform and lift effects on C^g 

1/rad 

K30 

Coefficient in nacelle drag coeffi- 
cient equation 

1/rad 

K31 

Coefficient in nacelle drag coeffi- 
cient equation 

l/rad^ 

K32 

Coefficient in nacelle lift coeffi- 
cient equation 

1/rad 

^34 

Coefficient in nacelle pitching momenu 
coefficient equation 

1/rad 

^35 

Coefficient in nacelle pitching moment 
coefficient equation 

1/rad- 

K36 

Coefficient in nacelle side force 
coefficient equation 

1/rad 

K37 

Coefficient in nacelle side force 
coefficient equation 

1/rad- 

K38 

Coefficient in nacelle yawing moment 
coefficient equation 

1/rad 

K39 

Coefficient in nacelle yawing moment 
coefficient equation 

1/rad- 

K40 

Coefficient in nacelle yawing moment 
coefficient equation 

1/rad 

K4I 

Coefficient in nacelle yawing moment 
coefficient equation 

l/rad2 

K42 

Coefficient in fuselage lift coeffi- 
cient equation 

ND 

i 

Unit vector in k direction 

— 

Ls 

Nacelle shaft length from pivot to 
spinner 

ft 
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Symbol 

'V 

I 




^•AC 


^0 

^PA 


M 

m 


M/? 

mf 

"*N 

N 

NF 


Nn 


N 


1 IMD 


N* 


Definition 


Rolling moment 

Distance from nacelle pivot to 
nacelle center of gravity 

Horizontal distance from nacelle pivot 
to aircraft component center of 
gravity positive - positive forward 
from pivot 

Horizontal distance from horizontal 
tail quarter chord tc wing aero- 
dynamic center 

Horizontal distance from pivot to 
center of gravity of fuselage mass 
element 

Wing root lift/foot 

Horizontal distance from pivot to 
center of gravity of pilot-*' station - 
positive forward from pivot 

Horizontal distance from pivot to wing 
mass element center of gravity 

Pitching moment 

Pitching moment, or aircraft mass 


Pitching moment/rotor thrust 

Mass of fuselage structure 

Mass of one nacelle 

Mass of wing 

Yawing moment 

Rotor normal force 

Engine gas generator speed 

Engine gas generator indicator 

Engine gas generator speed at sea 
level standard, static conditions 


Units 
ft- lb 
ft 

ft 

ft 

ft 

Ib/ft 

ft 


ft- lb 

ft- lb or 
slugs 

ft-ib/lb 

slugs 

slugs 

slugs 

ft- lb 

lb 

rev/min 

rev/min 
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Symbol 

Definition 

Units 

•^19 IND 

Referred engine gas generator speed 
indicator 

— 

N:i 

Engine power turbine speed 

rev/min 


Engine power turbine speed at sea 
level standard static conditions 

rev/min 

P 

Body axes roll rate 

rad/sec 

PC 

.lorizontal distance from wing leading 
edge to pivot location 

ft 

pN 

Nacelle axes rC 1 rate 

rad/sec 

pR 

Nacelle wind axes roll rate 

rad/sec 

P 

Body axes roll rate 

rad/sec 

Q 

Body axes pitch rate or rotor 
torque 

rad/sec or 
lb- ft 

^IND 

Torque indicator 

ND 

Qmax 

Maximum engine torque available 

lb- ft 

qN 

Nacelle axes pitch rate 

rad/sec 

qR 

Nacelle wind axes pitch rate 

rad/sec 

Q* 

Engine torque at sea level standard 
static condition 

Ib-f t 

q 

Body axes pitch rate or frees tream 
dynamic pressure 

rad/sec or 
lb/ft2 

•3s 

Dynamic pressure of rotor slipstream 

Ib/ft- 

R 

Body ■'xes yaw rate or rotor resultant 
force c' rocor radius 

rad/sec or 
lb or ft 

RHP 

Rotor horsepower 

— 

rN 

Nacelle axes yaw rate 

rad ''-’ec 

rR 

Nacelle wind axes yaw rate 

rad/sec 

r 

Body axes yaw rate 

rad/sec 

r 

Radius vector 
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Symbol 

Definition 

Units 


Landing gear tire radius 

ft 

s 

Surface area 

ft2 

SF 

Rotor side force 

lb 

SHP 

Shaft horsepower 

— 

SHP* 

Engine shaft horsepower at sea level 
standard static conditions 

— 

T 

Rotor thrust 

lb 

TEA 

Engine referred turbine inlet 
temperature 

deg 

('^ige/'^ogeJ 

Ratio of the rotor thrust in ground 
effect to the thrust out of ground 
effect 

■■ 

Ti->-T3 

Coefficients of curve fit equations 
for rotor/rotor interference 

ND 

t 

Time 

sec 

u 

Body axes longitudinal component of 
velocity at aircraft center of gravity 
or rotor hub, wing, horizontal and 
vertical tail velocities referred to 
rotor shaft and local surface chord 
axes, respectively 

ft/sec 

U' 

Body axes longitudinal component of 
velocity at rotor hub and wing aero- 
dyneunic center 

ft/sec 

UPA 

Body axes longitudinal component of 
velocity at pilot's station 

ft/sec 

V 

Total velocity 

ft/sec 

Vt 

Rotor tip speed 

ft/sec 

V 

Resultant flow through rotor disc 

ft/sec 

V* 

Non-dimensional rotor forward velocity 

ND 

V 

Total velocity vector 

ft/sec 
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Symbol 

V 


V 


'^PA 

V* 

W.L. 

W 

WDTIND 

W 


W 


WPA 


De finition 

Body axes lateral component of 
velocity at aircraft center of 
gravity or rotor hub wing, horizontal 
and vertical tail velocities referred 
to rotor shaft and local surface 
chord axes, respectively 

Body axes lateral compone’'it of 
velocity at rotor hub and wing 
aerodynamic center 

Rotor induced velocity 

Body axes lateral component of velo- 
city at pilot's station 
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Units 

ft/sec 


.ft/sec 

ft/sec 

ft/sec 


Non-dimensional rotor induced velocity ND 


Fuselage water line position 
Weight of aircraft conponents 
Fuel flew indicator 

Body axes vertical component of velo- 
city at aircraft center of gravity 
or rotor, hub, wing, horizontal and 
vertical tail velocities referred to 
rotor shaft and local surface chord 
axes, respectively 

Body axes vertical component of 
velocity at rotor hub and wing aero- 
dynaunic center 

Body axes vertical component of velo- 
city at pilot's station 


X 


aero 


in. 

lb 

ft/sec 


^subscript Longitudinal distance, measured 

positive forward from nacelle pivot 
along body axes 

Longitudinal force, measured positive lb 
forward along body axes 


ft/sec 


ft/sec 


ft 


Total longitudinal aerodynamic force 
at center of gravity measured positive 
forward along body axes 


lb 


XXX 
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Symool 


vsprscript 

subscript 


^North 

^subscript 


subscript 


aero 


ysprscript 

subscript 

^East 

^subscript 


'subscript 


'aero 


ySprscript 

^subscript 


^down 


Definition Units 

Longitudinal force, measured positive lb 
forward along body axes 


Longitudinal ground track velocity 

Lateral distance, measured positive 
along right wing along body axes 


ft/sec 


^w' fuel 


Lateral force, measured positive along lb 
right wing in body axes 

Total lateral aerodynamic force at lb 

center of gravity measured positive 
along right wing in body axes 

Lateral force, measured positive lb 

along right wing in body axes 

Lateral ground track velocity ft/sec 

Vertical distance, measured positive ft 

down nacelle pivot along body axes 

Vertical force, measured positive lb 

down along body axes 

Total vertical aerodynamic force at lb 

center of gravity, measured positive 
down along body axes 

Vertical force, measured positive lb 

down along body axes 

Vertical ground track velocity ft/sec 

Vertical distance from nacelle pivot ft 

to center of gravity of aircraft 
component, positive down from nacelle 
pivot along body axes 

Angle of attack rad 

Angle of sideslip rad 

Vertical distance between wing fuel ft 

center of gravity and wing mass 
element center of gravity 
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Symbol 

A ' 

w' fuel 


w'w 


6 


6 ' 


c 


6 • 


CR 


6 ' 


F 


I 


'ht 


*5 STEER 


e 

eo 


^iLR 


eiRL 


? 

?wl"*‘?w4 


Definition Units 

Vertical distance between fixed ft 

nacelle center of gravity eind wing 
mass element center of gravity 


Vertical distance between wing center ft 
of gravity and wing mass element 
center of gravity 

Control element (surface or stick) deg or in. 

eingular or linear displacement 


Vertical distance between cargo center ft 
of gravity and fuselage mass element 
center of gravity 

Vertical distance between crew center ft 
of gravity and fuselage mass element 
center of gravity 


Vertical distance between fuselage ft 

center of gravity and fuselage mass 
element center of gravity 

Vertical distance between horizontal ft 

tail center of gravity and fuselage 
mass element center of gravity 

Nose wheel steering angle, positive deg 

right 

Vertical distance between vertical ft 

tail center of gravity and fuselage 
mass element center of gravity 

Wing or rotor downwash angle rad 

Wing downwash angle at zero wing rad 

angle of attack 


Rotor/rotor interference angle, left rad 
rotor on right rotor 

Rotor/ rotor interference angle, right rad 
rotor on left rotor 


Wing on rotor interference 


rad 


Rotor sideslip angle or damping ratio rad or ND 
Wing damping ratio ND 
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Symbol 


Definition 


Units 


fuel 


H ' 

w'NF 




n ' 
c 


n 


CR 


n 


F 


tirt 


nvT 

'^VT 


n?R 

6 



Horizontal distance between wing fuel ft 
center of gravity and wing mass 
element center of gravity 

Horizontal distance between fixed ft 

nacelle center of gravity and wing 
mass element center of gravity 

Horizontal distance between wing ft 

center of gravity and wing mass 
element center of gravity 

Horizontal distance between cargo ft 

center of gravity and fuselage mass 
element center of gravity 

Horizontal distance between crew ft 

center of gravity and fuselage mass 
element center of gravity 

Horizontal distance between fuselage ft 
center of gravity and fuselage mass 
element center of gravity 

Horizontal tail efficiency ND 

Horizontal distance between horizontal lb 

tail center of gravity and fuselage 
mass element center of gravity 

Vertical tail efficiency factor 

Horizontal distance between vertical 
tail center of gravity and fuselage 
mass element center of gravity 

Transmission efficiency 

Aircraft pitch or Euler angle or 
temperature ratio 

Wing twist angle 

Rotor collective pitch angle at 
three-quarter radius 


ND 

ft 

ND 

rad or ND 

rad 

deg 


X Angle between the rotor shaft and a rad 

line drawn through the nacelle center 
of gravity from the pivot 
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Symbol 

Definition 

Units 

u 

Rotor advemce ratio = V/QR 

ND 


Tire sliding coefficient of friction 
when sliding sidewards (for concrete) 

ND 


Tire rolling coefficient of friction 
(for concrete) 

ND 


Coefficient of rolling friction for 
brakes 

ND 

^Rl'^'iR4 

Terms in wing immersed area calcula- 
tion 

— 

P 

Ambient air density 

slug/ ft 

0 

Fuselage sidewash angle 

rad 

"h 

Ambient density ratio 

ND 

T 

Angle between freestream velocity and 
rotor resultant force 

rad 

TD 

Engine response time constant 

sec 

d 

TE 

Engine response time constant 

sec 

■THT 

Horizontal tail effectiveness 

ND 

■^LAS 

Load alleviation system time constant 

sec 

tVT 

Vertical tail effectiveness 

ND 

TP 

Lateral directional SAS time constant 

sec 

Tr 

Lateral directional SAS time constant 

sec 


Lateral directional SAS time constant 

sec 

^4>6s 

Lateral directional SAS time constant 

sec 

T 

4' 

Lateral directional SAS time constant 

sec 


Lateral directional SAS time constant 

sec 

U 

Rotor thrust lesponse time constant 

sec 

T2 

Rotor thrust response time constant 

sec 


Aircraft roll angle or Euler angle 

rad 
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Symbol 

Definition 

Units 

(jip 

Rotor swashplate phase angle 

rad 


Functions in wing vertical bending 
equations 

— 

X 

Rotor wake skew cuigle 

rad 


Aircraft yaw angle or Euler angle 

rad 

Q 

Rotor or engine rotational speed 

rad/sec 

Q 

Angular velocity vector 

rad/sec 

01 

Natural frequency 

rad/sec 

“wl^“w3 

Wing natural frequencies 

rad/sec 


f 

(' 

/ 
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. Subscripts 


A 

AC 

ACT 

AERO 

a 

B 

c 

CG 

CR 

C/4 

DUM 

E 

EFF 

e 

F 

FAC 

FUEL 

FUELcg 

FUS 

F* 

f 

GLAS 

GYRO 

g 

HL 


Available 

Aerodynamic center 
Actuator 

Aerodynamic force 
Aileron 

Longitudinal stick 
Cargo 

Center of gravity 
Crew 

Quarter chord 
Dummy variable 
Engine 
Effective 

Elevator or effective 
Fuselage 

Fuselage aerodynamic center 
Fuel in wing 
Fuel center of gravity 
Fuselage 

Fuselage minus landing gear 
Flap 

Load alleviation system 
Gyroscopic 
Ground or gust 
Left rotor hub 
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Subscripts 

HR Right rotor hub 

HT Horizontal tail 

HTCG Horizontal tail center of gravity 

IGE In ground effect 

i Immersed 1 

L Left wing or rotor 

LAS Load alleviation system 

LE Left engine 

LG Landing gear 

L-L Rotor lead- lag 

LN Left nacelle 

LR Left rotor 

LRH Left rotor hub 

LT Left wing tip 

LW Left wing 

LWq Left wing referred to freestream 

MAX Maximum 

N Nacelle or natural frequency 

NF Fixed portion of nacelle 

NFCG Fixed portion of nacelle center of gravity 

NL Left nacelle 

NR Right nacelle 

NT Tilting portion of nacelle 

n Landing gear index, n*l left gear, n=2 right gear, 

n=3 nose gear 

OGE Out of ground effect 
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Subscripts 


P 

POWER 

PA 

R 

RE 

REQ 

RIGID 

RN 

RR 

RRH 

RT 

RUD 

RW 

RWq 

S 

SP 

STALL 

T 

TH 

VT 

VTCG 

W 

WAC 

WCG 

X 


Power, nacelle pivot, or rotor polar moment of 
inertia 

Power 

Pilot station 

Right wing, rotor or rudder pedal 

Right engine 

Required 

Rigid 

Right nacelle 
Right rotor 
Right rotor hub 
Right wing tip 
Rudder 
Right wing 

Right wing referred to freestream 
Rotor shaft, side, or lateral stick 
Spoiler 
Stall 

Tail, total or wing tip 

Throttle 

Vertical tail 

Vertical tail center of gravity 
Wing 

Wing aerodynamic center 
Wing center of gravity 

Along the longitudinal axis, positive forward 
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Subscripts 

y 


Along the lateral body axis, positive out 
right wing 


z Along the vertical body axis, positive down 

__ Denotes a vector quantity 

Superscripts 


(c) 

(CR) 

F 

F’ 

HT 

(HT) 

IGE 

LW 

N 

NL 

NR 

RW 

T 

VT 

(VT) 

W 

^^'fuel 

(Wf') 

(W'jjf) 

(W'w) 


Refers to cargo or payload weight 
Refers to aircraft crew weight 
Fuselage 

Fuselage less landing gear 
Horizontal tail 

Refers to horizontal tail weight component 
In ground effect 
Left wing 
Nacelle 

Left wing tip at pivot 
Right wing tip at pivot 
Right wing 

Total of horizontal and vertical tail 
Vertical tail 

Refers to vertical tail weight component 
Wing 

Refers to wing fuel weight 
Refers to fuselage weight component 
Refers to weight of fixed portion of nacelle 
Refers to wing weight component 
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Superscripts 

" Denotes an interim calculation or coefficient 

in local wind axes 

' ' ' Denotes an interim calculation 

Denotes average value 

* Denotes interim calculation or calculation in 

freestream wind axes 

' Denotes an interim calculation 

+ Denotes an interim calculation 

A Denotes a unit vector 
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1.0 INTRODUCTION 

The rising costs and diminishing availability of fossil 
fuels, the increasing congestion at major airports, and the 
growing need to reduce noise and air pollution are strong 
reasons for evaluating rotary -wing vehicles for the short-haux 
air travel market. 

The low disc loadings of rotor configurations allow 
vertical or short takeoff and landing for a relatively low in- 
stalled horsepower. This power economy results in improved 
fuel efficiency and reduced air pollution. The capability for 
V/STOL operation greatly reduces runway requirements and pro- 
vides a means to alleviate air traffic congestion at airports . 

The tilt rotor aircraft combines the V/STOL advantages of 
the helicopter with the speed and altitude advantages of fixed 
wing aircraft. In Reference 1 a study was made to define a tilt 
rotor aircraft capable of carrying 100 passengers over a 200 
nautical mile stage length at minimum direct operating cost. The 
configuration emerged as a 4-engined, 33905 kg (74749 lb) air- 
craft with a wing span of 25 meters (82 ft) , a rotor diameter 
of 17.16 meters (56.3 ft) and a cruise speed of 349 KTA3 . This 
represents a very large increase in size for a tilt-rotor 
vehicle compared to current (NASA-Army Bell XV-15) and past 
(XV-3) tilt rotor designs. The question then aris'^s as to the 
flying qualities of such a rotary wing vehicle and the impact 
of operating large tilt rotors in projected terminal area navi- 
gation systems. 

This report presents the development of a mathematical 
model which could be used in a piloted simulation of a large 
tilt rotor aircraft. This is a model of the 1985 Tilt Rotor 
Configuration originally planned for use on the NASA-Ames 
Flight Simulator for Advanced Aircraft. The model could also 
be used in research on advanced controllers for terminal area 
operation. 

The aircraft selected for modeling was the design point 
tilt rotor aircraft described in Reference 1 and detailed in 
Section 2.0. The math model is full force, with all inertial 
and aerodynamic terms included. The forces and moments 
generated by the large, hingeless rotors are calculated from a 
set of equations derived from a regression analysis of full- 
scale test data on a rotor of similar design. Direct calcula- 
tion of the rotor forces and moments in real time is not prac- 
ticable because of the complexity of the equations required to 
represent the flap- lag coupling effects of soft-in-plane hinge- 
less rotors. 

The aerodynamic in^erference effects of the rotor on wings 
and tails, the effect of the wing upwash on the rotor, and the 
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interference of one rotor on the other in edgewise flight/ are 
represented. Turbine engine performance, dynamic response and 
performance limits, both thermodynamic and mechanical, are in- 
cluded. 

The control system elements represented are pilot command, 
three axis stability augmentation and a thrust management/ 
governor system. Control system actuator dyneimics are included 
as first and second order lags. 

The effect of the tilting rotors and nacelles on the air- 
craft center of gravity and inertias are calculated. Forces 
and moments resulting from acceleration of the nacelles during 
rapid tilting maneuvers are included. 

The airframe c.g. and inertia representation permits the 
location, inertia, and c.g. of major components of the aircraft 
to be entered. All lengths, overall c.g. and inertias of the 
aircraft are then calculated. 

Wing and nacelle aeroelastic effects are treated on a 
quasi-static basis, i.e., coupling is through the aerodynamic 
terms only. 

The mathematical model is presented in detail in Appendix 
E and derivations of important equations are also included in 
the body of the report and the appendices. 
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2.0 DESCRIPTION OF AIRCRAFT 


The 1985 commercial 100 passenger tilt rotor aircraft is shown 
in Figures 2.1 and 2.2. Table 2.1 lists the major dimensions 
and characteristics of the aircraft. 

The 1985 tilt rotor has a takeoff gross weight of 33,905 kilo- 
grams (74,749 pounds) . The rotors are three-bladed and are of 
hingeless fiberglass composite construction. The rotor diam- 
eter IS 17.16 meters (56.3 feet) and the solidity ratio is 
0.089. In hover and low-speed flight, cyclic pitch control is 
applied to the rotor to provide control power and trim. These 
rotors are highly twisted (36 degrees) by comparison with 
helicopter blades to provide efficient operation at high advance 
ratio as well as in hover. 

The rotors and forward rotor transmission tilt; however, the 
engines, mounted outboard of the tilt package, remain stationary. 
This arrangement does not require the engines to be requalified 
for vertical operation and reduces the inertia of the tilt 
package . 

The aircraft has four engines, two on each wing tip. The rotors 
and engines are connected by means of a cross-shaft which provides 
the torque transmission across the aircraft in event of engine 
failure. The location of the engines outboard of the tilt package 
provides easy access to the engine bays for maintenance or engine 
removal. 

The span of the aircraft is 25 meters (82 feet) measured from 
outboard of one: nacelle to outboard of the other. Tne wing is 
p^.^uight and upt^pered with a NACA 634 221 section with a wing 
setting angle of 2° relative to the fuselage. The wing aspect 
ratio IS 7.14. 

The wing has full-span 30-percent-chord plain flaperons used as 
both flaps and ailerons. A leading edge umbrella flap is pro- 
vided which opens for hover and low-speed helicopter-type flight 
to alleviate the rotor download on the wing. This device is 
also used to ensure that wing unstalling at end of transition 
occurs simultaneously on both wings. 

The empennage T-tail configuration was selected tc reduce the 
impact of rotor downwash on the horizontal stabilizer in transi- 
tion flight. The horizontal tail volume ratio is 1.47, and the 
vertical tail volume ratio is 0.159. 

The tricycle landing gear configuration provides good ground 
handling characteristics is retractable. The undercarriage 
provides an overturning angle of 27®. 
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TABLE 2.1 


1985 COMMERCIAL TILT 

WING ; 

AREA (REFERENCE) 

ASPECT RATIO 
SPAN (BETWEEN ROTOR CJ, ) 
TAPER RATIO 
CHORDS : 

ROOT 

TIP 

MEAN AERODYNAMIC 
SWEEP BACK 
DIHEDRAL 
INCIDENCE: 

ROOT 

TIP 

AIRFOIL SECTION; 

ROOT 

TIP 

FUSELAGE ; 

LENGTH 

DEPTH 

WIDTH 

WETTED 


ROTOR - DIMENSIONAL DATA 


747.5 FT^ 
7.15 
73.1 FT 
1.00 

10.23 FT 
10.23 FT 
10.23 FT 
0 DEGREES 
0 DEGREES 

2.0 DEGREES 
2.0 DEGREES 

NACA 634221 
(MODIFIED) 

NACA 634221 
(MODIFIED) 

92.5 FT 

11.5 FT 
14.75 FT 
3563 FT- 


2-4 



D238-10002-1 


1985 COMMERCIAL TILT ROTOR - DIMENSIONAL DATA 

(Continued) 


NACELLES ; 

ENGINE : 

LENGTH 

DEPTH 

WIDTH 

WETTED AREA (PER NACELLE) 
TILTING: 

LENGTH 

DEPTH 

WIDTH 

WETTED AREA (PER NACELLE) 
HORIZONTAL TAIL : 

AREA (EXPOSED) 

AREA (REFERENCE) 

SPAN 

ASPECT RATIO 
TAPER RATIO 

DISTANCE (c/4)^g to (c/4)g,p 
CHORDS ; 

ROOT 

TIP 

MEAN AERODYNAMIC 
SWEEPBACK AT 0 PERCENT CHORD 
DIHEDRAL 


8.58 FT 

5.5 FT 

3.0 FT 
136. FT2 

15.5 FT 
5.08 FT 
3.42 FT 
122 FT2 

211.5 FT2 

227.5 Ft2 
35 FT 
5.38 
0.625 

50.75 FT 

8.0 FT 

5.0 FT 
6.62 FT 

9.75 DEGREES 
0 DEGREES 
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1985 COMMERCIAL TILT ROTOR - 

(Continued) 

INCIDENCE: 

ROOT 

TIP 

AIRFOIL SECTION: 

ROOT 

TIP 

\^RTICAL TAIL : 

AREA (EXPOSED, EXCLUDES DORSAL) 
AREA (REFERENCE) 

SPAN (REFERENCE) 

ASPECT RATIO 
TAPER RATIO 

DISTANCE (c/4)w to (c/4)vT 

CHORDS : 

ROOT 

TIP 

MEAN AERODYNAMIC 
SWEEPBACK AT 0 PERCENT CHORD 
AIRFOIL SECTION 
CONTROL SURFACES : 

FLAPERON: 

AREA (AFT OF HINGE) 

SPAN (LENGTH EACH SIDE) 

CHORD (% OF WING CHORD) 
SWEEPBACK OF HINGELINE 


D238-10002-1 

DIMENSIONAL DATA 


( 

0 DEGREES 
0 DEGREES 

NACA 64A010 (MODIFIED) 

NACA 64A010 (MODIFIED) 

210 FT2 
278 FT- 
20,67 FT 
1,536 
0.482 

40.5 FT 

13,33 FT 
8.83 FT 
14.1 FT 
39 DEGREES 

i 

NACA 64A008 (MODIFIED) 

150.4 FT- 

24.5 FT 
30 

0 DEGREES 
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1985 COMMERCIAL TILT ROTOR - DIMENSIONAL DATA 

(Continued) 


SPOILERS : 

AREA 

SPAN (LENGTH EACH SIDE) 

CHORD (% OF WING CHORD) 

POSITION OF LE (% OF WING CHORD) 
SWEEPBACK OF HINGE LINE 
LEADING EDGE UMBRELLA: 

AREA (PLAN) 

CHORD (% OF WING CHORD) 

SWEEPBACK OF HINGE LINE 
ELEVATORS : 

AREA (AFT OF HINGE LINE) 

CHORD (% HORIZONTAL TAIL CHORD) 
SWEEPBACK OF HINGE LINE 
RUDDER: 

AREA (AFT OF HINGE) 

CHORD (% VERTICAL TAIL CHORD) 
SWEEPBACK OF HINGE LINE 
ROTORS : 

NUMBER OF BLADES (PER ROTOR) 
DIAMETER 

BLADE AREA (PER BLADE) 

GEOMETRIC DISC AREA (TOTAL) 
SOLIDITY 


63.2 FT2 

24.5 FT 
12.65 
66.0 

0 DEGREES 

102.8 Ft2 

18.6 

0 DEGREES 

62.1 Ft2 
33 

3.0 DEGREES 

49.3 FT- 
25 

26.7 DEGREES 


56.3 FT 
70.1 FT^ 
4978 FT' 
.089 
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1985 COMMERCIAL TILT ROTOR - DIMENSIONAL DAT A 

(Continued) ~ 

AIRFOIL SECTION: 


ROOT 

VR-7 

70% RADIUS 

CD 

1 

> 

89% RADIUS 

< 

1 

TIP 

VR-9 
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3.0 EQUATIONS OF MOTION 


This section presents the derivation of the airrrame equations 
of motion and the siirplili cations that were made in order to 
obtain the final equations as presented in Appendix E. The 
treatment accounts for all six rigid-body degrees -of- freedom 
including the effects of the tilting nacelles and rotors. The 
principal features of the derivation are: 

o Assumption of X-Z plane of symmetry 

o The basic equations are derived about the 
instantaneous center of gravity of the air- 
craft since the center of gravity is strongly 
dependent on nacelle incidence. 

o Rotor and engine gyroscopic terms are included. 

o The wing elastic degrees of freedom do not 
couple inertially. The coupling occurs only 
through the aerodynamic terms. 

o Wing aeroelastic effects are not included in 
the center of gravity calculations. 

3 . 1 AXES SrS'iUM 

A set of right-handed orthogonal axes OXYZ is placed at the 
center of mass of tlie aircraft and is fixed in the aircraft 
such that G',{ lies in the lateral plane of symmetry and is 
positive iorward parallel to the fuselage water line zero. The 
remaining axes are placed as shown in Figure 3.1. 

The orientation of the aircraft is defined with respect to a 
set of earth-fixed axes C X'Y'Z'. With the axes OXYZ initially 
parallel to C X'Y'Z', the aircraft is yawed to the right about 
0 through an angle ip, then pitched up about OZ through the 
angle 0 , and finally rolled right about OX through the angle i . 

If V and ^ are the aircraft velocity and angular velocity 
vectors relative to the earth-fixed axes, the projections of 
these vectors on the moving axes are U, V, and W for the 
components along OX, OY, and OZ, and P, Q, and R for the angu- 
lar velocity components. 

Thus , 

V = Ui + V2 + Wk (3.1) 

£ = Pi + Qi + Rk (3.2) 
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where the unit vectors ]S 1^® along OX/ OY , and OZ . 

3.2 AIRCRAFT GROUND TRACK 

The components of V relative to the earth-fixed axes are ob- 
tained in terms of U, V, W and ip , 9, ^ as, (See Reference 2 ) , 

U cos 0 cos ^ + V (sin $ sin 9 cos - cos $ sin ’^) 

+ W tcos iji sin 9 cos ij; + sin 0 sin i|/) 

U cos 6 sin + V(sin (|i sin 9 sin i|) + cos ^ cos i|() (3.3) 

+ W (cos $ sin 9 sin - sin ^ cos 41) 

-U sin 9 + Vsin cos 0 + Wcos (fr cos 0 

"ntegration of these equations gives the aircraft ground track. 

A further relationship may be obtained between the rate of 

change of the Euler amgles (4, 0, i>) and the components of the 

angular velocity in the moving axes system, viz, 

^ “ (Rcos ()) + Qsin 4) sac 9 

6 « Q cos t “ Rsin $ (3.4) 

i ■» P + i sin 9 

3.3 FORCE EQUATION 

The total external force, F, acting at the aircraft center of 
mass is given by 

- “ df ^ 

where m is the mass of the aircraft and ^ is the rate of 

change of V witii respect to the moving reference frame OXYZ, 

X • e . 

• 5 =»Ui+vi + Wk (3.6) 

If F has components Fj^, Fy, and F^ along the respective axes 
then 


dt 


dt “ 


dt “ 
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P = 


fx i 


+ F 


yi 


+ F 


2 ^ 


III 


{Oi 


. A 

Vi + 


. A 

Wk 


P 

U 


A A 

j k 
Q R 
V W 


thus 


] 


Fjj * m (U + QW - RV) 

Fy = m (V + RU - PW) 

Fj = m (W + PV - QU) 

The forces Fj^, Fy and F^ are given by 


“ -'^AERO “ ^9 sin 6 
“ ^AERO ^9 sin ^ cos 0 
^2 “ ^AERO '*■ ™9 cos cos 0 


Where Xaero» etc.» are the components of the total aerodynamic 
force acting at the aircraft center of mass. 


Substituting equations (3.5) in equations (3.7), the follow- 
ing equations are obtained for the aircraft accelerations, 


U 

V 

w 


'AERO 

m 

^AERO 

m 

^AERO 

m 


- g 
+ 9 
+ 9 


sin 0- QW + RV 
cos 0 sin i(> - RU 
cos 0 cos + QU 


+ PW 
- PV 


( 3 . 9 ) 


3.4 MOMENT EQUATION 

The derivation of the equations for the total moment acting 
ad}out the aircraft center of mass is complicated by the fart 
that the center of mass changes position due to the tilting 
nacelles. Thus the centers of gravity of the principal air- 
craft con^onent masses of the wings (mw) , fuselage (including 
tails) (mf) , and nacelles (m^) , move with respect to the ref- 
erence axes OXYZ placed at the instantaneous overall center of 
gravity of the aircraft. The equation of motion for such a 
mass element will first be obtained and the total moment found 
by adding the contributions of all the elements.. 
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3.5 EQUATION OF MOTION FOR A MASS ELEMENT 


With reference to Figure (3.1) O'xyz is a right-handed set of 
axes placed at the center of gravity of the representative mass. 
The axes are parallel to the set OXYZ. The mass, m, rotates 
about its own center of gravity with angular velocity ^ which, 
in general, differs from ^ the angular velocity of the aircraft. 
If r is the radius vector from 0 to O' then the velocity of the 
center of mass of the element is 


V = X r 


(3.10) 


The angular momentum of this mass about 0 is 
h = m (r X V) + ho 


(3.11) 


where ho is the angular momentum of m about its own center of 
mass and is given by 

ho = r u) 


where 


^xx “^xy ~^xz 
“^yx ^yy ”^yz 
”^ix ”^zy ^2z 


(3.13) 


^xx» etc., are the moments and products of inertia of the 
mass about O'xyz. 

The total moment, G, about the aircraft center of mass is given 

Using equations (3.10), (3.11), and (3.12) in (3.14), the moment 
becomes 


° *(lf * ^ t) ^ It (il * ii ^ ij' 

+ m n X ^r X £ X rjj + ^ 


4-(l^) 

0 t — 


(3. IS) 
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which reduces to 


G = 2m?. /r. 

rl- 

5^r 

mr X — ^ + ni 

in , 

— (r.r) 

-mr / r . 

6n 1 

-v~ 

6t; 

- ot- 

5t 

-l- 

It J 

-2m l£ 

(£.r) 

-m(r..^<) (f;xr) 

^ 5 oj 

+ I ^ + 

n x/i 


6t 



6 1 

- 1 

-/ 


(3.16) 


The only masses that possess angular velocities different from 
that of the aircraft are the nacelles., which are free to pitch 
about O' with angular rate i = . Thus ^ may be written 

generally as dt 

^ + (Q + in) (3.17) 

Now, with r = Xi + y£ + Z^, where X, Y, and Z are the in- 
stantaneous coordinates of the individual mass center relative 
to the aircraft mass center, the \arious terms of equation 
(3.16) are, in component form. 


r. ^ = XX + YY + ZZ 
— 0 1 

rx = (YZ-ZY)i - (XZ-ZX)£ + (XY-YX)k 

(r.r) = (X^ + y 2 + z2) (Pi + Qi + Rk) 

6t - j. - 

6 n “ ♦ • 

r. ^ = XP + YQ + ZR 
6 1 

n.r = XP + YQ + ZR 

(r.fi) (£xr) = (XP+YQ+XR) (CZ-RY) l- (PZ-R:<) £+ (PY-XQ) k"] 

i = (Ixx^-^xz^^)i + ^yy (Q+^'n)! + (^zz^-^xz^)! 

6 

J2x(fw) = (QR Iza'Q^^xz-'^^OIyy-R^N^yy)! 

- (PR Izz“P^l xz -PR + R^Ixz)i 

+ (QR Ijj2tPQIyy t Pljqlyy “ PQ ^XX ) 


where, in the last two terms, the products of inertia Ixy and 
Jyz are zero fiom symmetry considerations. 

Substituting, the above relations into equation (3.16) and 
noting that Y- and Y are always zero (no lateral motion of thv 
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individual masses) the following expressions are obtained for 
the components of the moment G = + AM^ + ANk: 


AL = P[I^^ + m(Y2 + Z2)J - (R+PQ)II ^2 + m XZJ 


+ RQ[l 22 -Iyy+m(y 2 -z 2 )] + m YZ(r2-q 2) - 
+ m (YZ-2XYR - 2XZR + 2ZZP - XY (Q - PR) ) 

AM = Q tlyy+m(x2+z2)j - (R2-p2)[i^^ + mXZ] (3.19) 

+ PR Uxx"^z2 + m(z2-x2)] + lyylj^ 

+ra [XZ-XZ + 2Q(ZZ+XX)-XY (P+RQ) +YZ (PQ-R) ] 

AN = R (l22+m(-X^+Y2) ] - (P-RQ) XZ] (3.20) 

+ PQ [lyy-lxx + m(X2-Y2)] + lyyp 

+ m (2.XXR - YX - 2XZP - 2YZQ - YZ (Q+PR) +XY (Q2-p2 ) ] 
Summing the rolling moment equation: 


L - IxX ^”^XZ CR+PQ) + CIzz~^YY^^Q 

-2Q(iyi;-i^j,)Y^-2R(XjjjZyR + Xj,i.Z,.j,) + 2P(Z^K2„i, * 

Z!iLZNL)-(a-fi<HXNR-Xi,.i,) Yy]+ 2n, Zj (PZf - 
RX£ ) + 2m^Z^ (PZ^ “ RXt^) “R ^yy (^NL ^NR^ 


(3.21) 


where IxX/ ^XZ' ^ZZ» ®^d IyY are the inertias of the aircraft 
about its center of gravity, and the subscripts f, w, NL and 
NR stand for fuselage, wing, left nacelle and right nacelle. 

The remaining symbols are defined in the List of Symbols. 
Similar expressions are obtained for the pitching moment and 
yawing moment. In the interests of brevity the remainder of 
the discussion will be lindted to equation (3.21). 

Evaluation of the terms of the rolling moment equation indicate 
that this equation may be simplified considerably without a 
significant change in accuracy. For example, terms containing 
(Xnr-Xjjl) JAay be dropped because Xnr is normally identical to 
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• 

%L' i*®* the nacelles are raised or lowered together at. the 
same rate. Equation (3.21) may thus be written 

^“^XX^“^XZ t mN^N (Xnr-Znl) (3.22) 

where the last term has been retained in consideration of the 
high differential nacelle accelerations encountered during hover 
maneuvers . 

From the relationships presented in Appendix C the last term of 
Equation (3.22) may be rewritten as 

-imjgYjj cos (i^R-X) + i^L (3.2 3) 

-\nr sin (iNR-^) - ^nl 
which may be approximated to 

-imjjYjj ~ ^NL (3.24) 

since the nacelle rates appear as squared terms. 


0 

Similar treatment of the pitching moment and yawing moment 
equations results in the following final form of the moment 
equations . 


^AERO “ ~^Xi ^ ^^ZZ~^YY^^^ 

-imj^Yj^ [ijjR cos Ci^R"^) “ ^nL CiNL”^)^ 


^AERO “ ~ + (^XX“^ZZ)P^ 

^NR (^yy© '*■ 

^ ; N 

^NL {^yyQ 

^AERO “ ^2Z ^"^XZ ■*“ ^^YY"^XX^^® 


( 3 . 25 ) 


Z^ sin 
Zl sin 


(i>;L-N)3] 


+ imjjYj^ [Ij.jR sin (ij^R-X) - Ij,l sin 


where the moments L;s^ero» MaeRO» N^ero represent the sum of 
the aerodynamic moments and rotor/engine gyroscopic moments 


about the aircraft center of mass. i 


N 

yyo 


is the nacelle pitch 


inertia referred to the nacelle- fixed axes system described in 
Appendix C. Equations for the aircraft inertias are also 
presented in that Appendix. 


I 
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3.6 EQUATIONS OF MOTION FOR NACELLES 

The equation of motion for a nacelle is required in order to 
obtain the moment exerted by the nacelle on the wing tip at 
the pivot. This moment is then used in the equations for wing 
twist. 


The angular momentum of a nacelle about its pivot point is 
given by 

hp = X "^n3^ hoN 

(3.26) 


= (rxV) + - nin£p x V 

where r is the radius vector from aircraft c.g. to nacelle 

c. g. 

V is the velocity of the nacelle c.g. 

hojj is the angular momentum of the nacelle about its 
own c.g. 

m^g is the nacelle mass 

and £p is the radius vector from aircraft c.g. to nacelj.e 
pivot 

The term nij^ (£xV) + hQ^^ is the angular momentum of the nacelle 
about the aircraft c.g. (= • 

i.e. ^ - mfj(rp X V) 

The moment about the pivot is 


Go = - mn i- (rp X V) = G^g - ;.G 


-P dt 


dt 


dt 


(3.27) 


Since the quantity has already been obtained (equations 

(3.18), (3.19), and (3.20)), only the remaining term needs to 
be evaluated. 


II 

B 


mjj, _n2. 

X V 

+ £p X 

6V 

Jt 


/5r 
X _= 

+ nxr 1 


X 

J_ / 


£^xr 

1 it 

V6t 

- -/ 



6t \ 

, 

4 


■ 


\1 





n X 

rc X 

_= + n 

X r 






— P 

Ut 







(3.28) 
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Expansion of these terms results in the following expression 
rfir_ or / dr_\ /cr_ \ 6^r oQ 


( 6r„ dr / dr„\ /cr^ \ 6‘r on 

-zP. X _z + n /r -r|_i^ .n) + rp x — ~ + _z: (r.r^ 

dt dt ~\- 6t } “l*:t ~l ~ dt2 dt “ 

- r (r . fl .rj - 2 ^ (r . n) 

- \-P dt) \ dt -^7 dt -P ~ (3. 

Ep " (£p‘E) 

We require only the i component of this vector in order to ob- 
tain the nacelle pivot pitching moment. 

The components of the vectors tp, r and £ are 

Ep “ ^pi. ^n 1 ^pli “ " ^Cci. “ ^CG— 

A A A 

E. ~ ^N— ^N— 

A A A 

n = Pi + Qi+Rk . ^ 

Noting that the j components of -=£» _= are zero (since Y« is 

^ 6 t 6 t 

a constant ) , the above expression yields 


AH - + ZqqXj^ + ZjjXqc + PQ 


- RQ XnYn) 


(3.30) 


Combining this equation with Equation (3.19) and using the 
transformations given in Appendix C, the final equation for the 
right-hand nacelle pivot actuator pitching moment becomes, after 
some simplification, 

• r ^ ■'I ^ / ”'N\r* 

%R=-->-NR 2(ij,-^-,\) 

+ (R^-p2 ) sin (ij^f^- \ ) Cos (ii 4 R“-\)j -(R^-P^) ^zZq ^NR 

- lyy^ Q + i ^[xAEROsin(ii;R-X) + Zaero COS (i-jR-X)j 

- imjjYjj [(R-PQ) sin (iNR-X) - (P + RQ) cos (i^R-X)] 

4. M 

^%ERO 

where M^r^^rq includes the moment resulting from nacelle aero- 
dynamic loads and the rotor gyroscopic moments. The terms 
XaeRO snd Zaero are, respectively, the cotal aircraft aerody- 
namic X and Z forces. 
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The corresponding equation for the left nacelle actuator moment 
is obtained by substituting -Yn=Yn and changing the R subscript 
to L. 


3.7 D ETERMINATION OF ROTOR GYROSCOPIC MOMENTS 

The gyroscopic moments are most readily obtained as follows. 

A set of axes 0"x'y'z' is tak'^n at the rotor hub (rotor c.g.) 
parallel to the nacelle-fixed set of axes OxQyQZ^. Associated 
with each axis are the corresponding unit vectors i' i' and 
k ' . The angular velocity of the rotor v;ith respect to these 
axes is the vector 


~ ^pi * (3.32) 

where is the rotor rotational speed. 

The angular momentum of the rotor with respect to its c.g. is 




where ^ is the inertia matrix 


■Rxi 




^Rr • 


the off-diagonal terms being zero since the axes 0"x'y'z' are 
principal axes of inertia of the rotor and hub. 


In component form the angular momentum of the rotor is 
llo ** ^Ry I 


(3.34) 


With respect to the inertial axes OYXZ, the components of 
are 


ho = ^R^rCos ij-i - IpHRsia i^^k 
The hub moment is therefore given by 




dhp, 3 h„ 

— =5. = + D X ho 

6t 6t 


where = Pi + + R^ 


(3.35) 

(3.36) 

(3.37) 


Substitution of equations (3.35) and (3.37) into equation (3.36) 
results in the following equations for the rotor gyroscopic 
moments . 


✓ 
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(3.38) 

(3.39) 

(3.40) 

The above terms appear in the Computer Representation (Appendix 
E) as additions to the rotor aerodynamic forces and moments. 


- cos ij^ - (1 n+Q) sin ij^ 

» IrPJJr sin ijq + IrKOj^ cos ijj 


'•gyro 


“-IrOr sin ijj - (i^+Q) cos ijj 


i 

i 

f 

I 
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4.0 AIRFRAME AERODYNAMICS 

This section presents the mathematical equations and represen- 
tations of the aerodynamic data for the aircraft without rotors. 
The contribution of the rotors is described in Section 5 . The 
overall airframe aerodynamics are obtained from the following 
components : 

(a) Fuse.lage 

(b) Wings 

(c) Horizontal Tail 

(d) Vertical Tail 

(e) Nacelles 


The data and equations for each of the aerodynamic components 
are discussed below, together with the substantiating methods. 

The aerodynamic data are presented in local wind axes. Reso- 
lution to aircraft body axes is accomplished as described in 
the mathematical model (Appendix E) . Where required, the 
equations have been written so as to be applicable over the 
entire range of angle of attack +180 degrees. 

4.1 FUSELAGE 

The aerodynamic lift, drag, and pitching moment coefficients 
of the fuselage were estimated using the methods of Reference 3 . 
The forces and moments are referred to the point on the fuse- 
lage corresponding to the wing quarter chord position. This 
reference point was selected in order to minimize the number of 
force and moment transfer equations in the mathematical model. 
Wing-to-body carryover effects have been included in fuselage 
loads . 

The equations for the fuselage forces and moments are: 

Lift: CLp ** K^2 K3SinapCosap + K4SinapCosap [ 

SincipCosop I 

Drag: ^Dp * ^DOt?^^ I ^F I '*’^2 ^1 

£ 

I SinapCosap] + 

Side Force: Gyp = K7 SinSpCos3p + Ks SinSpCo? ' r | S: .i3pCos 2p | 
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Pitching Moment: » C» + K 5 SinopCosop + Kg SinapCosop] 


Yawing Moment: 
Rolling Moment: 


SinopCosapI + 

Qjp ~ ^Nqp KgSinSpCosSp + Kio^inSpCosSp | SinSpCos 6 p’| 


C-, =0 


where ep = Tan-l(2),c,^ = 


etc. 


'FUS 


Bp = Tan 


+ 




CMp = 


Mp 


i =''rvs 


etc. 


and ACq^q, aCmlg» landing gear contributions to fuse- 

lage drag and pitching moment coefficients, when the landing 
gear is extended. 

The fuselage forces and moments are then resolved into body 
axes at the aircraft C.G. 


4 . 2 NACELLES 


The forces and momf.nts acting on the nacelles were estimated 
using the cross-flow methods of Reference 4. For convenience 
the resulting forces and moments are referred to the rotor hub, 
so that they may be added directly to the rotor forces and 
moments . The following equations are for the forces cind moments 
on two nacelles: 






K 


32 



K 


36 



sin ajj cos 
+ K 3 o[ajjl + K31 
+ K 34 sin ajj cos ajj + K 35 
sin cos Bjj + sin 
+ K32sinSNCOs6j4 + K3QSinEjj 


sin ctj^cos aj^jsin ajjcos 
cos Bj^Isin cos 3^1 
cos 3j^-lsin£j^- cos | 


‘Ni 


V.. 


s 
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The nacelle forces and moments in nacelle axes are: 

“n = In cos sin - c sin 6^, cos a„]l 


AY' 

N 

N 


AH' 

N 


AN' 

N 


- % cos sin B.J 1 


N 


“•N 


N‘ 


- qN cos ON - cos Bn sin on -Cyj, sin SNsin on]^ 

= qN Sj^ i’wI-'B;;' C.^J^ sin Sjj cos o^, - sin 




'N 


N"2 


’n ®W^W^*^Nj, “ 


. 1 

N ” '^w' 


C sin B„ cos o„]i 
% N N 2 


4.3 HORIZONTAL TAIL 


Aerodynamics of the horizontal tail were obtained using the 

combination with test data. The 
horizontal tail includes a plain elevator. 

attack of the horizontal tail, including inter- 
ference effects, for zero elevator deflection, is 


“HT 


= Tan"^f!lHTl- , + i 
L’^HT 


HT 


where e is the total downwash at the tail due to wing, rotor 
and ground effects and Iht is the tail incidence angle. 

The effect o* elevator deflection on the effective tail angle 
of attack is introduced through the elevator effectiveness 
parameter, tht/ which is a function of the elevator and hori- 

areas. Thus the affective horizontal tail angle of 
attacic IS ^ 


**®HT * ^HT'^e 

where Bq is the elevator deflection. 

The tail downwash angle, e , depends on wing angle of attack and 
slipstream deflection. At a given rotor angle of 
attack, the slipstream deflection is a function of rotor thrust 
coefficient, Cts» where the coefficient is based on the slip- 
dynamic pressure. Figure 4.1 presents data on downwash 

® tunnel model 

’• seen, the downwash at low values of 

thrust coefficient is the same as the value of the power-off 
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ROTOR THRUST COUFFICIEIiT , Ct 

•^s 

Figure^ . 1 Variation of Horizontal Tail Downwash Angie with Thrust Coefficient 


wing downwa.jh (Cts " 0) • Above values of Cts in the neighbor- 
hood of Cts=. 5 the downwash increases with increasing thrust 
coefficient. The values in the increasing portion of t vs Cts 
were found to correspond approximately to the slipstream deflec 
tion angle cp. Therefore, the approach adopted in the mathe- 
matical model was to test if the rotor slipstream downwash 
(Fp) exceeded the wing downwash and, if so, to use the computed 
slipstream downwash value as the tail downwash angle. Other- 
wise the wing downwash value was used. 

Thus if 

fp > «o + If- (“w - ‘AC 

then t - ‘P > 
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Otherwise 
c « 




(a 


w “ ^AC 


W 


) 


(1-GEF) 

/1-m2 


In these expressions is the wing downwash angle at zero wing 

angle of attack, ^ is the downwash derivative, is the 

dcx 

distance from the wing to the tail aerodynamic cei ' ers, and 
Hac —7 is the fconiliar downwash lag term. In general, the 

^ Hr 

quantities cq and ^ depend on the average of the left and 

right flaperon deflections. The effect of differential deflec- 
tion of aileron/spoiler in producing an asymmetrical downwash 
field at the horizontal tail was not included because of the 
small contribution this makes to total aircraft rolling moment. 


The term (1-GEF) in the above equations is the ground effect 
factor. This quantity was obtained from Reference 2 and is a 
function of the wing span and height of the horizontal tail 
above the ground (Appendix E, Page E-42) . This factor, when 
multiplied by the downwash which would be found out of ground 
effect, yields the downwash in ground effect. Ground effect 
is discussed in more detail in Section 8 . 

The lift and drag forces acting on the horizontal tail are re- 
quired over the complete range of angle of attack -180® to 
+180®, since the tilt rotor can fly backwards. The following 
sketch shows the schematic variation of lift and drag coeffi- 
cients over this range plotted as a function of the effective 
horizontail tail angle of attack, ^eHT* 
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The angle of attack for Clktw. is denoted by anT+ and is the 
value of the effective angle or attack at the stall less 2 
degrees i.e. 

«HT+ = ^"“HTstalL 

Similarly the angle of attack for stall at negative angles of 
attack is 

^HT. - -'“KTstaLL 'HT‘e 

The slope of the lift curve within this range of positive and 
negative angles of attack is given by 

Cl “ Cl 

0 ttHTl ^ I 

✓i-M^ 

where ag/a is the ra tio c f tail lift-curve slopes in and out of 
ground effect, and ► 1-M“ is the Prandtl-Glauer ; correction 
factor for the effect of Mach number on lift-curve slope. 
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Within this region on the lift curve the value of lift coeffi- 
cient is given by Clht ~ “eKT ^he corresponding drag 
coefficient by 

c »c + 

°HT °Ojj.p ©HT 

After stall angle of attack is passed the lift is assumed to 
fall linearly to zero at Og = +90°. 


In these regions the lift is given by 

” Cl., ‘t 

where the appropriate signs are taken depending on the sign of 
'‘eHT* 


The corresponding drag is obtained by assuming a linear varia- 
tion of drag from the value of ^ value of Cq = 1.1 

(flat plate normal to stream) at “ 90°. Thus 



STALL 



“HT+ 


^^STALL 



^HTstaLL 

aarht ®HT 


and 


“stall ®ht 


*HT ^ l-“Cn 

+ “ht 


STALI 


(+90-aHT+) 


If the effective angle of attack of the hcritontal tail exceeds 
+90° the tail will point trailing-edge first into the relati’/e 
wind. Under this condition early stalling is precipitated be- 
cause of the sharp "leading edge" and blunt "trailing edge". 

In order to represent this, it was assumed that the attainable 
Clmax tail under these conditions is half that occurring 

in normal flight. 


Thus if 90° < Cl- < (180 - i Suit. ) 

®HT ~ 2 - 


1 ^ 


or (-180 + j oht+> i “ejjT ^ 
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then 


'Lht 


= .5Ct cx 


“HT. 


(ae^T - 90 °) 

(90“- 


or 


Cl 


HT 




(-90 + 1 ) 

2 


The corresponding drag coefficients are; 

for 90" < a- < (180 - ^ a ) ; 

GrT — ^ ill* 


'HT 


STALL 


^HT 


= 0.5 Cl ^ht 
a "" 

Cf + C„ 

= HTgijia^LL °HT 


STALL 


which gives C = C 

Cht “HTstaLL 


itARRT ®IIT 

j. too 0.5 ctR'T' “180") ) 

+ ®KT ^HTSTALL 


(0.5 aHT_“90‘ 


and for (-180 + ^ < -90"; 


Cl 


HT 


STALL 


= 0.5 Cl oht+ 


C£ 


+ Cr 


■'HT 


HT 


STALL 


STALL 


uAR Rij ®HT 


which gives 


= Cr 


®KT 


(a^ +180"-.5aHTJ (1.1 -Cd^. 


HT *^HT 


‘HT 


ktstall 


STALL 




In the range (180-.5aHT.) 1 <^eHT - 
unstalled 


Cl = Cl («e - 130°) 

^HT -^a ®HT 


CDrt = + Lrt 


n ARhTSht 
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and similarly for the range -180* ^ “©HT C“180 + .5 
^Lht " ^ ^ 

=Oht ' 

itARjjT®HT 

The cibove equations define the variation of tail lift and drag 
over the entire range of angle of attack. The tail pitching 
moment is not computed since it makes only a small contribu- 
tion to the total aircraft pitching moment. 

4.4 VERTICAL TAIL 

The aerodynamic forces and moments acting on the vertical tail 
were estimated using the methods of Reference 3. The angle of 
attack of the vertical tail is given by 


- Tan‘1 ’'VT 




(S) 


where u^ij, Vyiji, and Wy.p are the components of the velocity at 
the vertical tail aerodynamic center as given in Appendix C. 

The term 6p I is the sidewash correction for the presence 
' dS 

of the fuselage. 

As in the treatment of the horizontal tail, the effect of 
rudder deflection is obtained using a rudder effectiveness 
parameter xym. Thus the effective angle of attack of the 
vertical tail when the rudder is deflected is 

^evT ^ ^VT ■^VT '^RUD 

The treatment of the vertical tail aerodynamics through the 
complete angle of attack range -180° to +180° then follov" che 
same lines as that for the horizontal tail aerodynamics p 
viously described. 

The vertical tail forces and moments in body axes are then ob- 
tained from; 


VT - r 

^AERO COS ^ 




4-9 


I 


! 


D238-10002-1 


VT 


r 

AERO “ 

^^VT^VT 


VT 



'AERO “ 

^^VT^VT 

^“Cjj^^cos ( ) sin -C. 


'VT 


>VT 




VT 


VT 


VT 


“aero 

VT 

Naero 

VT 

LAERO 


^AERO ^^CG"^VT^ ^AERO ^^VT"^CG^ 
VT 

"'^AERO (Xcg-Xy^.) 

VT 

- yaero (zvt-zcg) 


4.5 WING AERODYNAMICS 

The treatment of the wing aerodynamics is the most complex of 
all the components. Because wing flexibility must be repre- 
sented, each wing panel required a separate treatment. The 
approach adopted for simulation purposes was first to obtain 
the aerodynamic forces and moments on the complete wing con- 
sidered as rigid and uninfluenced by slipstream interference 
effects. With this data as a basis, the effects of elastic 
deflection were introduced as an increment in the effective 
angle of attack of each wing panel and the rotor slipstream 
interference was then calculated. This approach is described 
in detail below. 


4.5.1 BASIC WING AERODYNAMICS 

The basic wing lift, drag and pitching moment coefficients for 
the wing in the presence of the fuselage rotors-off, were cal- 
culated using the methods of Reference 3 . This data is appli- 
cable to low speed flight. Corrections for Mach n umber effects 
are introduced through the Prandtl-Glauert factor >'1-M“ . Be- 
yond stall angle of attack, the lift, drag, and pitching 
moment curves are extended linearly to +90® angle of attack in 
order to provide a representation of wing behavior at low 
transition speeds when wing angles of attack approach 90® . The 
data was calculated for the complete range of flaperon settings. 


The complete wing basic lift, drag, and pitching moment data 
also applies to each individual wing panel provided the data 
is obtained at the appropriate panel angle of attack. This 
approximation is acceptable if the angles of attack of each 
wing panel are not substantially different. This condition is 
normally fulfilled. 
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In addition to the above data, the effects of spoiler deflec- 
tion on pcuiel lift, drag, and pitching moment are required. 

These were estimated using the data of Reference 3. As can be 
seen from the equations presented in Appendix E the spoiler 
effectiveness is strongly dependent upon flaperon deflection, 
a result of the spoilers being slot-lip spoilers. 

4.5.2 ROTOR SLIPSTREAM INTERFERENCE 

Before the basic wing aerodynamic data can be utilized in the 

calculation of the wing forces, the effects of the rotor slip- 

stream must be calculated. The calculation procedure presented / 

here has been developed and used at Boeing for some years , and 

gives acceptable agreement with wind tunnel test data on a wide 

variety of both tilt rotor and tilt wing configurations. 

The method uses momentum theory to obtain the direction and 
speed of the rotor slipstream in the neighborhood of the wing. 

From this the effective angle of attack of that part of the 
wing that is immersed in the slipstream is calculated. The 
lift, drag, and pitching moment on the wing are then calculated 
for this jjigle of attack as if the entire wing were immersed. 

The area of the wing immersed in the slipstream is now computed 
and, using the ratio of the immersed to total wing area, the 
forces acting on the immersed portion are approximated. 

At the angle of attack of the wing outside the slipstream, the 
wing forces and moments are obtained from the basic wing data 
as if no slipstream effects were present. These forces are 
then scaled by the ratio of unimmersed to total wing areas to 
obtain approximately the forces acting on the unimmersed wing. 

The sum of the approximations to immersed and unimmersed wing 
forces is now formed. This Siam is then multiplied by a correc- 
tion factor to obtain the final forces . 

This correction factor is obtained from a consideration of the 
mass flows associated with the rotor-wing combination. In the 
following outline of the method only one rotor is considered. 

From the following sketch, which shows tne forces acting on 
the nts^r, the inclination of 




0238-10002-1 

the resultant force on the rotor to the freestream direction 
is given by 

tr = Or + Tan“l 

The resultant force on the rotor is 
R = /t2 + NP‘ + SF-^ 

where T, NF and SF are the thrust, normal force and sideforce, 
respectively. 

The mass flow through the disc is 
m = p A V 

where A is the disc area and V is obtained from the induced 
velocity triangle at the disc plane. 

V’ « / (Vq + Vj^ cos c)‘ + (•^ sin t)^ 

The resultant force on the rotor is related to the mass flow 
by 

R = 2m = 2pA V v^ 

From these equations the following quartic equation is obtained 
for the induced velocity at the dir.c. 

vi + 2V*v| cos T + v| V| « 1 
where the nondimens ional notations 



have been introduced. 
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This equation is then solved for v* and the direction of tne 
slipstream just behind the rotor disc is calculated from 


c 


P 


Tan 


v*sin T 
v*cos T+V* 


The rotor thrust coefficient Ct is defined as 

s 


^Tg - (q + T)A 
A 

with T = R cos (t - ctj^) 

and q = j PV^ = ^ v 2 r 

then CTg = cos (r-aR) 

cos ( T + V* 

T" 


NOTE: Because the rotor diameter 

to wing chord is large the 
slipstrecun is considered 
to be uncontracted in the 
vicinity of the wing. 


T^e aspect ratio of the slipstream- immersed wing area is given 



where is the immersed area calculated by the method de- 
scribed in Appendix D, and c is the wing chord. 


The lift on the wing/ if the slipstream were absent/ is obtain- 
ed by calculating the effective angle of attack of the wing 


Oq = sin 


-1 


ww 


+ 


/u2+w2 

L w w J 

where w^/ u^ are the velocities at the wing aerodynamic center 
and is the elastic twist at the point. The lift coefficient 
(C*) for this angle of attack is obtained from the aerodynamic 
data for the appropriate f laperon/spoiler deflection. 


Similarly the lift (C£) and drag (Cp) coefficients of the wing 
in the slipstream (assuming wing is completely immersed) are 
obtained fromthe aerodynamic data at the angle of attack 


a 


s 


lift coefficient of the wing with slipstream is 


The total 
therefore 
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= K' 


where 


A 

_ L 


i (C£ cos e - C£ sin e)+ C£ (1-C^^) (1-^) 


q^Sw 


in which qg is the nominal slipstream dynamic pressure, defined 

by qs = q + A * 

The factor is a correction factor to account for the fact 
that the lift-sharing between the immersed and unimmersed por- 
tions of the wing is not simply proportional to the respective 
areas . 

From considerations of the mass flows associated with the 
ving— rotor combination the factor was obtained in the form 

V* + 

=• Ct 


V* 


V* + V* 

where, from wing theory. 

Cl ^ 


^gj = 1 + Cl U_ " L.' 

Ct — ± AR i AR 

a n L J 


The drag and pitching moments for the wing with slipstream are 
obtained similarly and are given by: 


Cr. = K 

Dg A 


CMs “ ^A 


^ s • s * 1 

-i (c;; sin s + CJ cos o + CJ a-c^^) a-i)j 
^ cji t cji a-cj^) (1-^1 ) 


V. 

The rolling moment and yawing moment coefficients for the 
wing are given by: 


■ »20+ ^<21 <l-=Ts> V ‘ 


"LW 




+ 

•^SpoWER 

=ns " ■'22 ^ ^ 

+AC. 


(1-=T< 


(C* - C* ) 

AC Drw i^LW 


SpoWER 
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where the increment in rolling moment due to power is 


I' 


-IC - (1-C ) Cf 


.^‘2 ( S 


RW 


and the increment in yawing moment is 


AC. 


^POWER 


“HP 


Sr„ 


5 J H (W„ 


-|CDc - (l-CTs)Cg iri-ify ' 
I ®LW ^LwJ 1 2 \ S / LW _ 


Figure 4.2 shows a correlation between the wing-in-s Upstream 
method described above and experimental results for the Boeing 
Model 160 tilt rotor aircraft. As may be seen the simple 
treatment gives acceptable predictions of wing forces and 
moments . 
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The mathematical representation of the rotor for the 1985 
tramsport is based upon full scale test data obtained on the 
26 foot rotor for the Boeing Model 222 tilt rotor. The test 
data was curve-fitted and scaled by solidity to yield equations 
suitable for representing the 1985 transport rotor. Where test 
data was not availadale, the rotor performance was calculated 
using Boeing rotor performance computer programs. The use of 
mathematical expressions for the rotor forces and moments re- 
sults in maximum computational efficiency and minimum cycle 
times. This method is preferrable to table look-up procedures. 

5 . 1 Sign Convention 

The sign convention for rotor forces and moments is defined in 
Figures. 1» which shows the rotors under combined pitch 
^“t.L ” sideslip 3. The resultant rotor angle 

of attack is given by 


aj^ = cos“^ (cos ^t.L. 


and the rocor disc "sideslip" angle 


Cjj = Tan 1 


Tan 3 

sin a^.L. 


- 


The resulting rotor forces and moments are defined with respect 
to the plane containing the resultant rotor angle of attack 
e.g. normal force lies in this plane while rotor side force is 
perpendicular to it. 


5 . 2 Isolated Rotor Aerodynamics 

The equations used to represent the isolated rotor aerodynamics 
are presented below. The equations are used to compute the 
rotor wind- axes forces which are then resolved through the 
rotor sideslip angle into nacelle axes and hence transferred 
to aircraft body axes for use in the equations of motion. 


5.2.1 Thrust Vs 675 


The thrust produced by the rotor at any 
obtained from the following equations 


0 



6. 3015u 


flight condition is 
+ 5.5816u" 


- 8 u sin a + 1.115 


( 1 ) 


i 


i 

1 


4 
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and C.J. is given by 

Ct = 0.000679 $ + 0.000015 ^2 

( 2 ) 

+ 0.0022 + 0 . 000211y 2 iji 

5.2.2 Thrust Vs Power 

Once thrust has been established the power coef^'lcient is given 
by 

Cp = 0.00006 + 0.00057 u + 0.000085 + 1.12 

- 0.024075 Ct + vCt (0.53 + 0.456 u -39.937 uC,p 
+ 31.79 Ct) + [0.0115U - 0.03u2 - Ct(3.4u -8^2) 

7T 

- 0.22064ti (Ct + 0.001971) sin a (3) 

+ (0.3082u - 2.18 u2) Ct sin a 

5.2.3 Normal Force 

Normal force is obtained as the sum of three terms 

Cnf = F(u,a,CT) + + ^*^NF Bi (4) 

3Bj_ 

where the cyclic pitch derivatives are functions of ct, u, and 

Cip . 

In performing the analysis the cyclic derivatives were first 
defined as: 


3C^F 

TT — = 0.0000217“^ + 0.0014483u2 - 0.0000734u 

(5) 

- 0.0006U sin 2ct + 0.00425 Ct 


and 


3Cf,jF 

TT — = 0.0000425 - 0.0010492U - 0.0017028u2 

^ ( 6 ) 

+ 0.00l7S92y sin a - 0.0245 Cp 

The following expressions may be used to calculate normal force 
with zero cyclic pitch. 
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For 0 ^ V ^0.6 

Cnf = CNFi = 0.068u^ sin 2a + [0.133695uCt 
+ se.iiip (1 -u) ]k 
where K = sin a for a > 20® 


(7) 


and K * sin a(10-0.45a°) for 0 £ a ^ 20 
For 0.6 < y 

Cnf = (CnFi) ( l-0.8(v-0.6) ) (8) 

5.2.4 Side Force 

Side force is defined in a similar manner to normal force 
CsF = F(u, Ct, a) + ^CSF Ai + ^CsF 3^^ (9) 

3A]_ 


where the cyclic derivatives are given by: 


SCgF 

3Ai 


=-0.0000328 + 0.0003119y + 0.0013178y2 
+ 0.0189 Cij - 0.001342 y sin a 


( 10 ) 


and 

3Csf 

= 0.00001683 + 0.0011208y^ - 0.0000568y ( 11 ) 

+ 0.00328 Cij. - 0.00052438 y siu 2a 

The side force at zero cyclic is given by the following equa- 
tions : 


C 5 F = 0.00430 y sin a - 0.0028827 y (^rad^“ 



+ 0.012 

y sin a Crp (90-’ 

y°) + 2«19u^ sin a Cip 

where 

y® = tan’^ 

\ y-yj cos a I 


[ Vi^ sin a J 


and 


+ Ct2)1/2 . ^2) 



( 12 ) 

(13) 

(14) 
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5.2.5 Hub Pitching Moment 


Pitching moment is computed in the same manner as normal force 
and side force. 

CpM = F(a, M, RPM, C^) + (15) 

3A]_ 

where the cyclic pitch derivatives are functions of a, , RPM 
and Cm. 


3CpM 


= 0.0001620 + 0.00036652|i 

- 0.00056151y^ - 0.00000591ti (RPM-298) 
+ 0.0002826P sin 2a + 0.0015 Crj. 


= -0.00008609 36 + 0.000056^- 


+O.OOOJ385u^ -0.0019 Ct 
-0.00000551U (RPM-298) 


+0. 00048791m sin a 


CpM = 0.009950 M sin a -0.010960 m“ sin a 

+ 0.0028126 M sin 2» - O. 0057743 u sin a 


I RPMl 
I 29 8 J 


+ (1.802 u sin a - 7.56 (u sin a) 

5.2.6 Hub Yawing Moment 

The yawing moment d ivatives due to cyclic pitch are similar 
to the pitching moment derivativer- and are given by 


= -0.000086093 + 0.0000612 
0.0003385 M- -0.0019 C^. 
-0.00000551 u (RPM-298) 


+0.0003 u sin a 
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and 


3Cym 

3bY“ = -0.0001620 -0.00086652 ;a 

+0.00056151 + 0.00000591 y (RPM-298) 

-0.001 Ct - 0.0003638 p sin 2a 


( 20 ) 


The yaw moment at zero cyclic pitch is given by the following 
equations 

For 0 p ^0.37 

Cyf/j — (0.018369 p —0.0007) P sin a —1,2 p^ Cm sin a 
+ 0.00631 -0.002604P -0.004877 ||^ -ijl 1 "1 

and for p > 0.37 

Cym = (0.01916 - 0.15321 (u -0.5435)=) sin a 
- 1.2 Crp sin a 

5.2.7 Pitching Moment due to Pit h R...te 

dCpjij 


srn a 
( 21 ) 

v22) 


1000 


dQ 


- 1.5 + p 


0 < p < . 2 


= 0.25 + 7.26 p 


.2 < p < .39 


= 4.1681 -2.79 p p > .39 

5.2.8 Yawing Moment due to Yaw Rate 

dCYM dCpfi 

dR ^ ~ dQ 

5 . 3 Ro tor/Rotor Interference 

A procedure for calculatj.ng rotor-on-rotor interference effects 
is included in the mathematical model. Rotor-on-rotor inter- 
ference arises during sidev/ard flight at low airspeeds with 
the rotors up and, to a lesser extent, during slipped flight in 
the transition configurations. The basis for the method is as 
follows . 
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The above sketch depicts the tilt rotor aircraft flying side- 
wards at low speeci. The wake of the upwind rotor interferes 
with the inflow to the downwind rotor producing a change in 
this rotor's forces and moments. 

Reference 6 presents calculated values of the normal component 
of the induced velocity near a rotor having a triangular disc 
loading, for different wake skew angles, This data is used 

to compute an interference angle at the downwind rotor. The 
interference angle is subtracted from the isolated rotor angle 
of attack and the resulting angle of attack is used in the cal- 
culation of the forces and moments. The rotor/rotor interference 
effect is washed out with nacelle angle and sideslip angle so 
that there is no interference at the high end of transition and 
in cruise. The equations used to calculate interference are 
presented in Appendix E under the rotor/rotor interference section. 

5 . 4 Effect of Wing Upwash on Rotor Performance 

The rotor operates in the upwash field associated with the lift- 
ing wing. Thus the rotor behaves as if it were operating at 
an increased angle of attack. The effective upwash angles were 
calculated using lifting line theory. In the mathematical 
model the upwash angles are input in the form of a table of up- 
wash angles as a function of wing lift coefficient, and nacelle 
incidence angle. 
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Control of the 1985 tilt rotor aircraft is accomplished by uti- 
lization of longitudinal cyclic/ differential longitudinal 
cyclic, collective and differential collective pitch, in con- 
junction with the airplane control surfaces. The airplane 
control surfaces consist of conventional elevator and rudder 
and a flaperon and spoiler arrangement. The primary controls 
in each axis for each regime of flight are shown in Table 6.1. 

• 

The rotor controls provide a major portion of the control capa- 
bility from hover through the low transition speed range, but 
airplane surface controls are operative in all regimes of 
flight, including hover. The rotor controls are phased out 
during transition as nacelle incidence decreases, speed in- 
creases, and the airplane controls become more effective. 


6 . 2 Longitudinal Control 

Longitudinal control in hover is provided by longitudinal cyclic 
pitch. This is phased out through transition as the elevator 
becomes more effective. The elevator provides longitudinal 
control in the cruise mode . 


6 . 3 Lateral/Directional Control 

Roll control in hover is accomplished by differential collec- 
tive (thrust) and yaw control by differential longitudinal 
cyclic (thrust vector tilt) . Differential engine power is pro- 
vided (via the governor) to ensure maintenance of roll control 
in the event of cross-shaft failure and also to minimize cross- 
shaft torque. 

In transition, differential collective and differential cyclic 
per inch of control movement are scheduled as functxons of 
nacelle incidence. Longitudinal and lateral cyclic, elevator 
angle and flap deflection are also scheduled with nacelle angle 
to provide a • :)ntrols- fixed trimmed condition through transi- 
tion. 


In cruise, lateral control is provided by f laperon/spoilers 
and rudder. The flaps are full-span, single-slotted of 30 per- 
cent chord wxth a fixed hinge point 14.6 percent below the wing 
chord line. The flaps act as flaperons for roll control and 
deflect downward to a maximum of 20 degrees from the nominal 
flap setting. Maximum incremental lift from the flaps is at- 
tained at approximately 35 degrees deflection and the maximum 
rolling moment occurs at the same time, so the flaperon de- 
flection for roll control is limited to a maximum total flap 
deflection of 35 degrees. If, for example, the flaps are 
symmetrically deflected 30 degrees, only 5 degrees additional 
deflection is utilized for roll control. Full-span spoil ors 
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of 12.7 percent chord are located >rward of the flaps and 
hinged to the rear spar. The spoilers are "slot-lipped", i.e., 
they open up the slot forward of the flap with the flaps ex- 
tended resulting in a large increase in roll control as com- 
pared to the control power with flaps closed. Maximum deflec- 
tion of the spoilers for roll control is 45 degrees from the 
closed position. 

Maximum spoiler rolling moment coefficient is also attained 
with flaps deflected approximately 35 degrees. Spoiler effec- 
tiveness with the flaps retracted is approximately one-third 
that attainable with the flaps extended. 

The spoiler? and flaps are also used in conjunction with down- 
load alleviation devices, referred to as umbrellas, mounted on 
the leading edge of the wing for download relief in the hover 
and low speed range. The umbrellas are IS. 6 percent chord on 
the upper and lower wing surfaces. Maximum deflections of the 
surfaces for download alleviation are; flaps 70 degrees, 
spoilers 110 degrees from closed, and umbrellas aft-edge-of-the- 
upper surface up to 20 degrees from vertical and aft-edge-of- 
lower-surface down to 10 degrees from vertical. The umbrellas 
and spoilers retract at 50 knots automatically. 

6 . 4 Thrust/Collective Control 


In hover, forward motion of the thrust/collective lever mechan- 
ically commands both -increased collective pitch and increased 
power. The governor provides a fine adjustment to the collec- 
tive pitch to maintain rpm. Over-travel of the pilot's lev.-*r, 
beyond the normal max power position, provides a collective 
pitch landing flare capability. The over-travel is entered by 
going through a "gate", which shuts down the rotor governor and 
leaves the pilot's lever directly connected to collective pitch, 
as in a helicopter collective pitch lever. 

The collective pitch is also scheduled through transition as a 
function of nacelle incidence, minimizing the adjustment need- 
ed from the governor. 

In cruise the mechanical interconnection of the thrust/collec- 
tive lever with collective pitch is phased out completely so 
that a pure power demand system with governed pitch, like a 
conventional fixed wing airolane, is provided. The control 
system block diagrams are shown in Appendix E. 

6 . 5 Control Feel 

Control force gradient variation with dynamic pressure prevents 
excessive sensitivity of control at high speed. The force 
gradients of the primary controls (longitudinal and lateral 
stick, and pedals) are varied linearly with dynamic pressure. 
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The rudder and elevator deflections vary linearly with pilot's 
rudder pedal and longitudinal stick travel. Aileron deflection 
is programmed linearly and spoiler deflection non-linearly with 
lateral stick deflection, to provide near-linear rolling moment 
effectiveness to near cruise speed. As mentioned earlier, 
spoiler deflection is limited at high speed by limiting the 
actuator capacity. The control breakout forces and force 
gradients are shown in Appendix E. 

6 . 6 StsJaility Augmentation Systems 

Stability augmentation systems are provided to enhance aircraft 
flying qualities. The system consists of longitudinal, lateral 
and directional SAS. The longitudinal stability augmentation 
system incorporates a pitch rate feedback and a longitudinal 
stick pickoff. In addition, a pitch attitude signal is in- 
corporated to provide some degree of attitude stabilization 
without the autopilot. (An autopilot is not represented in 
this simulation.) These signals are shaped and put through 
an authority limit. The longitudinal SAS commands longitudinal 
cyclic pitch to provide the required damping in hover and 
transition. It is not required in the cruise mode and is 
phased out at 175 knots. The block diagram of the longitudinal 
SAS in given in Appendix e. 

The lateral stability augmentation system is operative in all 
flight modes- It consists of roll rate feedback for increased 
damping in roll, a roll attitude feedback, to provide roll atti- 
tude stability, and a lateral stick pickoff. In addition, a 
sideslip feedback is incorporated to compensate for dihedral 
effect. A lateral SAS authority limit is incorporated in the 
circuit. The output of the lateral stability augmentation sys- 
tem is input to the control system in terms of equivalent 
lateral stick, since the drive actuator is in series with, and 
commands the same control as, the pilots lateral stick control 
linkage. The lateral 5AS never opposes the pilots' command. 

The block diagram of this system is shown in Appendix E. 

A directional stability augmentation system is pr-vided and 
operates in all flight regimes. The yaw channel consists of 
yaw rate feedback for increased directional damping in hover 
and low speed flight modes, yaw attitude feedback to provide 
yaw attitude stability, and a rudder pedal pickoff for quicken- 
ing. Directional damping provided by the rotors is quite high 
in the higher transition and cruise speed ranges. No addition- 
al yaw rate damping is therefore needed in cruise. A feedback 
is provided to modify the effective yawing moment due to roll 
rate which exists in the basic unaugmented aircraft configura- 
tion in the cruise speed range. A directional SA£- authority 
limit is incorporated. The SAS command is input to the control 
system in terms of equivalent inches of rudder pedal. The 
block diagram for the directional stability augmentation system 
is shown .’n Appendix E» 


6-4 



I 


D238-10002-1 


6 . 7 Thrust Management System 

The thrust and power management system for a tilt rotor air- 
craft must be compatible with both the helicopter and airplane 
configurations. Thrust control for the hover task, rpm control, 
gust response (especially in the cruise flight regime) , and 
effect on aircraft flying qualities must all be considered. 
Classically, helicopters have used collective pitch demand to 
control thrust and fuel governing to control rpm while fixed- 
wing aircraft have used fuel flow demand to control thrust and 
collective pitch governing to control rpm. Each system has its 
advantages. For a tilt rotor aircraft it is desirable from a 
practical viewpoint to have one type of governing for both the 
helicopter and fixed-wing flight regimes. Collective pitch 
governing was chosen for the 1985 tilt rotor for several rea- 
sons : 

o It is more readily adapted to the hover flight 
regime than the fuel governor is to cruise 

o It has better gust response characteristics 

o It is fast acting and has high accuracy 

o Thrust response to pilot control can be easily 
shaped with feed forward loops 

o It has been demonstrated successfully in hover, 
transition and cruise in the CL-84 aircraft 

With collective pitch governing there are two areas in the 
thrust management system to be considered: (1) design of the 

collective pitch governor; and (2) the feed forward loops for 
shaping pilot thrust control. The block diagram for this sys- 
tem is shown in Appendix E. 

The governor was designed to meet the following objectives: 

(1) 0.3 percent steady state error in 2.5 to 3 seconds; (2) 

2 percent rpm overshoot; and (3) satisfactory effect on air- 
craft flying qualities in the all-operational mode (i.e., all 
aircraft components operational and performing as designed) . 

A single governor reference that uses the rpm signal from each 
rotor and averages them satisfies the design criteria. To 
achieve the required accuracy and transient response goals , 
integral as well as propc.tional feedback of rpm is necessary 
in both the hover and cruise regimes. Governor gain is schedul- 
ed with nacelle incidence to maintain a near optimum level of 
governing throughout the flight envelope. Gains are varied 
linearly as the rotor rpm is changed from hover to cruise. The 
second requirement of the governor system is shaping the rotor 
thrust output for a pilot throttle input. Considerations ir. 
determining the proper shaping include: 
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(1) throttle sensitivity 

(2) time constant to reach 63% of steady-state thrust 

(3) allowable thrust overshoot 

Variable pilot's control sensitivity is employed to give the 
optimum sensitivity in the hover power range yet maintain full 
power control within a reasonable throttle throw (8 inches) . 
Shaping of the pilot command with collective quickening is 
used to improve the thrust time constant and thrust response 
transient shaping so that the pilot may perfomn the precision 
hover task with a minimum of difficulty. In the cruise regime, 
shaping of the thrust output is unnecessary and is phased out 
during transition. 

The thrust/collective pitch control system is designed in such 
a manner that, during hover, when the pilot moves his control, 
he commands both a change in engine fuel setting and, mechanic- 
ally, a change in collective setting. The governor then oper- 
ates with a time lag to trim the collective to the value re- 
quired to maintain rpm. The mechanical collective change 
feature is washed out as a function of nacelle incidence so 
that when nacelle incidence is decreased Lo zero, the pilot 
commands only engine fuel. 
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7.0 ENGINE MODEL 


This section describes the representation of engine performance 
and dynamics. The basic engine cycle performance data consists 
of tabulated values of four variables; power, fuel flow, gas 
generator shaft rpm, and power turbine shaft rpm. These para- 
meters are a function of Mach number and turbine inlet tempera- 
ture. All data are in referred, normalized format as shown in 
Table 7.1. Because of the normalized, referred format, all 
data are valid for any ambient conditions. The effects on 
engine performance of operating at non-optimum power turbine 
speed are included in the model. The referred format also 
facilitates the inclusion of engine thermodynamic and mechanical 
limits. Limitations on engine cycle operation may be input in 
any combination of the following; fuel flow, torque, gas gen- 
erator speed, gas generator referred rpm or output shaft speed. 
The flow charts which describe this routine mathematically are 
shown in Appendix E, and the performance data in Appendix F. 

A simplified dynamic model of the Lycoming LTC4V-1 engine was 
formulated for use in the tilt rotor mathematical model. The 
model basically consists of two first-order lags in series with 
variable time constants and gains. The output of the model is 
rate-limited to reflect actual engine performance. This simpli- 
fied model gives satisfactory results for both large and small 
power transients. The block diagram for this system is shown 
as part of the thrust management system block diagram shown in 
Appendix E. 


TABLE 7.1 

ENGINE CYCLE 

DATA FORMAT 

VARIABLE 

SYMBOL 

REFERRED, 
NORMALIZED FORM 

Thrust 

^N 


Power 

SHP 

SHP/ 5 /TsHP* 

Gas Generator rpm 

Ni 

Ni/ v^onj 

Power Turbine rpm 

Nil 

Nfi/ V 

Fuel Flow 

Wf 

W^/3 V^F* 


Wf/o vTsHP* 


Turbine Inlet T T/e i 

Temperature j 

Wher¥l * = Max. Power Setting, Static , Sea Level, j 

Standard Day j 

0 = Ambient Temperature (°R) Divided by 518. 69° R ' 

5 = Ambient Pressure (psia) Divided by 14.696 psial 
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8 . 0 GROutJD EFFECTS 

The effects of operating near the ground on the rotors and air- 
frame are included in this model. The presence of the ground 
on the airframe imposes a boundary condition which inhibits the 
downward flow of air normally associated with the lifting action 
of the wing and tail. The reduced downwash has three main 
effects ; 

o A reduction in the downwash angle at the tail 

o An increase in the wing lift-curve slope 

o An increase in the tail lift-curve slope 

These have been accounted for by the methods given in Reference 
2, Appendix B-7. The data given in the reference for the 
change in wing and tail lift-curve slope has been used directly. 
The equation specified for the change in downwash angle at the 
tail due to ground proximity was modified for convenience. The 
equation as stated is: 

(Ae) _ bl + 4(h-H)^ 

; b| + 4(h+H)- 

where (Ac)g = the change in tail downwash angle due to ground 
proximity 

e = the downwash remote from ground 

h = the height of the tail root quarter-chord 
point above the ground 

H = the height of the wing root quarter-chord 
point above the ground 

b^ = a function of wing lift and wing flap geometry 

For this mathematical model/ the b^ in the above equation was 
taken to be equal to the wing span, b^. This results in a 
small error in the change in horizontal tail downwash. It is, 
however, sufficiently accurate for this simulation. 

Ground effects on the rotor are difficult to predict analytic- 
ally, especially in forward flight. Wind tunnel test data for 
the Boeing Model 160 powered model. Reference 5, was plotted as 
a thrust ratio versus effective rotor height/diameter ratio, 
for two rotor advance ratios. This data, shown in Figure 8.1, 
was curve fitted and linearly interpolated for advance ratio. 

The resulting equation is as follows: - (for the right rotor. 
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The left rotor is identical except for subscripts) 


“^IGE 

TqGE, 


RR 



(.1741 - .6216ppp) 


EFF 

RR 

+ 1.2479 - .8806 y RR] 


(-1 

\d 1 


(1.4779 uj^j^-.4143) 


EFF 

RR 


where 



hRR 

2R[|sin(9 + if^j^) cos (ji | + .0174] 


- - ZpowN - Xqq) sin 6 


+ t (Ls sin iNj^ + Zcg) cos ^ sin cji] cos C 


= Rotor hub height above the ground 
L3 = Distance from the nacelle pivot to the rotor hub 

Xqq = Longitudinal distance from the pivot to the CG 

Zcg ~ Vertical distance from the pivot to the CG 

0 = Aircraft pitch attitude 

= Aircraft roll attitude 
^Nr ~ Right rotor nacelle angle 

Yf^ = Wing semispan 

The equation for the effective rotor height to diameter ratio 
(h/D)gpp was derived by dividing the rotor hub height by 
[sin(0+ifj) cos ^1. This yields the rotor height along the 
shaft. For the cruise condition the hub height is infinite, 
(h/D)EFF is infinite and the augmentation ratio due to ground 
effect is unity. Some special conditions which must be observ- 
ed when using these equations are noted in Figure 8,1. 
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9.0 AIRFRAME REPRESENTATION 


An airframe representation/preprocessor calculation is included 
in the mathematical model that enables the user to input the 
location of major structural elements of the aircraft in terms 
of water line, butt line and station line location. All lengths, 
center of gravity distances and inertias used in the equations 
are then calculated. This feature enables the user to quickly 
change the location of major structural elements to assess their 
impact on vehicle response. 


In the derivation of the basic equations of .notion, the aircraft 
was divided into three principal mass elements. The fuselage 
mass element (mf) , the wing mass element {m ^) , and the tilting 
nacelle mass element (m^) . The components of the three mass 
elements are shown below. 






fuselage mass 
element (mf) 


win'-’ mass element 
( 1 %) 

tilting nacelle mass 
element (mj^) 


r Fuselage and contents 

Horizontal tail and contents 
< Vertical tail and contents 
I Crew and trapped liquids 
Cargo 

r Wing and contents 
Fuel carried in wing 
( Fixed nacelles ai.d/or engines 

I Tilting nacelle (including 
] rotors) 


These three mass elements a'nng with their respective distances 
from the nacelle pivot to the center of each mass element are 
used to compute the aircraft center of gravity distances with 
respect to the nacelle pivot. The equations for these center 
of gravity distances, including the effects of nacelle tilt are: 


V .. I'f'f^f + ^wW 

ACG 

m 

s mfhf + myhw 

m 

The masses and distances 
in the sketch below. 


^|^j[cos + cos (iNR-X)j 

^ sin (i^iL"^' + sin (iNR“ ')| 


used in these equations are defined 
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are indicated by the 
direction of the 
arrows 


The quantities required to compute mg, ’.f, m„, m, ■' , m^, ' , 

hf, hw are available from an aircraft three-view dxaxi/ing and a 
standard mass properties buildup. The quantities i and ■ (de- 
fined in the sketch) are easily obtainable from a draw.ig. The 
mass quantities (m, m^^ mf , rtv) are computed from a me prop>^r- 
tles buiiciap by adding up the components of each mass element 
as described in the previous paragraph. The lengths f, , 
hf and hw are computed by summing thv. weignt moments of the 
components of each mass element about the nacelle pivot. The 
equations for these operations have been derived and are pre- 
sented in Appendix B under the prepri^jessor equations. The 
input data to these equations include the weight of each com- 
ponent, and its location in terms of water line, fuselage 
station line, and butt line. 

When the center of gravity distance of each mass element has 
been determined, the component and total aircraft mass moments 
of inertia can be computed. The equations for the total air- 
craft mass moments of inertia are presented in Aopendix C. 

The moments of inertia of each mass element are computed by 
application of the parallel axis theorem. The moments of in- 
ertia of each component about its own center oi gravity must be 
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; 


known. The parallel axis tlieorem states: 
N 

^xx ^ ^ tlxxoi ^ ^i ^ 

i=l 


N 


^yy " ^^yyoi 


i=l 


N 


Izz ” t^zzoi * ^ 


N 


•xz ~ ^ t^xzoi '*' 

i=l ^ 
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whare N represents the number of component masses. 

These equations have been expanded to compute the moments of 
inertia of each mass element and are shown in Appendix E under 
the preprocessor section. 

Other lengths required for the mathematical model are computed 
in this section. The input data for these computations are in 
terms of the water line, butt line, and fuselage station line 
locations of the elements in question. 
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10.0 AEROELASTIC R E PRESENTATION 

Two aeroelastic degrees of freedom are included in the tilt 
rotor mathematical model. These are first mode wing vertical 
bending and first mode wing torsion. The stability and control 
characteristics of flexible airplanes may be significantly in- 
fluenced by distortions of the structure under transient load- 
ing conditions. When the separation in frequency between the 
elastic degrees of freedom and the rigid body motions is not 
large, then significant aerodynamic and inertial coupling can 
occur between the two. Many of the important effects of elastic 
distortion, however, can be accounted for simply by modifying 
the aerodynamic equations. The assumption is made that the 
changes in aerodyneunic loading take place so slowly that the 
structure is at all times in static equilibrium. This is 
equivalent to assuming that the natural frequencies of vibra- 
tion of the structure are much higher than the frequencies of 
the rigid body motions. Thus a change in load produces a pro- 
portional change in the shape of the airplane, which in turn 
influences the load. This is known as the method of "quasi- 
static” deflections where all the coupling occurs in the aero- 
dynamic equations. 

Since for the 1985 transport, the rigid-body short-period modes 
are separated from the elastic modes by a substantial margin, 
the method of "quasistatic" deflection is used to represent the 
wing bending and torsion modes, with the only coupling in the 
aerodynamic terms (through angle, of attack). The wing twists 
and bends instantcineously when subjected to an applied load. 

The assumptions made in deriving the wing bending and torsion 
relationships are as fellows: 

• No coupling between bending and torsion modes 

• Wings are cantilevered from the fuselage 

• Elliptical loading assumed for the rigid untwisted 
wing 

• Aerodyncimic loads act at the wing quarter chord 

• Wing elastic axis coincident with cross shaft 

• Wing center of mass assumed to lie on the elastic 
axis 

• First wing torsional mode assumed linear from tip 
to root 

In the mathematical model, wing twist at the tip is calculated 
using the following equation: 
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Kf 


't " ^ACT “ <3 


w 


bw 


■'ITIq 


+ q 


2 

Cw 


dCmc/4 

dC^ 


+ ^WAC ^Lg ^ 

Cw I I 6z 



3 tt 


“RIGID 


where: K 0 ^ 

0 t 

%CT 

^E 

^E 

R 

q 

bw 

Citio 


= Wing torsional spring constant 
= Wing twist angle in degrees 
= Nacelle actuator pitching moment 

* Engine inertia 
= Engine speed 

* Body yaw rate 

* Dynamic pressure 

= Wing reference chord 
= Wing reference span 

= Wing zero lift pitching moment coefficient 


dCjrtc/4 

dC^ 

“rigid 


wing section pitching moment slope with 
section lift coefficient 

Wing lift curve slope 

Wing angle of attack without twist 


Assuming a linear mode shape from the wing tip to the root and 
a cantilevered wing (zero twist at root) , the wing twist at the 
aerodynamic center location of the wing is obtained by linear 
interpolation. The wing twist represents the change in angle 
of attack of the wing tip and aerodynamic center and are used 
in the aerodynamic equations. 

Wing vertical bending deflection is also treated on a quasi- 
static basis. The form of the equation used in the mathemati- 
cal model for the wing tip deflection is as follows: 
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^Wi Z^RO %2 ^AERO " *^W3 ^aerO " " %g a^^(. 


where: hj^ 

zW 

AERO 

ZN 

AERO 

lN 

^AERO 


Wing tip deflection 
Wing lift 

Total wing lift 

Nacelle rolling moment 


^WAC 

^W]_^^W5 


Vertical acceleration of the nacelle 

Vertical acceleration of the wing aero- 
dynamic center 

Constants for 1985 transport wing 


The form of the equation for the wing deflection at the aero 
dynamic center is written similarly: 


^^WAC ^^5 Z^j^Q - %g a^^^. 

The symbols represent the same quantities as the tip deflections 
except the quantities to Kw^^g are different from to Kg. 

These equations are derived in Appendix A . Since the wings 
are assumed cantilevered, these equations may be written for 
the left and right sides. The equations as used in the mathe- 
matical model are written in Appendix E. 

The wing tip and aerodynamic center vertical bending velocities 
are computed by dividing the change in vertical bending deflec- 
tion by the simulation time frame. The vertical bending deflec- 
tions and velocities are then added to the velocity components 
at the wing tip and aerodynamic center. These velocity com- 
ponents are then used in the calculation of the aerodynamic 
angle of attack. 

In addition to the aerodynamic coupling via angle of attack, as 
discussed above, the wing tip vertical forces and moments act 
as the driving functions to a set of second order equations 
that are forced at the wing vertical bending frequency. This 
results in giving a pilot a "seat of the pants" feel for the 
vibratory aspects of the wing vertical bending mode. 
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11.0 CONCLUSIONS 


A full-force, 11 degree-of-freedom mathematical model of a 
1985-era Tilt Rotor 100 passenger aircraft has been formulated. 
The mathematical model is complete and includes full-force 
nonlinear aerodyncunics of all components , turbine engine per- 
formance and response, thrust management system and governor, 
pilot's controls and automatic stabilization systems. The 
model is suitable for implementation on the Ames Flight 
Simulator for Advanced Aircraft. 
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j> ,PPENDIX A - TREATMENT OF WING FLEXIBILITY 

As described in Section 10 the large separation which exists 
between the natural frequencies of vibration of the wing struc- 
ture and the aircraft rigid body motions, enables the elastic 
deformations of the wing structure to be calculated on a quasi- 
static basis. 

In the simple treatment presented below, the bending and tor- 
sion modes are considered to be uncoupled. The wing is treated 
as a cantilever with a built-in root end. The wing is free to 
twist about the elastic axis which is assumed to coincide with 
the nacelle pivot line. The center of mass of each chordwise 
strip is also taken to lie on the pivot line. The unloaded 
wing has neither geometric nor aerodynamic twist. 

WING TWIST 


Spanwise twisting of the wing takes place under the action of 
the nacelle aerodynamic and inertial moments, the wing lift 
distribution, and the spanwise distribution of aerodynamic 
pitching moment. The nacelle aerodynamic moments consist of 
rotor hub loads, transferred to the pivot, together with the 
aerodynamic loads on the nacelle itself. Nacelle inertial 
moments include the gyroscopic effects of the rotor drive 
system. 

With reference to Figure A.l, Mn is the moment supplied or 
absorbed by the nacelle tilt actuator. If K 0 is the wing 
stiffness as seen by the wing tip, then 


Mjj — K0 0tji (A-1) 

The total moment about the elastic axis due to wing aerody- 
namics, nacelle loads and engine gyroscopic torque is 
b /2 

T = J m dy + Mn + Mgyj-o (A-2) 

o 

The aerodynamic moment about the elastic axis at any station 
y is given by 

M = M ^^4 + ilx (A-3) 

where i is the section lift and x is the distance from the 
quarter chord to the elastic axis. In terms of the section 
aerodynamic coefficients, 

m(y) = jpV^c^ ^mc /4 + ^A-4) 

o 


A-1 
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The section lift coefficient, C^, is given by 

dC 

Co = k (a - aj 




(A-5) 


where 

e 


= k ao (aR - - «o ®t(y^ ^ 

is the wing root section angle of attack 
P 


is the rotor induced downwash, assumed constant 
spanwise 




is the section zero-lift angle 
is the structural twist at station y 


The factor k /l- is introduced so that, for the 

wing, the lift distribution is elliptical. The value 


untwisted 
value of k is 

obtained from the rigid wing elliptical loading as 


k = - Cl 
IT -^a 


(A-6) 


Thus the equation for Cj;^ becomes, with cirjqid = UR-Sp- Oq , 

(A-7) 

In equation (A-4) we can write, for low angles of attack, 

dCm, 


- A ^ 




Co — — Ct 
^ 7T 

“RIGID ^ 

'1-1 b 1 + 9p 

/I- 1 b 1 


Cmc/4 “ ^ 
and therefore 




dC)l 


]c/± Cj 


m(y) = f pV^c^ ^ + 


dC 


mc/4 ^ X 

c 


(A-8) 


(A-9) 


The equation for the total wing twisting moment, equation 
(A-2) , can now be written as, 

T = Mactuator+ ^YRO + J PV^c^Cn,^ b +1 pV^c' 

+ xl Cjidy (A-10) 

\ dCo cl o 


A-2 


/ 
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Figure A. 1. Wing Geometry for Derivation of FiexibiHty 
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Using equation (A-7) , assuming a linear structural twist from 
root to tip and performing the indicated integrations , the 
equation for total wing twisting moment becomes 


T = K00.I. = M 


actuator 


+ i^a^rTn + 4 pV^bC^Cjn + 


"ab 


6ir 


•lira 


RIGID 


gyro 4 
+ 40„ 


2 \ dCi 


cZi+l 


(A-11) 


The equation for the actuator moment is given in the equations 
of motion. Section 5.0. 


_ 1 


Rearranging, and writing q = qe (l”C«r ) = 5- pV 

s * 

0 T = Mjj+Mgyj.Q+'^s (l“CTg) 'tn +i|+by,C 

L \ dCr I 


CA-12) 


■ §? ‘^sJ^wCw Cl„ ( 1 -ct^) 


a 


<^^m + X 
dCL c 


where CMq / the zero-lift wing section pitching moment coeffi- 
cient, is a function of flap deflection: 


-mo 


= Cl + C2fif + C35f' 


(A-13) 


Knowing the tip value of twist, the twist at any other span- 
wise station is obtained by assuming a linear variation of 
twist from zero at the root to the tip value. 

WING VERTICAL BENDING 

The spanwise bending moraent at any spanwise station y , on the 
wing is the sum of the bending moments due to wing aerodynamic 
lift, wing weight, nacelle lift, nacelle weight and net torque 
on the nacelle. The expressions for each contribution to the 
bending moments are derived below. 

o Bending moment due to wing loading . 

Assuming an elliptical distribution of lift the bending 
moment is given by 

a 

M (yj^) = J i(y) (y-yi)dy (A-14) 

yi 


_ Hob' 


-m If ■ 1^1 


] ; 


1 '' 
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where Iq is the lift per unit length at the wing root. Intro- 
ducing the spanwise variable 6=cos~i |2yj making the required 

substitutions and integrating, the beifeing moment at any point 
y is : 


M (y) = o 
4 


■L (sin 9-0 cos 0)- i sin*0 (A-15) 
2 6 


o Bending due to nacelle net vertical load. 


The net vertical force on nacelle is 


- nWN 


where F® is the aerodynamic force and nWjj is the inertial 
load on the nacelle. The bending moment due to nacelle 
force is 


M^(y) = ^ (1-cos 0) 


(A-16) 


o Bending due to wing weight . 

Assxaming a uniform distribution of wing weight 

w 

M (y^) = -n J w(y) (y-yi)dy 

yi 

and w(y) = 2W/b where W is the weight of one wing panel 

b/2 

.*. M'*'(y^) = J (y-y^)dy (A-17) 

^1 

« nWb 1 

i.e. M'^(y) = - “J” (1-cos 0 - J sin*9) 

o Bending due to nacelle torque! (rolling moment) 


T(y) = constant = T 


(A-18) 


Total bending moment at station y is therefore 
M(y) = M^(y) + M^(y) + M'*'(y) + T 


(A-19) 


Assuming a linear variation of El from root to tip given 
by 


El(y) = EIo Tl-a (^) = 


= El, 


(1-a cos 9), (A-20) 
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the curvature of the wing due to bending is 


M{y) _ d*z _ J-ob* 

(sin 9-9 cos 0)-l/3 sin®0) 

+ 

l-COS 0 

EI(y) dy^ 8 EIq 

1-a cos 9 

2EIq 

1-a cos 0 


nWfjb 1-cos 0 1 _ nW^^b [ 1-cos 9 - j sin^0 

2EIq 1-a cos 0 2 EIq 1-a cos 0 


+ ^ r. 

EI^ I 


[l-a cos 0 


(A-21) 


Double integration of this equation yields the following ex- 
pression for the bending deflection of the wing at any point 
y on the span:- 

, s Lb ^ i. , b ^ j nW^b ^ j. 

z(y) = (|)t + i>2 ~ —2— 

8TTEI„ ^ 8EI„ 8EI« 


nWwb^ 

8EIo 


+ Tb^ 
4 El, 


(A-22) 


where <(>1 = j | J 


(sin 0-9 cos 9)- j sin*0 
1-a cos 9 


dy > dy 


= ^3 Vj [ f 


*4 if 


17 / 

o L o 


1-cos 9 
1-a cos 0 


dy ( dy 


1-cos 0 - j sin^9 
1-a cos 0 

dy 7 


dy ) dy 


j / 1-a cos 9^ 

o L o ) 

and where the wing lift (2 wing panels) L = j ^gb. The function 
through ({i- were obtained numerically and are°presented in 
Figure A. 2. 


Since L = -2 Z 
F^= - Z 
T = - L 


W 

AERO 

N 

AERO 

N 

AERO 
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where 


II 

1 _ '^AERO 1 _ r 

2 w m 2 ^w ^AC 


II 


a„ 

T 



m 

w 

is 

the 

mass of two wing panels 

m 

is 

the 

total aircraft mass 


®WAC 

is 

the 

acceleration of the 

wing 

center 




is 

the 

acceleration of the 

wing 


and since the values of ((>, through (ji^ are constant for any 
given station y on the wifig we can write the final equation 
for wing bending in the form 

^1 “ ^W, ^AERO ^AERO " ^AERO " 

12 3 4 


- K, 


■w. 


‘WAC 


where h, 


K, 


W, 




K, 


W, 


K, 


W, 




b*4>2 

8El“ 

o 

41TEI 

o 

8EIq 

8EI^ 

o 


Thii? is the form given in the computer representation. The 
bending deflection at the aerodynamic center and at the wing 
tip are obtained using the values of (J).-*- de appropriate to 
those stations. ^ ^ 


I 
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Figure A.2. Wing Bending Functions 
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APPENDIX B - DERIVATION OF LANDING GEAR EQUATIONS 


Presented below are the equations for landing gear forces and 
moments arising from ground contact. The derivation accounts 
for brake and friction forces together with a simplified re- 
presentation of the oleo dynamics. Nose wheex steering is not 
included. 

With reference to Figure B-1 the distance from the center of 
gravity to the bottom of the right main wheel following a 
positive pitch rotation is 

h- = X sin 0 - Z cos 0 - r (B-1) 

where X and Z are the coordinates of the hub of the wheel 
relative to the C.G. and r is the tire radius. If the air- 
craft is now rolled right, through the angle <p, the bottom of 
the right gear moves through a distance. 

h^ = sin <f) + (Z+r) (cos j cos 9 (B-2) 

The height of the bottom of the wheel above the ground is 
therefore 


» = "CG + hg - 


(B-3) 


and the oleo deflection during ground contact is given by 

(B-4) 


N = 


”CG ~*~^9 -^(j) 


cos <j> cos 6 


By differentiation of equation B-4 and making small angle 
assumptions regarding the aircraft pitch and roll angles during 
touchdown, the rate of change of oleo strut deflection is ob- 
tained as 


H, 


h„ = 


CG 

cos (1> cos 0 


+ XQ - YP 


(B-5) 


Assuming that the oleo response is that of a second order 
system, the equation of motion for the landing gear is 


^G ^ST ^T °ST ^T 


(B-6) 


where Kg^ and Dg^ are the equivalent spring rates and damping 
for the oleo, and F_ is the force on the landing gear strut. 
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Tire Friction and Side Force 

The friction force acting on each tire during ground cor tact 
is resolved into a force F along the line of intersection of 
the plane of the wheel and^the ground plane, positive forward, 
and a side force F at right angles to F lying in the ground 
plane and positive^to starboard. The friction force F is 
assumed to be proportional to oleo force and the amount of 
braking exerted by the pilot. The side force is proportional 
to the oleo force. 


1/ 


The components of tire friction are: 


''u ' ‘“o * >'1 ®G> V 


F- = p 


_ F V 

s ^GZ 

V 


(B-7) 

(B-8) 


where p©/ Vi Vi are the coefficients for rolling friction, 
brake friction and®sliding friction. Bq is expressed as a 
percentage of full brake pedal deflection. The signs of the 
forward and sidewards velocity are introduced to properly 
orient the tire forces. 

The force and moment contributions of each landing gear to the 
aircraft total forces and moments are, assuming small angles; 


“n = - ^G2 » 

n n 

^GZ 

n n 

AZ = F 9-F 4> + F-_ 

" “n =n 

“n ^ * V> 

n 

AL * AZ„ Y - AY (Z + r + h_) 
n n n n n n t 


AN =-AX Y X AY 
n n n n n 


(B-o) 

(B-10) 

(B-11) 

(B-12) 

(B-13) 

(B-14) 


where n=l, 2 and 3 denote the left main gear, right main gear 
and nose gear, respectively. 


The total contribution of the landing gear forces to the forces 
and moments at the center of gravity of the aircraft are; 


i 

i ^ 
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'lg - I 


'lg - I 

n=l 


'LG = I 
n-1 


^LG = I 
n=l 


‘"lg I “n 


“lG i ''“n 
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APPENDIX C - VELOCITY AND ACCELERATION TRANSFORMATIONS 
AND CENTER OF GRAVITY/INERTIA EQUATIONS 

C . 1 Velocity Transformations 

The calculation of aerodynamic forces on wings, fuselage, 
nacelles, and tail surfaces requires that the angle of attack 
and relative wind velocity at these surfaces be known. These 
velocities are obtained most conveniently in terms of the velo- 
city of the pivot reference point. 

With reference to Figure C.l, the velocity of a general point 
in the aircraft relative to the airplane center of gravity is 


V - U 


+ Q X r 


(C-1) 


where r is the radius vector from the c.g. to the point and n 
is the angular velocity of the aircraft. Thus, expanding 
equation C-1, the velocity of the pivot relative to the c.g. is 



II 

+ QZp 

- YpR 

v' 

• 

- PZ,^ 

- XpR 

p 

P 

p 


w' 

II 

+ PYp 

- QXp 


(C-2) 


where Xp, Yp and Z_ are the distances of the pivot from the 
c.g., measured positively forward, to the right and downwards, 
respectively. If we measure all distances from the pivot loca- 
tion then Xp = ~XgG' “0, Zp = -ZcG and the velocity 

of the pivot relative to inertial space can be written. 


Up=U+u^ =U-XcG- QZcg 


Vp » V f v^ = V + PZcG - XcqR 


(C-3) 


Wp = W + w^ - W + QX(|,q - Z^Q 


where U, V, and W are the con¥>onents of the velocity of the air- 
plane center of gravity. 

The velocity of a point in the aircraft relative to the pivot 


C-1 







r 




« m 

u = 

X 

+ 

QZ 

- YR 


V = 

• 

Y 

+ 

RX 

- PZ 

- - 

w = 

• 

z 

+ 

PY 

- QX 


where X, 

Y 

$ 

and 

Z are 
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(C-4) 


aU’UX jiiv^ L.J. \v-» - / ^ 

ing components are obtained relative to inertial space. The 
indicated distances are measured relative to the pivot. 


Velocity of Horizontal Tail Aerodynamic Center 



= Up 

+ 

ZrtQ 

^HT 

II 

+ 

XrtR “ Zjj,pP 

WjjT 

- Wp 

- 

XrtQ 


(C-5) 


Velocity of Vertical Tail Aerodynamic Center 


^VT ~ '^P 

Vyip = Up + XyipR - ZyrpP 
'"^VT ~ ^ 

Velocity of Left Wing Aerodynamic Center - Body Axes 

u£w * Up + Q (ZwAC ^ILwaC^ ^ ^WAC^ 


(C-6) 


7 ^ ar 

^LW 


%AC^ " ^ 1 l (C-7) 

^WAC 


w. 


I - 


LW 


"p - ^WAC'’ - W * ^ 


■’WAC 


4 


where ^Il^ac elastic deflection of the left wing aero- 

dynamic center. The equations for the right wing are obtained 
by substituting 

^Rwac “ “^bWAC 

^1 Rwac ^ AC 


C-3 
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Velocity of Left Wing Aerodynamic Center-Chord Axes 

In order to compute wing angle-of- attack the velocity components 
are required relative to the wing chord line. If the wing 
chord makes an angle i^, with the body centerline then 

ULW = uiw cos iv, - wf^W sin iv, 

Vlw * ’^fiW (C-8) 

''lW = ''LW cos iw + w£w sin i„ 

The equations for the right wing are obtained by changing the 
subscript. 

Velocity of Left Rotor Hub - Body Axes 


u' 

RL 


+ RY„ - L 
N s 

^k 

H 

+ Lg (R cos 

«k 

= Wp 

- PYn - Lg 


<iNL ^NL 

ifjL + P sin %l) - Phi^ 
(iNL + Q) cos i^L + 


(C-9) 


where Lg is the distance from the rotor pivot point to the 
rotor hub and hij^ is the deflection of the wing tip. The equa- 
tions for the right hub are obtained by changing subscripts and 
substituting = -Yn. 


Velocity of Left Rotor Hub - Shaft Axes 

Since the rotor aerodynamic forces and moments are functions of 
the shaft angle of attack and sideslip, the velocity components 
are required relative to shaft axes. 


'IRL = 'iRL cos ifiL - ^NL 

^RL “ '^RL (C-10) 

'^RL * ''k ^NL 

The corresponding equations for the right hub are obtained by 
changing the subscript. 

C .2 Center of Gravity and Inertia Equations 

Equations are required that express the overall aircraft center 
of gravity position and inertias in terms of the centers of 
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gravity and inertias of the individual mass components. In 
order to do this a fixed reference point is chosen in the air- 
craft defined by the intersection of the line joining the 
nacelle pivots and the vertical plane of symmetry of the air- 
craft, see Figure C.l. A set of axes is taken at this 

pivot reference point, parallel to the axes OXYZ at the air- 
craft center of gravity. If the location of the aircraft 
center of gravity with respect to the pivot reference axes is 
(Xqg» ^cg» 2'cg^ hy,) are the x and z 

coordinates or the fuselage and wing masses measured from the 
pivot, then the following relationships are obtained between 
the centers of mass of the components and the aircraft center 
of gravity. 


Fuselage CG Relative to Ai r craft CG 

Xf « if - X^2Q 


(C-11) 


Xf = hf - Z^G 


Wing CG I^elative to Aircraft CG 


^ ” ®'w ” ^CG 
^w “ N# “ ^CG 


(C-12) 


Nacelle CG Relative to Aircraft CG 



i 

COS 

Unr - 

- ^CG 

^NL “ 

1 

cos 


- 4g 

^NR “ 

1 

sin 


" ^CG 

*NL * 

1 

sin 

(i^L ” 

“ ^CG 


where i is the distance from the nacelle pivot point to the 
nacelle c.g., and > is the angular depression of the nacelle 
center of mass below the nacelle pivot, when the nacelle is in 
the down position, see Figure C.l. 
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Aircraft Center of Gravity Position 

By taking moments about the pivot, the aircraft center of 
gravity is given by 


» _ ”*f ^f r , 7 

m \-nrj cos(ijjjj-X) + cos (ij^j^-X)J 


ZcG — 


mf hf + m^ 


m 


K) 


(C-14) 


sin(ijjj^— X) + sin(ijjj^— X) 


The equations of motion (Section 3 ) require the first and 
second time derivatives of the center of gravity position. 
They are as follows: 

Center of Gravity V elocity Relative to Pivot Point 


• I 


^CG ^ 


•I 


m 


iNRSin(ijjj^-x) + iuLsinCi^L-X) 




(C-15) 


|ijjRCOs (ijjj^-X) + 

Center of Gravity Acceleration Relative to Pivot Point 


^CG 


= -z 


ir^ 


m 


,2 


iNRSin(ijjj^-X) + iNLsin(iuL-X) + Lj- cos 


(iNL~^) + ifiR cos 




(C-IC) 


■! ■' 
i : 

i ' • 


'CG “ 


•2 / — 
m 


iNRCOs(ijjR-x) + iNLCOs(ifjL-X) - 


ijjLSin 




NR- 


‘•NR"^^ j 


Pilot Station Velocities - Body Axes 

The velocities at the pilot's station are required in order 
to drive the visual display. From Equatiwns (C-3) and (C-4) 
the components of velocity of the pilot's station in body 
axes are : 


C-6 


) 
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u 


PA 




- ^^PA 


^PA “ ^P ^*'PA " ^^PA 

WpA » Wp + PYp^ - Q4p^ 


C-3 Pilot Station Acceleration - Body Axes 

The pilot station acceleration is also required to drive the 
visual display. These accelerations are derived here. 

The velocity at the pilot's station is 


VpA = VcG + a X rpA + 

where rp^ is the vector from the aircraft CG to the pilot's 
station and *^£PA is the rate of change of the pilot's station 

with respect to the aircraft CG. 

The pilot's station acceleration is 


^A 


^G 


= f^A = ^ZCG + 1_ (Q X £p^) + 1. 

dt dt dt dt \ fit / 

+ (£ X £p^) + n X (n X £p^) + + £ 


fit 

fin 


fit 


^ £pa 
fit 


a + !M X tj,, + 2i) X + a (r^ .£) - il'r + 1 !|pa 


with £ = Pi_ + q5 + 

II . Pi + Q3 + 4 


^A 


SrpA 

fit 


Sa - ='CG> i * '"'PA - ’'cg' 2 * 'Sa - ^CG> i- 
(Xp, - X-.) 1 (Yf , - Y„) X ■ ^CO* - 


C-7 




D238-10002-1 

and noting that Yqq and the time derivatives of Xp^, Yp^, 

Zp^ are always zero^ the cdjove equation yields the pilot's 
station accelerations as: 


~ — (Q “^CG^ + (Q + R ) (Xqq - ^PA^ 

+ Yp^ (PQ - R) - 2 QZcg - X cg 

ayp^ = ^AERO- f (P - QR) (Zcq “ ZpA> + (R + PQ) UpA " *cg) 
m 


- Ypj^ (R + P ) + 2 (PZgg " PXcg) 


^AERO 


‘ZPA 


t (6 - PR) (XcG “ J^Pa) + + Q^) ‘ 2pa) 


m 


+ Yp^ (P + QR) + ZQX^g ” Z^G 


where a 


^CG 


=‘AERQ 

m 


etc. 


t 


and Xp, = UpA, the distance from the pivot tc the pilot's 

station 

C . 4 Aircraft Inertias 

The aircraft roll inertia about the aircraft center of gravity i 

is, from the parallel axis theorem, 


I =1^ 

XX XX 


+ l''' + 

XX 


.NL 

^XX 


+ I^R .K 

XX 


ra^Zj + 


™w^w 


2mj^Yjj 


^2 

"'n^nl 




(C-3 


where 1^^^, etc., are the inertias of the various components 
about their individual centers of gravity. 
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In the case of the nacelles the inertias are dependent 

on the nacelle tilt angle, ifj* These inertias are related to 
the inertics of the nacelle with respect to a set of nacelle- 
fixed axes 0"xyz placed as shown in Figure 3.1. The re.'', ution- 
ships are 


t-N _ _N . ,xN _N , . 2. _ 

ijjx ~ ^xxr J Sin i^j ~ I.»„ sin 2i. 


ZZq -XXq 


"N ^xzo 


N 


[N = iN 

yy yy© 


'■22 


-I** + (I** 

^ZZq + UxXq 


. 2 

^zz_' ^xz 

o o 


(C-18) 


^xz “ ^xZq 2ijj + i (^xXq “ ^zz©^ 

Using equations (C-18) together with (C-13) , (C-11) , and (C-12) , 
in equation (C-17) , the roll inertia becomes 


xx 


^xx + ^xx + + (iN, - iN^ ) (sin^W + sin^i^To) 


‘ XX , 


■ZZq XXq 


NL 


‘NR' 


- I«o + 2 . mjhjZj 

+ - Vw4g 


’ “n^NL^CG ” "^^NR^CG 


in(ijjL - ^)j 


« + l'' + 2l” + ) (sin^i + sir^i ) 

XX XX XX - 22- XX -' NL nr' 


N 2 

- (sin 2i„T.-^An2i„B) + 2 m„Y„ + ro^zh^rZ, 


■•NL 


‘NR' 


“N"N 


‘f“f^f 


"Stf^w^w ” ^"'n I^^NR ^^NR ~ sin (ifj^ " 

since the terms containing sun to zero. 
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Similarly 


I ss (cos 2i-_ + cos 2i™) 

xz xz xz xz NL NR 


+ (I^ - ) (sin 2i + sin 2i ) + 


2 XXn zz 


NL 


NR f f f 


+ m^Z^£^, + imjj j^ZjjR cos (i^R - 
t^zz” lyy) = ^zz “ ^yy ■*■ ^zz ~ ^yy ^ (3 ^zZq “ ^yyo^ 




+ - I^Q^^sin %L + sin iRR)+ Ixz^^^^in 2i^L 

+ sin “(m^h^Z^ + m^h^Z^) + nijjJl sin (ijjj^ - X) 

1 ^ 

^NR sin (ijjR - X)J + 2mj,Yj, 


Similar expressions are obtained for lyy and l^z and these are 
presented in Appendix E . 
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APPENDIX D - CALCULATION OF SLIPSTREAM-IMMERSED WING AREAS 

The wing areas washed by the rotor slipstreams are required in 
the calculation of wing lift and drag. These immersed areas 
depend on rotor shaft inclination, wing angle of attack and side- 
slip, and rotor thrust. The equations presented in Appendix E 
for the immersed areas Sij^ and Sij^ were obtained as follows. 



The above sketch shows a rotor under conditions of combined 
angle of attack (a<r.L.^ sideslip (S) . The resultant angle 

of attack of the shaft is given by 


aR = cos'Mcos otrj, L. cos e) (D-1) 

If the rotor shaft is inclined to the fuselage centerline at 
angle ijj and the fuselage is at angle of attack aj.- then 


“T.L. ^f ^N 


(D-2) 


The rotor "sideslip" angle, is defined by 




Tan"^ 


Tan S 
Sin 


(D-3) 


and is the angle shown in the sketch. 

Figure D.l presents four views of the geometry of rotor slip- 
stream/wing planform interaction. 

Figure D.l [a] is a view of the plane taken through the rotor 
shaft parallel to the aircraft vertical plane of symmetry. The 
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line PT is the wing chords the distances PC and hp are the hori- 
zontal and vertical coordinates of the pivot measured from the 
wing leading edge# and i is the spinner- to-pivot shaft length. 

Figure D.l[b] is a view teOcen normal to the rotor disc plane. 

In this view, the traces of the slipstream on planes taken 
through the wing leading and trailing edges parall .il to the 
disc plane appear as circles. This assumes that the slipstream 
is a sheared circular cylinder. 

Figure D.l[c] is a section taken in the plane containing the 
rotor shaft and the freestream velocity vector V^. The angle 
e is the deflection of the slipstream relative to the freestream 
direction. Planes are taken through the wing leading and trail- 
ing edges parallel to the rotor disc. These intersect the rotor 
shaftline at the points 0 and T, and intersect the slipstream 
centerline at the points 0* and 0". These points enable the 
slipstream traces shown in (b) to be constructed. 

Figure D.l[d] is a view taken perpendicular to the wing surface 
showing the areas washed by tne slipstream. For convenience, 
this view combines the immersed areas of both left and right 
wings. In general, the imprint of the slipstream on the wing 
will be bounded in the chordwise direction by curves lines; 
however, the approximation is made that these lines are straight. 

The immersed area of the right wing panel is (assuming that the 
tip is immersed) , 



"iR- 

1(PM + 
2 

TN)c 





s 

i(PR + 
2 

RM + TS 

+ SN)c 


(D-4) 

From 

Figure 

D.l[b] 

PR » 

00' sin 

C 

(D-5) 

From 

Figure 

D.llc) 

00' * 

* (i-OD) 

Tan (oj^-e) 

(D-6) 

Prom 

Figure 

D.l (a] 

OD »» 

PC COS 


(D-7) 

Prom 

Figure 

D.l(b] 

RM = 

R'M' « J 

f°s - O'R'2 

4 

(D-8) 

From 

Figure 

D,l(b] 

O'R' 

“ 00' COS ; + OP 

(D-9) 

Prom 

Figure 

D.lta] 

OP = 

PC sin 

(ijj- iw) + hpcos (ij,- 

iw) (0 
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These equations define the leading edge intersection PM. If 
RM is zero or negative, the slipstream does not intersect the 
leading edge and the wing is considered to be miaffected by the 
slipstream. 


For the trailing edge intersection, TN; 
TS = 00" sin q 

00"= (A. + c cos (ijj-i^^) -OD) Tan 

SN = S'N' = ^ “ 0"S'^ 

4 


(D-11) 

(D-12) 

(D-13) 


(D-14) 

(D-15) 


SN 


0"s' = 00" cos ? + TT' 

TT' *» OP - c sin (ijj-ip^) 

If we write 

= PR, = RMf ?3 * TS, and 54 

then, using the above equations, 

[i-PC cos (ijj-ij^)+ hp sin(i^-i^^)] Tan (aj^-e) sin 5 (D-16) 

and 

« 2 - 


V - 


{[i-PC cos(iy-i^^)+hpSin(ijj-i^^-) ]Tan(a^-e)cos ; 


+ PC sin(ijj-iw) + hpcos (ij^-i^) } 


(D-17) 


The corresponding equations for 53 and 54 are obtained by re- 
placing PC in (D-16) and (D-17) and (PC-c) 

Thus the immersed area of the right wing panel is given by 

SIr = i c (Cl + C 2 + 53 + 54) (D- 18 ) 

From the symmetry of Figure D.l[d], SN=BS and RM*AR. The total 
immersed area of both wing panels is 

Si^ = i c (AM + BN) = j c (2C2+2?4)=C (C 2 +^) (D-19) 

and therefore the immersed area of the left wing is obtained 
from 


Si^ ” 


(D-20) 


The cdDOve equations correspond to those presented in Appendix 
E for calculating immersed wing area. 
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APPENDIX E 

The equations, data, and control diagrams required for FSAA 
simulation of the Boeing Vertol 1985 Tilt Rotor Transport are 
presented in the following pages. The simulation block diagram 
is shown on page E-6. Each element of this diagram is numbered. 
The reference table on page E-2 lists the block diagram element 
number, the function of the element, and the starting number 
of the pages containing the equations for the element. 

Data for the simulation is provided in Appendix F. 
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DENSITY C3\IiCl)LATI0NS 



INPUT: h 

ATMIND 0 STD ATMOS 
1 HCT ATMOS 

} 2 TROPICAL ATMOS 

r 

j OUTPUT: 6, 0, a, M, p 

/ 





































THRUST MANAGEMENT 
SYSTEM 


1 , 

I 

» 


1 
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ROTOR CONTROL COORDINATE AXIS TRANSFORM 


""ICL = ""ICL S ^ ®1CL ^ 


Bt =-Ai^ sin (J>„ + Bi cos 

-L PT J-fT ^ p -L n 


^ICL = ^HL - ®1CL ^1 


Bi = Ai <sin 5 + Bi cos 5 

CL CL HL CL HL 


())p is the control phase angle. 4)p is positive 
for the control axis moved opposite to rotor 


rotation. 


Al^„ = Ai „ cos <|i„ + sin <t> 

CR ^CR P CR P 


BIcr ^IcR "^P ®lcR '^1 


'"IcR * '"J-CR 'hr* Sr 


®1CR =-''1cr““’' Sr* ®1cr 'i 
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CENTER OF GRAVITY CALCULATION 


C.G. LOCATION RELATIVE TO PIVOT 


V _ mfif + ^ ^ 

CG I ^ 

m 




cos + cos 


Z = 

CG 

m 




in + sin (ijjR-^) 


C.G. VELOCITY RELATIVE TO PIVOT 


* /^N 

~ 

CG im 


H.L * "nr *"nR-"’ 


Kr’- <■ (-1 

CG Im i 


‘■n '"ni.-"> * "nr '"nr-"> 


C.G. ACCELERATION RELATIVE TO PIVOT 


CG 


^ (!!n 

'm 


^CG " " ^ 


m. 


N 


m 


NL 


■NR 


“NL 


NR 


sin 

1 

:z 

•H 



cos 

(i ~X ) 
NL 

sin 

<^NR-"> 

+ 

9 

i2 

■^NR 

cos 





. 2 



cos 

<"nl •'> 


^NL 

sin 





2 



cos 

(i -X) 
NR 

— 

^NR 

sin 

(i „-x) 

NR 
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FUSELAGE PIVOT VELOCITY 


U 


U 


Z q 
CG^ 



D238-10002-1 


V = V + 
P 


Z p 
CG^ 


- 

CG 


Wp = W + 




- z 


CG 
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y' 
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VELOCITIES OF AIRCRAFT COMPONENTS 
LEFT WING A.C. VELOCITY - BODY AXES 

^LW “ ^WAC^ ^ ^^LWAC 

^LW “ ^WAC^ “ ^WACP “ P ^^LWAC 

^LW “ " "^WACP ^WAC^ ^IlWAC 

ROTOR WING A.C, VELOCITY - BODY AXES 


%W " ^WAC^ * ^WAC^ ^ 


RWAC 


^RW ^ ^iJAr^ ^ ^tJArP P 


"WAC 


WAC^ 


i-RWAC 


^RW ^ '^WAfP ““ ^WAC^ ^ 


'WAC 


WAC 


LRWAC 


LEFT ROTOR HUB VELOCITY ■ BODY AXES 

^RL ^ ^ ^NL ^ ^1 l 

^RL = ^P ^NL P P 

. * 

”rL = ”p - P ’'n - ^<i»L * ■J> =°= ^81. * >'1 l 

RIGHT ROTOR HUB VELOCITY - BODY AXES 


U 


RR 


'RR 


W. 


RR 


^P “ ^NR ^’nR '*’ ^^R 


Vp + Lg (r COS ijjp + p sin i^jp) - p hip 

^ • 

Wp + p Yj, - Lg (ijjp + q) cos ijjp + hlj^ 
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E-20 







I 


r 

^ • 

LEFT ROTOR HUB VELOCITY - SHAFT AXES 



Url 

= Kl 

COS 

i»L - "rL 

^NL 


Vrl 

= <L 





Wrl 

= Url 

sin 

isL ■" "k 

^NL 

RIGHT 

ROTOR HUB VELOCITY 

- SHAFT AXES 



^RR 

= U^R 

cos 

^NR " ^RR 

iNR 


Vrr 

“ ^RR 





Wrr 

“ ^RR 

sin 

■^NR %R 

^NR 

LEFT WING A.C. 

VELOCITY - CHORD AXES 




= \w 

cos 

- "lw 

^W 


Vlw 

= v{w 





Wlw 

= Ulw 

sin 

iw + ^LW cos 

iw 

RIGHT 

WING A.C. 

. VELOCITY 

- CHORD AXES 



Urw 

= Urw 

cos 

iw - Wrw sin 

iw 


Vrw 

= 





Wrw 

= '^^RW 

sin 

iw Wrw cos 

iw 
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D238-10002-1 ? 

i-^ 

HORIZONTAL STABILIZER A.:. VELOCITY ~ BODY AX ES " | 

-• I ; 

“ i ^ 

^HT " I 

i ' 

^HT " ” ^HtP I 

J 

I 

\ \ 

'^HT ” i ;■ 

-i 

VERTICAL FIN A.C. VELOCITY: RIGHT FIN (BODY AXES) i 

UvTr* ^P " ^VT^ 

VvT^= Vp + 

W^T = Wp - X^,^q + 

R 1 

LEFT FIN 

^VTl “ Zyipq + 

VvT^ = Vp + * ZyjpP 

^VTl ~ ^ ” ^vtP 


1 ’ 

4. : 

i : 

i 


f 



] 


I 


WING AERODYNAMICS 


CALCULAIE ROTOR INTERFERENCE TERMS: 


D238-10002-1 


+ Tan 

RR RR 


-1 


NF, 


^R I 


RR 




T2 + NF2 + SF^ 
R R R 


V 


V* = 
R 


RR 




Rj^r +10 


2pA 


4 2 V v3 cos T + vj '--I = 1 (Solve for v ) 

R R R R R R 


'^p = Tan"^ 

^RR 


^RR 


RR 


'TS 


RR 


V* + V. cos T 
L R R 

<^RR - °8R> 

<"rR - “RR> * ’^R 


LR 

\r 

V* 


= ct. , + Tan 
lR 


-1 


[=ft] 


= "\/t" + NF2 + SF2 

V L L L 


'LR 




Rlr1+io 
2pA 


4 2 v,_ v3_ cos 4 V3^ V3, . 


■L "L 


L L 


1 (Solve for ) 
L 


^Plr • 

[- 



n - COS 

'^’■Sl.R ■ 

^ 'LR 


<’lr - 'lr> * ''*1. 
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£ = Tan”"^ V /(W„„+ 


U ) 


RR RR WRR RR 


^HL ■ + 


/(W + e U ) 

RL'^ RL ^WRL RL 


USED IN ROTOR 
CONTROL TRANS- 
FORMATIONS 


5 


w ■= 




i -■ 

^R “ 

'“RR + 

“lr^ 




i 


^^prr'*' 

^PLR^ 

1 

USED IN TAIL AERO- 



^TS 

^^TSRR 

■*■ ^TSLR^ 

1 

DYNAMICS . 



II 

1 

'Sl ^ 

"nr^ 

1 

USED IN WING/ROTOR 



^LW 

(O-LSRW 

+ ‘^LSLW) 


INTERFERENCE. 



(1-Cts) 

J 

1 



PC cos 

- 

i ) + h 
w p 

sin (i^, 
N 

- i )] tan (a_ - e_) sin C 

V,' K r 




[Lg - PC cos - i^) + hp sin (i^ - i^)] tan 

- - 7p) cos ? + PC sin + hp cos 


(i-M - i ) 
N w 


ir: 0 or Imaginary; - 0 and - 0, 

also - 0 and ■ 0.0 

Form by replacing PC in equation with (PC - c^) 

Form by replacing PC in equation with (PC - ) 

IF:Ej ^4 = 0 or Imaginary; = 0 and * 0, 

also - 0 and ■ 0,0 

IF: UNBRELLAS OPEN; SET C = 0.0 

LU 

UMBRELLA LOGIC: 

IF ij^j^gp -■ F.j^ or qp > 8.479 Ibs/ft^ sec umbrellas 
closed (hysteresis + 1®; qp + *1 Ib/ft'). 
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- Cw/2 (Cri + Cr2 ^3 

(Si/S)R„ = 2 (Si^/S„) 


% = =w Ur2 + ^R4^ 


^LW 


(Sj^/S)Ly — 2 (Sj^/S^,) 


(ARi)LW = (Si /c2 ) 


I“i'rw = '=1.^=3 ) 


AR^ = S„/C^ (FROM PREPROCESSOR) 


(^Lcsi/'^La^LW ~ 


" "Law fl/(AR.)L,, - 1/AR„] 


(*^Lai/^La^RW 


^ (1/(AR.) - 1/AR ] 

L^w 1 Rw w 


Hl 


Kar 


V* + (C /C ) V* 

L Lai Let LW L 

V£ + v£ 

^R ^ (^Lai/^Lg)RW 

V* + V* 

R R 
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✓ 



n 


[1/2 p (u2 + v2 + w2) + (T^ + Tj^)/2A] 


[1/2 p + v|y + w|^^) + Tj^/A] 


[1/2 p (u2 + V- + w2 ) + T /A] 

LW LW LW L 
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WING ANGLE OF ATTACK AND SIDESLIP 
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-1 


'LWO 


'RWO 


= Sin 


. -1 
= sin 


W 


LW 


. ^ LW LW 
‘ ^RW 

^RW 


+ 0 


+ 9 


tLWAC 


tRWAC 


^LWO 


’RWO 


, -1 
= Sin 


^LW 


= sin 


-1 r 


/u£w 

+ + 


i 


'^RW 

- 

i 

'u^ 

+ + 

W2 


RW 

RW 

RW - 



^LWSSO ^LWO ^PLR 


°^RWSSO *” ^RWO ^PRR 


^^LWO ^RWO^'^^ 


a 


LW RIGID 


'rw rigid 


= Sin 


= sin 


-II 


W 


LW 


Vu£w + w£w 


w 


RW 


OJKv * "rw j 


'PLR 


- e 


PRR 


^LWO 


"^LWO ^W " ^tLWAC 


RWO 


a - i - 6 
RWO W tRWAC 
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CALCULATION OF INCREMENTAL 


' I j. 
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DRAG AND MOMENT COEFFICIENTS 

= '^aLw ^ ^f' '^SP = 

= + «f' «SP ' 

= 6 + 6/6 =6 

aRw f SP SPr 

= 15- ■'■<5-, 6„p = 6-p 

apy f SP SPp 


^LW f 

SP 

S^L 

+ < 5 . 

, 6 

= 6 

^LW ^ 

SP 

SPl 

+ 5 

/ 6 

= 6 

^RW ^ 

SP 

SPr 

+ 

» 6 

= 6 

^RW f 

SP 

SPr 

f 

' Sp 

= 0 

(0°l 5 £ 

62) 


(o 2 < 5 <_ 

63) 



(5 > <SJ 

j 

(0^5^ 55 ) 
(5 > 65 ) 


28 
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" ^14^SP 


“ ^15 ^16Sp ^17^|p 


(0 I'SP I'SP ) 


(6 >6 ) 
SP SPj^ 


AC^ = + b,6 


SP 


O^'SP "1"SP 


Fj^ + F2A + F362 


= F + F,6 + F 6' 


(0°< 6 < 6.) 


(6 > 6^) 


“nl 


ao + a^A 


(0°^ A ^ Aj^) 


NL 


— 5 L #■ 


= S + »4« 


= a. 


(A > Aj^) 
(0°1 6 1 


(A > Aj^) 


If a“j^ 1 “ 1 “nl calculate: 

C. = a^ + C' a + AC ■ + F AC, 


'La 


W 


SP 


C,„ = C- - F AC 

LWi L Lgp 


*^LW2 ■ ^^L^ 


(a ^ a^) 

= a, + ci^ a, + a,, + a,^ a + a^a^ < a 3 ) 


^LW, = ^6 ■" '-La ‘‘S " ^23 - “24 " ' “25 
2 w 


= + a^£ C^„ + a,7 C^ ,j + a,j, ^^t.u " ^lw ^ 


"'"'2 


"27 ''LW^ ' 28 '“LW^ 


'SP 
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• • 


If < a < a"*”, + a,- calculate: 
NL ““ NL lo 


C = a + C' a"*" + AC + F AC 


Ljjl 6 NL 


L6 


'SP 


a = a - a + a 
DUM NL 0 


AC = a„^ + a„ a + a„ 

L„, 20 21 DUM 22 DUM 

NL 


C, = CJ + AC. 

^ ^NL ^NL 


C, „ = C, - F AC, 

LWi L Lgp 


*^LW2 ^6 “ 


(a 1 a^) 


a^ t C. .. a« + t a/,, ot t a^,, ct 


'La “0 “20 “21 

W 


22 


(a > aQ) 


^DO 1.7 ^0 7 Ct t.T t 3.»5(j (C, „ “ C, ,, ) 


w 


‘26 ''LW 


‘27 ^LW]_ “28 ""LWj^ 


LW. 


5 SP 


and print stall warning. 


If + a^^g < a <_ 90® calculate: 


C, = (a, + C/ + AC, + F AC, ) (90° - a)/(90° - - 


La,, NL 

w 


La 


^SP 


NL 


^8^ 


+ a „ 

2 0 18 


J 


C,„ = C, - F AC, (90° - a)/(90° - a+ - 


'LWj^ L 


'SP 


NL 18' 


(a ^ a 2 ) 


= C, - F AC, (90° - a.,)/(90° - - a^^g) (a > a^) 


SP 


*^LW2 “ ^6 *^Lay ^0 + ^20 + ^21 ^22 


(a < a, ) 


“i * ^Lc,, ^0 

w 


(a > a 2 ) 
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'• ^DOy ^26 ^ 1^2 ^27 ^LWj^ ■•■ ®28 


+ AC ^ + AC^ 

DO^ Dgp 


C = C 
D Di 


(a < a^) 


C =C +(l-C)(a-c( )/(90° - a ) (a > a ) 
D 2 2 2 


and print stall warning 
If 33 - a^9 £ a < a^L calculate: 


CLk,t = ^6 


+ Ct 


'^NL '' 

<^DUM = a “ «NL ■'■ ^3 


“NL ■’■ ^^15 ^ ^^LsP 


®23 '*’ ^24 “duM ^25 "^DUM 
NL 

Ct = Cf + ACt 
^ ^NL ^NL 

Ct = Ct - F AC, 

^Wi ^ ^SP 


W 2 


“NL ■*■ ®23 ■*■ ^24 ■*■ ^25 

“w 


Cd - + 326 + 327 + 328 


- Ct ) 
^Wo 


+ ACj^ + AC_ 
*^05 °SP 


and print stall warning. 
If - 90 ® ^ ct £ 33 - aj^g calculate: 


“1 “ *^NL ■ ^19 


“nl ■*■ ^ 


^NL ^ ^6 "SP 
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“dUM “ “ “nl ^3 


AC, 


^NL 


” ^23 ^24 “dUM ®25 °‘dUM 


= C; + AC, 
^NL ^ 


NL 




“DUM 

= 2^2 “ ^19 

AC, 

‘'NL 

= 323 ■*■ ^24 


'Wi 


-W 2 


= Ct' + ACt - F ACt 
‘'NL ‘'NL ‘'SP 


“ ^6 ‘^Lq, ®3 


(a^a^) 

(a<a^) 




+ ACn + ACn 

*^06 SP 


— ~ ~ ®3 ^19^^ (90 + ^2 “ ^19^ 

and print stall warning. 
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CALCULATE ; 
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i 




® - “lW / 6 - + 6 ^ 

•> LW 




“ ^Vf ® “ “ “orr / 6 = 6, + 6 

“ '"sSj, f »RW 


"RW 


* I 


= Cm @ a = , 6 = 6 f + 6 . 


LW 


* » 


^MrWq = Cm @ a = aRWo , 6 = 6f + 6 ^^^ 


AS FOLLOWS: 


If a < a < a 
1 ~ ~ 2 

Calculate C = b + b o 
M 2 3 


If a > a, 


ACj^ — b^ + b 5 6 + bg 6 ^ 

6 




C^ - b2 + b3 ct2+ AC^j 




If 


ot < ot . 


^ ^ ^ ^ Si, 


Cm = (90 + a)/(90 + a^) 
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CALCULATE : 


t f f 

^L 

= C ; 

*'LW0 

1 t 1 

^LW 



t f 1 

Ct 

^RW 

~ Ci ; 

*"RW0 

f t 1 

^Drw 

^^RW' 

- 

‘^»RW 

* 

^^LWMAX 

^^MAX 

+ ACl^ 




^^RWMAX 

^^MAX 

+ ACl^ 

+ ACl^„ 
SP 



p* 1 

^^RW 

= c? 
lrwo 

ct' 

°RW 

= Cq' ; 

°RWO 

C* — 

^RW 

“^^RWO 

Cl' 

^LW 

= Cl 
‘'LWO 

Cd' 

*^LW 

n* * 

““ C n f 

‘^LWO 

Cm = 

LW 

Cm' 

LWO 

c* 

^LWMAX 

= c£ 

^LWMAX 




c* 

^RWMAX 

= Cl 

^RWMAX 




“ '=1-0 

*"LW 

*IGE 

(ag/^)v' 

f 1 t 

/ VT 

- m2 



= CLj^y (ag/i)w/ Vl-M- ; 

1 

^"IGE_ 

'^RW 

f f 1 

^’"RW 

= Clj^j. (ag/a)w'' 

yi-M^ ; 

P*IGE 

^Lrw = 

*' 

*“^RW 


(a.g /a ) w/ 

(ag/a)w/'\/l-M2 
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AC’ 


;ge 


IGE 



IGE 

(C* 

LLW 

- c’" )2/ttAR ; 

LLW w 

IGE 

AC” 

DRW 

= ^9 

IGE 

(C" 

LRW 

^99 

IGE 

(C* 

LLW 

- C*' )2/7tAR ; 

LLW w 

IGE 

AC* 

DRW 

= ^9 

IGE 

(C* 

LRW 


XGE IGE 

1 C£lw max ' SET AC{5lW = 0.0 & C£l„ = C£l„ maX 


IGE 

TP. C” > C 

• ^LRW - LRW MAX 


IGE IGE 

; SET AC£j^y = 0.0 & C£j^^ = 


IGE lOE 

IF: C* > C* ; SET AC* „ 

LLW ~ LEW MAX DLW 


IGE 

0.0 & C* 

LLW 


C* 

LLW MAX 


IGE IGE 

IF: C* > C* ; SET AC* 

LRW ~ LRW MAX DRW 


IGE 

0.0 & C* 

LRW 


C* 

LRW MAX 


IF: (ag/a) > 1.0 

; SET 

K„ = -1.0 

(ag/a) _< ]. .0 

; SET 

K„ = +1.0 

CALCULATE 

IGE 


c” 

C” 


LLW 

LLW 



1 1 1 

IGE 

C” = 

C + AC" 

DLW 

DLW 

DLW 


IGE 


C” 

C" 


LRW 

LRW 

t % 1 

IGE 

C" = 

C + 

AC” 

DRW 

DRW 

DRW 


IGE 


C* 

C* 


LLW 

LLW 

IGE 

C* = 

C* ' + 


DLW 

DLW 

DLW 


IGE 


C* 

C-‘ 


LRW 

LRW 

IGE 

C* = 

C* ' + 

AC* 

DRW 

DRW 

DRW 
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'LSLW 


^At 1 (s 1 Ly ^^LLW ^PLR " ^DLW 


^LLW ■ ^TSLR^ 


PLR ^ DLW 
S 


-P) 

'S ‘ LW. 


^LSRW ^ 


^ jif") 


RW 


• .'■ COS e - C" sin e ) 
LRW PRR DRW PRR' 


t- C* (1 - C ) 
LRW TSRR 




RW, 




DSLW 


Al 1 's 


K’ < 1“; — I (C" sin e + C" cos e ) 

* ' ' LLW PLR DLW PLR 

LW 


+ C* (1 - C ) 
DLW TSLR^ 


1 - -I 

's 


LW J 


"DSKW 




RW 


^^LRW ^PRR ^DRW ^PRR^ 


+ C* (1 - C ) 
DRW TSRR 


RW 


^MSLW 

“'ij 

1^1 
's ' 

LW 

<=5l«i 

■'■ ^MLW “ “"TSLR^ 

■-(M 1 

LWj 

C 

MSRW 


l^il 
's ' 

RW 

<=HRW> 

^MRW ” *^TSRR^ 

ri-(!ii 1 
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i r ^ C* 1 i 1 -1/2 (---) ,,i 

"‘'Xs fOUER “ 1 V 

-IC,3R„ -U - ?TS> C£rw' !l-V2(|i)Rwl| 

AC = 1/4'Vr« 

NS POWER I^SRW 

-'=DSEW - Ss’ "Slw' 


c = (K + K c ) ( 1 - 

^;^SW 20 2l L TS - 


\ 2bw 1 


'NSW 


'^*L« ■ ‘^LRw’ 


= (K.. C,^ (1 - C^j) Bj + 


22 L 


1 1-Cts1 

(K„> [ 

\ 2b„ 1 

" 1 


DRW DLW 


■^^JrW ^““rWO ■ ^'LLW ( AC 


+ AC 


NS POWER 
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sp:jcial conditions (for umbrellas open) 

IF: UMBRELLAS CLOSED; GO THROUGH KING EQUATIONS 

IF: UMBRELLAS OPEN; CALCULATE THE WING FORCES AND MOMENTS AS 

— FOLLOWS : 


= fe„ q 
AERO ^ 


vx-v- / r ~^LW 
SLW ''tSLr' [] 


^AERO ■ ^SRW ^TSRR^ 


-U 


RW 




LW 

‘aero 

RW 

‘aero 

LW' 

^AERO 


0.0 

0.0 

Tl (D/TIl ( 


^.AERO '^R 


GO TO WING BENDING 


M' 


.LW 

AERO 


.LW 


M 


.RW 

AERO 


= 2 aero ^ ^L 

2 

=-^c ^aLo + ^r J 


LW . RW 
^AERO ^AERO 

WING BENDING 




.W 

'AERO 


Z 


RW 

AERO 


1-(S /S) 

i RW 


- I 


LW 

AERO 


l-(Si/S)Lw] 


W 

^'aero 


= 0.0 


IF: (h/D]j,ppl 1.3; (D/T)^ = [h/D]£pp + [h/D]ppp + K 

T" LR T LR 


'EFF- 

LR 


& (M/T)j^ = [h/D]ppp + Kj^ 2 t^/2i]£pp + ^ m3 

t“ LR ~ LR T 


D2' 

T 



! 

i 
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1.3; (D/T), = K_.; & (M/T) = K 


IF: lh/D]gpp> 1.^; "'^04' “ "M 

LR T T 


[h/DJpr.p+ l'n2 ^D3 


IF: [h/DlEFFl <D/T)^ = ih/uJeff" “‘/-'EFF _ 

FR <t. RR T T 


S- (M/T)j^ = [h/DlEFF + [h/Dlgpp + 

t” RR T RR T 


IF: [h/D]ppp > 1-3; (D/T)p - & (M/T)p - 

RR T T 


/ 

I 

; 
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FORCE AND MOMENT TPANSFORMATIONS 
FROM WING A. C. TO ELASTIC AXIS 

PITCHING MOMENT 


M^'^ = C q W c„ 

AERO MSRW ^SRW 2 ~ 


- X 


2RW 

’WAC AERO 


+ z 

WAC AERO 


wLW — r> ^ Sm 

\eF0 “ ^MSLW ^SLW 2 “ 


X 

WAC AERO 


+ Z 


XLW 

'WAC AERO 


VERTICAL FORCES 


,RW' ^ 

'aero 

jLW’ _ 

'aep.o 


-C 


LSRW 


C CL ’ 

DSRW RWO 


-C 


LSLW 


- C 


'DSLW ^LWO 


Sw 


^SRW 2 




^SLW T 


NOTE: Z^'^' & Z^'^' ARE USED IN VERPICAL BENDING EQUATIONS 

AERO AERO 
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WINr FORCE & MOMENT RESOLUTION - BODY AXES 0 C.G, 


X 


LW 

AERO 


AERO 


vRW 

■'^AERO 



r 

1 


: -C + C 

a * 

q 

i DSLW LSLW LWO 

^SLW 

L 

r 

J 

1 


"•^DSRW + ^LSRW “RWO 

J 

^SRW 

— 

1 

n 

i 



L DSLW LWO 

^SLW 

2 

."'"DSRW ^RWO j 


Sw 

2 


FROM VERTICAL BENDING 




+ X 


’W 


2 




AERO AERO CG (AERO AERO 
- Z 


CG 


LW RW 

AERO AERO 


N 


W 

'AERO 


C q S b 
NSW W W 


NOTE: OBSERVE WING SPECIAL CONDITIONS. 


f 
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HORIZONTAL 


TAIL AERODYNAMICS 


WING AND TAIL ALTITUDE - GROUND EFFECT 
h 


Wc/4 ~^D0WN 

^^WAC 

Tc/4 ~~^D0WN ■*■ 

«HT 


- Sin 0 + - Z„^^) cos e 


CG' 


HORIZONTAL TAIL ANGLE OF ATTACK 


^AC ^WAC “ ^HT 


CG “WAG' 
'CG “ 


(FROM PREPROCESSOR) 

GEF = [fa2 + 4(h , - h , )-]/[h- + 4 (h + h )2] 

W Tc/4 Wc/4 W Tc/4 Wc/4 


IF: 


> e + de/da (a - I W/U^) 
o W AC 


then ^ (1 - GEF)/"\/*l - ll- 


IF: £p < + de/da (a^ - I W/U^) 

then e = + de/da W/U^ ) (1 - GEF)/"\/l " M‘ 


WHERE dt/da = de/da0 + OfR)/2 




(TT<0) 

= Tan-1 fi'^nx/UHT^+ Ihi 

This form for aj,^ is to be used for resolution of forces only. 
IF: I j I > 180° then calculate a from 

c<HT "(sign aj^.p) 360°+ 

and use this value to obtain the forces and moments. 
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, TAIL 

LIFT AND DRAG 


^HT 

= Or^ + Tr^ 5g 



= (ot^j'P “ 2° ) 

"^STALL 

+ XR.J. 

^HT_ 

(cirt ” 2 / 

^^STALL 

+ Trt 

^La 

= Clort (ag/a)j^^/ 

VI^ 


Where (ag/a)^^^ = ^^^TC/4^ 


IF; 


“hT_- - ''•HT+ 


^Lht 


COrt - ‘-DOrt + CLr^/^ARrjEr^ 


IF; 


a < a < 90 ‘ 

HT+ 6rt - 


Cl 

HT 

” ^La 

«HT^ 

(90 

Clrt 

= 

Cl a 

“HT, 

STALL 



CDrt 

STALL 

CdOrt + 

^°HT 


STALL 



•HT 


/w - c„ ) 




'HT 


'HT 


) 


STALL 


90° - a 


HT+ 
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HORIZONTAL TAIL LIFT AND DRAG (CONTINUED) 


IF: 90“ < <(180“ - .5 otuT-) 

®HT “ 


C, = .5 C, (X 


^HT 


La HT_ 


(a^ - 90“)/(90“ - .5 ) 


= HT 


'HT 


= .5 Clqi a^T 


STALL 


03 


HT 


STALL 


'HT /’’ARh^Eht + Cq 

'^^STALL ^HT 


CDrt ^DfjT 


STALL 


+ ( + *5 qht_ ■ 180°|(1.1-Cd 


II 


^HT 


STALL 


(.5 anx - 90“) 


IF: (180“ - .5 (x„„ ) < a„ < 180‘ 

HT- - _ 


CLjj^ = Cl« - 180“) 


Cn = Cnn + Ct /mAR„„E„_ 


IF: 


"“90 ^ ^HT 

_ Hi. 


Ct _ = Clo, ^^hT- (-50“ - )/(-50° - Cht ) 


'HT 


■HT 


HT 


Cl 


HT 


= CLa «HT 


STALL 


'Dht “ ^DOhT ^LHT T 

STALL "stall 


O'D — O'D 

^HT HT 


STALL 


HT 


‘^HT_ 11^*^ " 


HT 


STALL) 


(-90“ - anx ) 
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HORIZONTAL TAIL LIFT AND DRAG (CONTINUED) 


IF: (-180° + <-90‘ 


Cl — .5 Clq <^hT ^*^6 90°)/(-90° + .5 cxu't' ) 

i^Ht ^ ®hT + 


^Lut “ *5 Clu aHT_L 

*^^STALL + 


= 


'^HTstALL ^ ''^HTstalL 


Cn.._ = C 


^HT ^HT STALL "" ^ 


- (a,,_ + 180" - .5 aHT.. ) (1.1 - C 


HT 


-D 


H T c T 


) 


STALL 


(.5 - 90°) 


IF: -180° <_ “ejj-p "^(-180° + .5 anT^) 


Cl = Cl« (cie + 180°) 

Lri Lot 6jj^ 


CDht " ^DOrt Cl^^/^ARhtErt 
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VERTICAL TAIL AERODYNAMICS 

VERTICAL TAIL ANGLE OF ATTACK AND SIDESLIP 


6 


VT 


Tan 



+ W 


2 

VT 


+ Sf (do/dB) {NOTE; 


THIS VALUE OF aVT IS USED 
IN RESOLUTION OF FORCES 
AND MOMENTS) 


IF: 


1 1 >180° -(sign a^^)(360°) 


a 

evT 

^'^vt'*’ ^VT 

RUD^ 


“VT+ 

^“'^ISTALL 

2 )‘^'^VT 

'^RUD 

P > 
< 

H 

1 

STALL 

2°)+Xy^ 

'^RUD 

^Ya 

’ Vl-M" 




(NOTE: THIS VALUE OF 
avT ONLY USED 
IN CALCULATION 
OF FORCE AND 
MOMENT COEFFI- 
CIENTS) 


TAIL DYNAMIC PRESSURE AND SIDEWASH 


q = p/2 (U“ + + w2) 

a =(do/dB) Bp 
VERTICAL TAIL LIFT AND DRAG 


IF: ^ 


Cy - Oyp, Cl- 
iVT *=VT 


Cn — Cnn 0^ /tt AR E 

^VT ^'-'VT ^VT VT VT 
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VERTICAL TAIL LIFT AND DRAG (CONTINUED) 


IF: ■' 

VT+ - 


'VT 


c = c 

^VT 

“VT^ 

(90 

Cy 

VT 

STALL 

“ ^Ya 

“VT. 

^DvT 

^^STALL 

” ^DOvT 


vt; 


Cd ~ Cp 
Uyi i^VT 




+ (a. ^ - CLyrj* )( 1.1 “ Cy 


STALL 


'VT 


'VT 


) 


STALL 


(90° - '^VTi.) 


IF: 


90 ° < a. < ( 180 ° 

- ,5 

VT 


^Yvt ~ ^Ya “VT_ 

(a 

evT 

S = -5 Sa 

1 ym 1 U 

STALL 

Ct,Tm 

VT_ 

Cy = Ct)Q 

'“^STALL 

9 

+ Cy 

Cd = Cn + 

^ ^STALL 

(otp 

^VT 


"VT 


‘VT- 


VT 


/ ir ARy'j-Ey'p 


STALL 


- Cd 


VT 


(,5 - 90°) 


IF: (180° - .5 avT^) 1 < 180° 

C\ 


= ^Ya ^ 

itTfTl ^TT«-P 


' VT 


'VT 


CDvt ^DOvt Cy^^^/ttARvt^VT 




VT' 


) 

STALL 
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VERTICAL TAIL LIFT AND DRAG (CONTINUED) 


IF: -90° ^ ^ 

- VT 


Cv = C,. )/(-90° - ) 


^Ya ^VT 


-VT 


Cy 


VT 


= Cya otVT 


STALL 


^VT 




STALL 


VT 


^VT 


STALL 


/ TT ARy y E y .j. 


Ct^ ““ L-rv 

^VT ^VT 


“ Q^VT-^ ^ 


STALL 


*VT 


>VT 


) 


STALL 


(“90 “ ^VT ^ 


IF: (-180® + .5 < -90‘ 

Vi_j, 6 vt 


^Yvt “ ^VT^. ^^evT 90^^)/ (-90® + .5 


Cy 


VT 


= ,5 Cya avT 


STALL 


'VT 


- C|)0-^ “*■ C'^ 


STALL 


VT 


■VT 


/ ARy yEy -p 


STALL 


= C. - (a^ + 180® - .5 a ) (1,1 

"^TstalL VT VT^ 


- c. 


'VT 


( • 5 ^ V T^ ^ 90®) 


IF: -180®< cte < (-180® + ,5 ayy ) 


- '^VT 


Cy — Cy^ (ap + 180® ) 
^VT ^VT 


CY^^/iTARvT^VT 


) 

STALL 
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TAIL EQUATIONS LOG I C 


HORIZONTAL TAIL 


1. If h„ > 100 feet; set GEF =0.0 
Wc/4 


2- If the umbrellas open; set e = 


Ep (1 


VT- 


3. If a > a„„ print stall warning 
®HT 


4. If a ® T print stall warning 

®HT 


VERTICAL TAIL 


1. If a„ > a,,_, print stall warning 
®VT 


2. If a- < oivT- print stall warning 
VT 


- GEF) 
m2 


E-4<? 
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TAIL FORCE AND MOMENT RESOLUTION TO C.G. 


HORIZONTAL TAIL - NOTE ; - IF UMBRELLAS OPEN AND U>0; SET n, 

= .5 ’ 


^AERO ^ ^nuT ~ COS (6„„ ~ o) + C_ sin 


DHT HT HT 

(a — i ) 1 Q S n 
HT HT^ ^ HT HT 


AERO 


sin (B - a) q S n 
DHT VT ^ HT ^HT 


'AERO 


Clht cos cos (6^^- a) sin 


HT 'HT 


AERO 


""^AERO 


AERO 


HT u 

^AERO^^CG “ ^HT^ ^AERO^^HT " 


(X —X ^ 
AERO ^AERO ^HT^ 


VERTICAL TAII 


AERO 


aero 


AERO 


jj"*"DVT ^“^VT ~ (^)cos(ay^ - ij^^) -CYy^sln ( - a 

^^HT " ^HT^ j ^ ^VT "^VT 

= [‘'vvT <^T - -SVT <^T - ^ ^'T \- 

jj*^DVT ^^VT " ^ ‘hT " ^HT^ " *“YVT 

sin (i - i ) q S n 

HT HT J ^ VT VT 


n 

T VT 
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VERTICAL TAIL (CONTINUED) 


(Z - Z ) 
'AERO AERO VT CG 


M 


VT 

AERO 


VT 

7 (y 

^AERO 


VT 

Xyi) + ^2ER0 




VT VT 

^AERO ^ “’^AERO 


TOTAL TAIL CONTRIBUTION 


yT ^ VT „HT 

^AERO '^AERO ‘'-.fiRO 


"AERO 


„VT HT 

"^AERO "AERO 


M 


AERO 


VT ,.HT 

-^ERO '"aero 


T ^ VT HT 

AERO AERO ^AERO 


AERO -^AERO '^‘^AERO 


,^T 


^AERO -^AERO AERO 


E“5 » 
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NACELLE AERODYNAMICS 


NACELLE ANGLE OF ATTACK AND SIDESLIP 


'^RN 

= Tan ^ ; 

; ^RR'^^RRj 

' ^RN 

= 1/2 p 

^RR 

’^LN 

= Tan"^l 

1 

[^^RL-^^RL j 


= 1/2 P 

Sr 

“RN 

Tan U|j^ 

'’^rr] 



^LN 

= Tan'H 

_^rl/ v'^^rl 

* «l J 




NACELLE WIND AXIS FORCE & MOMENT COEFFICIENTS 

+ K^o i 1 + ^31 1 I NOTE : CHECK RANGE OF 

otRN ctLN TO 
DETERMINE VALUES 
CONSTANTS. 


*^DRN Son 
^OLN Son 


SrN - ^32 

Sin 

'•‘RN 

cos 

'^LLN = ^32 

Sin 

-LN 

COS L \ 


S,RN 

Sion 


Sin 

COS 

■‘rn 

^ S 5 

(sin .= j., 

CO -5 

''rn^ 




^RN 

cos 

'^RN ^ 





S.LN 

II 

O 

o 

+ ^ 3 ^ 

5 ^ n ^ "L ^ 

cos 

‘ln 

+ K33 

(sin 

cos 

‘LN^ 




sin 

cos 

^LN 

! 





SPECIAL CONDITIONS 

1. IF: 0(FT/SIC)S- RIGHT NACELLE AERO i 0.0 & 

HOLD Vf,LUE OF 

2. IF: (FT/SEC)"; LEFT N.ACELLE AERO E 0.0 & 

HOLD VALUE OF 
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^•RN ^36 ^RN + K 37 (Sin 0 ^^^ Cos Cj^^MSin 3j^j^ Cos 

'YLN " ^36 ^.N ^ "ln ®LN ®LN • 

'NRN " SoRN K 3 gSir. Sj^^^Cos + K 35 (Sin Bj^j^Cos&j^j^J 1 SinSj^^j Cos 


^NLN SoLN ^ ^40^^^ ^LN^°® ^LN '*’ 3, „CosS, „) I SinS. „ Cos 


LN LN 


5^'RN ^LN ° ° 

NACELLE FORCES & MOtlENTS - HACELLE AXES 



^RN^W^ ~^DRN^°®^RN 

'*' ^LRN 

■ 1 

RN 

^RN^W^^YRN ^°®^RN 

^DRN' 

RN 

‘^RN^W^“^LRN^°®=‘r>^ 

^drm' 

RN 

= a S b f ^ W p 
^R.'rw^w K ^ ^MRN 

sinS-,., 

RX 


RN 


‘RN ^NkN 


"^'"rN = ‘5r."/w f^RN ^RN 

''^RN ^ °RN^w\ ^^XRN "‘rN ~ ~ ^MRN =R.N 

bw ‘ 


i-N 

‘^LN^W ^“SlX 

=°= ‘LS * =U.S 

‘LN ■ ^YLN 

^-LN 

COS 11/2 


‘^LM^W ^*"YLN 

"LN* ‘ '"I.I.^^ 

^Ln11/2 



"2ln 


‘ln~^yln^^*^“ 

LN®^‘‘^L 

;]l/2 


” "'iLN^W^v’f" — 
bw 

^MLN ^LN 

‘LN * ^NLN' 

sin -xln 

11/2 


AN[jj = <^lnSw*^W cos ct|^J - 


_1 sin .-Lv cos, L^'] 1/2 
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LANDING GEAR EQUATIONS 


PERFORiM THE FOLLOWING CALCULATIONS FOR EACH WHEEL OF THE 
LANDING GEAR WHERE - n = 1 LEFT MAIN GEAR 

n = 2 RIGHT MAIN GEAR 
n = 3 NOSE GEAR 

LANDING GEAR - A/C LOCATION 


= - XCG + ^Gn 
Yn = 

~ " ^CG ^Gn 
STRUT DEFLECTION 


^G6n 

= sin S 

- Zj^ COS 9 - 




= [y^ sin ^ 

+ + r ) 

n n 

(cos ♦- 

1)J cos 9 

^Tn 

^■“DOWN 

h ^ h ^ ^ 

G 0 n G n 

)/(cos J 

cos 5 ) 


RATE OF STRUT DEFLECTION 

^Tn “ “^DOWN ^ ; cos i) + Xj,q - P 

VERTICAL FORCE 

fGZn = ^STn !"'Tn + °STn ^Tn 

NOTE : COMPUTE Fczn ONLY IF h-pn < 0; 

IF h-pn ■■ 0; F^^n = 0.0 & 

REMAINING CALCULATIONS MAY BE SET 
TO ZERO. 
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LONGITUDINAL FORCE: 


= '“o * “l FgZ„ “/i“l 


NOTE: is per. -'nt brake pedal deflect! 


on 


SIDE FORCE: 


F = y F v/!v| 

Sn S GZn ' ' ' 


FORCE AND MOMENT CONTRIBUTION OF EACH WHEEL 


X ^ 
n 

- 

F 

^GZn 

6 (n = 1,2) : 



AX 3 

= ^U3 

COS 

- F^^ sin 
STEER S3 

^ STEER 

p a 

‘G7 3 ' 

It 

c 

^Sn ^ 

F 

^GZn 

V (n = 1,2) ; 



^^3 

= ^S3 

cos 

^ STEER 3 

STEER 

^GZ3 • 


^yn “ ~ ^Sn ^ ^GZn 


nM = 


1N_ = 


= - ^z 

“ n 

X 

n 

+ 

^x 

n 


+ 

r + 
n 

’■‘Tn) 

— r ^ 

n 

Y 

n 

- 

n 



r + 
n 


It 

t 

X 




V 

n 
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AX 


LG 


3 

Z A 

1 


X 

n 


AY 


LG 


3 



n 


AZ 


LG 


3 

Z AZ 

1 


n 




3 



3 

AM. „ = E AM 
LG ^ n 

3 

ANlg = ^ 


E-5 6 
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FUSELAGE AERODYNAMICS 


FUSELAGE INPUT EQUATIONS 
-1 


= Tan 
F 


|w/u) 


ap = sin otp cos 


y = f[j2- + y2 ^^2 

F V 


gp = Tan 


-1 


V/Vu^ + W2 


6p = Sin Bp cos Bj 


qp = l/2p V2 

^FUS ^ 

FUSELAGE WIND AXIS COEFFICIENTS 


(1 + KqI&pP) 

+ 1 a 

. + ■ 


■ + Kg 3’ |S' 

1 

+ Kj 0.' + Kg 

,, 1 1 ^ I i + P 

f‘ "Silg 

+ Kg Sp + KpQ 

.5 H 0 » i 
"F ‘ F ' 

IS UP; 

t> 

n 

Ml 

O 

o 


'DLG 


SPECIAL CONDITIONS 

1, IF ^ 1 (ft/sec) = FUSELAGE AERO = 0.0 & 
' HOLD VALUE OF ip & Bp 


✓ 
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FUSELAGE FORCES AND MOMENT ABOUT A/C C.G. 


F’ 

^AERO 

t ‘-OF 

cos Op + Cj^p sin Op - Cyp sin Sp cosop] qpS^ 

^AERO 

■ I'^YF 

cos Bp - Cpp sin Bp 1 qp S 

w 


^AERO 

■ '-'^LF 

COS Op - cos Sp sin 

“f 




C sin S sin a ] q S 

YF F F F W 



"^AERO 

= [- (c, 

w^’^w^Sf ^F “f " 

^NF ® ^ 

^F ^’W 



F ’ 

'^AERO ^^CG " ^FAC^ 



F’ 

M 

'aero 


, F ’ 

cos Bp] qp Sy c^ + 

^^CG “ ^FAC 

] 



F ’ 

“^AERO ^^CG ~ ^FAC^ 



N^’ 

‘aero 

= [Sf 

cos ap - (c^/b^) C^p sin 

S p sin a p ] 

Sy by 



Y ’ 

^AERO ^^CG ” ^FAC 

] 


F 

''^AERO 

= x"’ 
■'aero 




AERO 

= Y^' 
AERO 

+ -\g 



“^AERO 

= Z^’ 
^AERO 




-^AERO 

-^AERO 

+ a£ 

■'■^LG 



^AERO 

F ' 

= M 

^aero 




‘^AERO 

F ’ 

= N 

'\\ero 

“lc 
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WING ON ROTOR INTERFERENCE 


AVERAGE NACELLE INCIDENCE 




AVERAGE LIFT COEFFICIENT 


Cj^y =0.5 


^^LSRW ^LSLW^ 

(1 — t!rr, q ) 


=MRR ‘ 'WRL ® is ‘ <^LW 


WING INTERFERENCE LOGIC 


1. IF: Umbrellas open, set C^y = 0.0 & £ = -p (1~GEF ) 

/1-M^ 
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ROTOR/ROTOR INTERFERENCE 

POSITIVE SIDESLIP, I.E., V > 0.0 (Logic Required) 
X = 1.5708 - 


^ '^RL 

V* , 

RR/ 


= + T, X + T^x^Jx 


6 V, 


L‘l ■ ‘2 

\ 


RL 


^Irl 


iRL 


V, 


= -tan 


’*R 

V* 

RR 
-1 




^RR 


2p ttR^ 

‘5'^RL 


Vlr + 1.0 


= (iBpl) (.40528 i^J 


^'iLR = 0.0 

NEGATIVE SIDESLIP, I.E., V < 0.0 
X = 1.5708 - 


■PLR 




V* 

LR 


6v, 


-It + t 

L ^ ^ 

X + ' 


^ r \ 



V* 

^ RR^ 

V 

= - 

•tan"^ 

«''LR 



L^RR ^ 


^LR 


2p ttR^ 


LR 

I 

“iLR 
'iLR 

^'iRL “ 

NOTE: v*j^ & v*L FROM WING EQUATIONS 


= (l 3 pl) (.40528 i^i^) ■■■Vlr 



I 
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ROTOR EQUATIONS 


RIGHT ROTOR 


^RR ■ 


-1- I + (Wpp + U 


1 


RR '"RR "^RR 'WRR^ 


+ "iLR 


U 


RR 


RR 


RR 


VU^ + V- + W2 ; p = ^ . 

'' RR RR RR RR |J^r!R 


LEFT ROTOR 


-ia/v^ + (W + U £ ) 

1 \ V PT. ' PT PT MPT/ 


= tan “ RL RL RL WRL 


U 


I 

J 


RL 


iRL 


V 


LR 




KL ^ ''^L ''^RL ' ^LR 


V 


LR 


^ 


ROTOR ANGULAR RATE TRANSFORMS 


RIGHT-NACELLE AXES 


LEFT-NACELLE AXES 


NR 


P COS t r sin 


= P COS - r sin i„^ 

NL NL NL 


= 

^NR 

q + 

^NR 


«NL 

= q + 

^NL 





II 

- r 

COS ij^j^ - p sin 

^NR 


= r cos ipjL 

+ p 

sin 

^NL 


RIGHT 

WIND 

AXES 


LEFT 

WxND 

AXES 





"nr ’ 

pN 

NR 



pR = 
NL 

: pN 
NL 






= 

NR 

qN 

NR 

M 

COS +R‘ Sin ", 

HR NR 

HR 

NL 

= 

NL 

COS 

HL 

R^ 

NL 

sin 

HL 

rH = 

NR 

NR 

cos .% -Q^’ sin ■ 

HR NR 

HR 

rR 

NL 

= R^ 
NL 

cos ' 

+ 

HL 

oN 

NL 

sin 

HL 

NOTE: 

USE 

WIND AXIS RATES 

IN 

ROTOR 

ROUTINE. 
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RIGHT ROTOR 
THRUST 


T ^S + l“ 

C Cos At Cos Bi 

T 2 S +1 

"^ORR ^R 


WHERE: = 0.000679 <(> + 0.000015 *2 

+ 0.0022 + 0.000211 


and 



tan~^ 


u cos a 
0.75 


6.3015U + 5.5816u2 


I 


- 8 y sin a + 1.115 

GROUND EFFECT 


~ “^DOWN ^NR ® 


(L Sin i + Z „) Cos - Y Sin 
S NR CG N 


Cos e 


Dl EFF 
RR 


‘RR 


2R|^|sin{9 + i^^) Cos >1 + .0174j 


IGE 


'■OGE 


RR 


-1 ' 

dJ 

EFF 

RR 


(.1741 - .6216 


(- 

Id 


EFF 

RR 


(1.4779 yRR - .4143) 


+ 1.2479 - .8806 y 


RR 


'Trr 


= c • 

Trr 


"^IGe I 

"OGE* 


RR 


SPECIAL CONDITIONS : IF > 0.283; 


IGE 


'OGEi 


= 1.0 


RR 


or IF 


1^1 - 1 

I'^ige' 

'Oi gpp 1 

i'^OGE, 


= 1.0 


RR 


RR 


9 


✓ 
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^prr ^ ^pqrr 

= 0.00006 + 0.00057 u + 0.000085 u- + 1.12 

- 0.024075 Ct + ’hiCt (0.53 + 0.456 y -39.937 LiC,p 

+ 31.79 Ct) + [0.0115U - 0.03y2 - 0^(3. 4y -Ou^)] 

- 0.22064U (Ct + 0.001971) sin a 


+ (0.3082y - 2.18y2) sin a 
NORMAL FORCE 

dC 


dC, 


c = r ^ _^RR A + ^^NFrr 

+ -^ICR + — ®1C, 




dA 


1 Cl 


dB. 


ICt 


WHERE; C\jFq = Cnfj_ = 0.068y3 sin 2j. + [0.133695-Ct 

+ 73.444b Ct" (1-b) ]K 0 l ■- 


0.6 


For 0.6 

yir ^N] 
dCjaF 
'^^RR 


dA 


ICR 


dCj^F 

~dBT 


RR 

CR 


' sin jc 

for : 

'1 

> 20 ° 





= sin a 

( 10 - 0 , 

.45a°) 

for 0 

1 < 

-j. < 

20 

' = (CisiFj^) ( 1 - 

-0 

. 8 ( b- 

0 . 6 ) ) 





Ct 

^RR 


°NF 2 

^RR 

+ D 

NF 3 

“RR 



+ 

°NF 5 

^RR 

sin 

2 i 

RR 

" ^NF^ 

*“T 

^RR 

+ 

^NF 2 

’■^RR 

+ E 

NF 3 

“RR 



4- 

^NFs 

'■‘rr 

sin 

■‘rr 



+ D.. 


( °‘RAp ) 
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SIDE FORCE 


dCsFpp 

C = C + A 

SF 

RR ORR dA 


dCsF 


RR 


ICR 


B 


dBi 


ICR 


‘ICR ‘^^ICR 

WHERE; ^sFq = 0.00430 y s±n ct - 0.0028827 y " 

+ 0.012 y sin a Ct (90-ii)°) + 2.19y^ “s: a 

„ - 1 r u~yi cos a 1 

Wile re il» = tan ‘ ; =- 

yj_ sin a I 


and pj_ - 


’(u** + c^2) - y" /2 ' 


i/' 


dC 


SF 


RR = 


dA 


ICR 




+ Dsf- ^JRR "^RR 

J 


dCsj,. 


RR = 


dB 


ICR 


E C + E y“+E y +E„„ 
^SF^ Tj^p_ SF 2 RR SF 3 RR SF^ 


■*■ '^RR 2 
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HUB PITCHING MOMENT 


C 


PM 


RR 


dCpM 


RR 


■'PM 


ORR 


dA. 


ICR 


CR 




^®1CR 


B 


IfR 


dCpMpR 

+ 

dQ 


Q 


R 

NR 


WHERE : 


'PMr 


= 0.009950 ;i sin a -0.010960 l;2 sin 


+ 0.0028126 VI sin 2a - 0.0057743 u sin a 


RPM 

298 


+ (1. 802 VI sin a - 7.56 (m sina) ^ Cr 


1000 


dCpM = 1^5 + 


dQ 


0 < M 


0.25 + 7.26 u .2 <'v ^ . 39 
4.1681 -2.79 V V > .39 


dC 


PM 


RR 


dA 


ICR 


" ^'‘^RR ^ ^ 

'^PMs ^'rR ^ *^RR °Pr!c '^RR ^'”R^"“‘o^ 
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hub pitching moment (CONTINUED) 


d C 


PMRR 


= E. 


d B 


ICR 


PM^ TRR 


+E„ y2+E y+E 


PM 2 RR 


PM, 


PM. 


^PM^ ^RR sin aRR + UpMg 


HUB YAWING MOMENT 


'YM 


= C 


RR 


YM 


+ B + 

1C8 dB,,. IM dR 


ORR 


CR 


'ICR 


*^^YMrr R 

~ ^R 


Where : 


For 0 ^ y ^ 0. 37 

Cym “ (0.018369 y -0 . 0007) ysincx -1.2 y^ C.p sin a 


0.00631 -0.002604y -0.004877 (|~ -l) 

\29 8 / 


/ 


RPM _T „ 

1 Trrrr sm ct 

^293 , 


and for y > 0.37 

Cym = (0.01916 - 0.15321 (y -0.5435)2) gi^ a 
- 1.2 y2 Cip sin a 

«^CyM = _ dCpM 


dR 


dQ 
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ROTOR FORCE & MOMENT CALCULATION 
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! 

m0 


= f c puR'+n^ 

Tj^ TRR R 


NF, 


= f, 


4o2 


NFp ^NFRR ^R 


SF 


R 


^SF„ ^SFRR 

IV 


^ ~ ^PMj^ ^PMRR 


N, 


^YMj^ SmRR 


■'RREQ 


*^PRR P^^^^R 


RHPRR “ Iqrreq 


550 { 


LEFT ROTOR FOLLOWS SIMILAR FORMAT WITH SUBSCRIPTS CHANGED , 
THE LEFT ROTOR ALTITUDE EQUATION IS AS FOLLOWS ; 


^LR =-2dOWN + ^NL - ® 

+ [ (L„ sin i„, + Z„) cos 41 + Y„ sin <j> ] cos 0 
S NL CG N 


r 

*! 

iw 


or; 


\r “ ^RR 2 * °°® ® 
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ROTOR FORCE & MOMENT RESOLUTION 


HUB MOMENTS - NACELLE AXES 


LEFT 


- *^LREQ ~ 


\rH = ^HL “ \ ^HL 


- (p sin iNL + r cos Inl) ( KIp^L + Nel Ki Ig ^el) 

^LRH - COS 5hl ” ^L sin 5 hl (KIpi^L ^el Kj^Ig^EL) (<3+iNL) 


RIGHT 


“^RH ®RREQ ■*■ ^R 
*^RRH “ ^R ^HR ^R ^HR 

+ (P Sin i^R + r COS i^R) (KIp^R - NgR KiIeS^er) 

**RRH " ^R ^HR ~ % siri - (KIpHR - NerKiIe^^er) (q + iNR) 

NOTE ; NACELLE AXES ARE RIGHT HANDED SYSTEMS 
* 0 if nor-tilting engines 
= 1 if tilting engines 


>1 ^ 
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RESOLUTION OF ROTOR/NACELLE FORCES TO BODY AXES AT PIVOTS 
LEFT ROTOR 

NL t 

XaERO = (Tl + AXln )cos i^L- sin Inl (NFl cos Chl+ sin ?hl 


,NL 

AERO 

jNL' 

'aero 


- AZ ) 
LN ' 


'hL ■ «HL * “’'ln 


-<’'l * ‘\n> Sl ^'■l 'hL * 




NL’ 


AERO '^LRH 


.NL 


'hl - ‘\l.> 

I 

DU ^ cos ^xiT ^ s xn X 


NL 


NL ^^LRH 


’ NL 

* “m + I's’^AERo’ 


^AERO \rH ■*■ ^-"^LN ■*■ ^^L ^HL ''' ^ 


HL 


• - ^ - I^ r K 2 

’^AERO “ ^NL ^\rh ■*■ ^^LN ^AERO ^ " sin ^*^RH ^^LN^ 

■'■ ^E ^EL ^'eL ^2 
NACELLE EQUATION INPUT - LEFT 

NL 

"nlaero “ “aero ^ Ie "ei “el “2 

GLAS INPUTS - LEFT 


“nLAEEO - «LRH * '■s 'hl * 'hL> 

GLAS 

N - N,,„„ + (SPj^ cos - NFj^ sin 5^j^) 


K 


NLAERO LRH 
GLAS 


2 = 0 if tilting 
engines 

= 1 if non- 
tilting engines 


I i 


i I 

I 
1 


I 
I 
I 
I 

1 1 ' 


T 


i 1 
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RIGHT ROTOR 


= (T + AX' ) cos i + sin i (-NF cos C 
AERO R RN NR NR U HR 


+ SF sin 5 + AZ* ) 

R HR RN 


Y = ^SFpCOs 5 - NF sin C + AY' 

AERO HR R HR RN 


i ■ 


AERO 


- (T + AX' ) sin i + cos i (-NF cos C 
R RN NR NR R HR 


+ SF sin ^ + AZ' ) 

R HR RN 


-AERO 


(X + AJ?' ) cos i + sin i (N + L y’^^ + ANi„) 

RRH RN NR NR RRH s AERO RN' 


^AERO = \rH ^ ^^^N * ^ ^HR ‘ "^R ''s ^HR 


's ^^RN " ^E ^ER ^ %R ^2 


^AERO ^NR ^’^RRH '^^RN ^AERO^ “ ^NR ^^RRH ^^RN^ 


^E ^ER ^ ^ER ^2 


NACELLE EQUATION INPUT - RIGHT 


M = + I n r N K 

NRAERO AERO E ER ER 2 


GLAS INPUTS - RIGHT 


'^NRAERO ” ^RRH ^HR “ ®^R ^R^ 

GLAS 


*^NRAERO “ ^RRH ~ ^®^R ^HR ^R^ 

GLAS 
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WING VERTICAL BENDING 


RIGHT WING TIP DEFLECTION 


- , iAERO + Y 

RT m N ^ 

- _ ^IRO + Y p 

RWAC m WAC 

h, = K + K + K - K a - K a 

Wj^ AERO W2 AERO W3 «<-AERO RT W5 RWAC 




Where Ah-, Is the difference of h-, between time frames and 

^R 1r 

At Is the time frame* 

RIGHT WING A.C. DEFLECTION 


h, = K„ z”p' + K„ Z^",- + K„ - K„ a - K„ 

IrWac ”6 AERO Wy AERO Wg ^AERO Wg RT W 


,RW' 


■NR' 


10 


RWAC 


h., = Ah, /At 

■‘■RWAC ■‘■RWAC 

Where: Ah 

and At Is -vhe time frame. 

FORCE AND MOMENT EFFECTS 


1 Is the difference of hi between time frames 

■‘RWAC ■‘RWAC 


;nr 

'AERO 


,RW 

'aero 


VNR 

AAERO 


= -2C„, 0), 


-2C. 


= -2C u) 


NR 

0)2 

z>" 

+ 

0 )^ 

7 NR' 

AERO 

W 1 

AERO 


VI 

AERO 

RW 

(1) ^ 

ZRW 

+ 

03 ^ 


'aero “ 

W 2 

AERO 


W2 

AERO 

- 

0)2 

^R 

+ 

03 ^ 

yNR' 

■AERO 

W 3 

AERO 


V 3 

*^AER0 



:#> 




|j 
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LEFT WING DSFLECTION 

- _ ^aero . V a 

ar ^T, tr 


"LT 


m 


N 


t... 


ZaERO 


aWAC 


m 


^WAC P 


h - K + K 

"iL W, ^AERO '^AERO 

1 2 


K - K a - K a 

^W^'^AERO LT LWAC 


/ 


hiL - Ah]^L/At 


Where: Ahj^j^ is the difference of between time frames and 

At is the time frame. 

LEFT WING A.C. DEFLECTION 

. _ NL' LW’ _ y,NL' _ - _ - 

IL W, AERO ^W, ^AERO ^W„ '^AERO LT ^W, „ LWAC 

WAC 6 7 8 9 10 


h. 


= Ah.,T ■ /At 
"WAC 


i; 


i 

i 


Where: Ahii is the difference of hii between time frames 
■^^WAC ^‘"WAC 

and At is the time frame. 

FORCE AND MOMENT EFFECTS 


,NL 

"aero 

•LW 

AERO 


= -2C.„ CO. 


■25,.,. CO,,, Z 


^NL 


AERO 


= -2 5 0) 




NL 

- 0)2 

2NL 

•f 

OJ ^ 

ZNL- 

AERO 

W1 

AERO 


W1 

AERO 

LW 

0 

LW 


0 

LW' 

AERO 

'■^W2 

^AERO 

+ 

^W2 

^AERO 

NL 

- 0)2 

^!TL 

+ 

9 

0) *" 

.j^L' 


W3 -^ERO 


W3 -^ERO 


FORM , z’"'^ , and 

AERO AERO "^AERO 
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WING TORSION 


LEFT WING TWIST ?iT TIP 


^0t ®tLW “ *^NLACT “ ^EL ^ 


+ q 




SLW 


^MO ■ ^TSLR^ 
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* - =ISLR> 9s.w 

\ dc^ <H//\ 6n / 


^^®tLW ■*■ ^““lWRIGID^ 


RIGHT WING TWIST AT TIP 


• r 

I V 


i \ 
1 
l 
l 
1 
I 


^9t ®tRW “ *^NRACT " ^E ^ER ^ 




WHERE; 


+ a 


(1 

^SRW 

MO 

2 

+ (1 - 

^TSRR^ ‘^SRW 


Uq 

tRW 


'TSRR 


) 


^2y^^HwcM ^ ^WAC 


dc. 


w 



RWRIGID 


C =C+C6+C62 
MO 1 2 F 3 F 


') 


^Lg^w \ 

6n J 


fl - WAC „ 

“tLWAC " ~Y^ ”tLW 

A - ^WAC „ 

°tRWAC " ”tRW 

NOTE; If umbrellas are open; set terms containing 
qg (1 - equal to zero. 


E-74 





0 


L 


f 

[ I 


TOTAL FORCE AND MOI4ENT SUMJ4ATION ABOUT C.G. 


X 


.NL 


.NR 


AERO 


AERO 


X — + X — + X^ + + x^^ + x'^ 

AERO AERO AERO ' AERO AERO AERO 


,NL 


,NR 


,LW 


Y"" + Y“" + Y* + y“" + + Y^ 

AERO AERO AERO AERO AERO AERO 


.NL 


.NR 


.LW 


,RW 


2aER0 “ ZaEro + 2Xero + Zaero + ^AiRO + Zaero + ZaERO 
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1 . 


I ; 
f ■; 


I ^ 


i 


: I 


_ ^nl ^ .t 

AERO “-^AERO +*^AER0 '^‘^AERO ''■«>^AER0 ■’’-^AERO 


NR _ NL , NL NR , 

+ ^N '^IaerO ^IaerqJ + ZcQ (Yaero + ^AER0‘ 


I L 


.NL 


.NR 


-W 


^AERO “ ^AERO ■*■ ^AERO + ‘UeRO ^AERO *^AERO 

+ Xcg(Zaerp + “ ^CG (^AERO + ^AERO^ 


^AERO “ Naero + Naero + ^aerO ^aero ^aeRO 

+ Yr (^AERC “ ^AERO) " ^CG <^AERO ^AERO^ 


f ' ■ 


f 

i 


/ < 

^ 1 


I i 
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BASIC EQUATIONS OF MOTION 


PRELIMINARY CALCULAT IONS 


FUSE 


W.R.T. A/C C.G. 


Xf ~ ^ 


CG 


Zf - hf - Zqq 


WING C.G. W.R.T. A/C C.G. 


X = Jl - 
w w CG 


Z = h - 
"w w CG 


NACELLE C.G. 's W.R.T. A/C C.G. 


\ ^ COS - X) -X 


2r = sin - \)-ZcG 


\ = i cos (ijjL " '^CG 


2l = -^CG 


' .1 

f 

1 \ 

D238-10002-1 

1 i 

I 

(T# 

I 
1 
I 
I 
1 
i 
1 
I 
I 
I 
I 


I 



L 


L 
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PRELIMINARY CALCULATIONS 


INERTIA TERMS 


^ ?! I^^) + I^''^ + 21' 

^ ij ij ij ij 


^xx ” ^ ^xx^ ~ ^xx^ (sin + sin 

"^xz ^ ^NR '*' ^ ^ ^ 


+ m^ hf Zj + m„ h„ 


- ^.m^jIZj^ sin (iMt> - ^) + sin (i^T " 


"NR 


NL 


'XX 


•xz 


^zz “ ^yy 

' ^XZ^ ■*■ Ixz^ 1/2 dxx 


- 1 ^ 2 ) (sin 2 iNR+ sin 2 i^L/ 


L 


+ I* (cos 2 i„_ + cos 2 i^ )+(m l Z + m l Z ] 
XZ NR NL iff vjww 

+ nijjJ. cos (i(jR " ■*■ Z^^ cos (i^L ~ 
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L 



4 -a,, h , 






4 

* 

I 

5 ^ 
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INERTIA TERMS 

"yy ' ^ 'w* * ”f “f ’‘t * '*f * "w “» ’‘w * hw *w> 

+ nijj n t Xj^ cos (a.jijj) ~ X) “ Zj^ sin * ')! 

“^yy “ ^xx ■ ^zz 
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INERTIA TERI4S 


I =Z + (I* - I’ ) (sin^ + sin^ i ) 

^ -y 9 - 7 . XX Z Z 


. 2 ^ 


z z . z z 
k 


XX zz 


+ I- (sin 2 i + sin 2 + 2 rOj, Yj, 


+ nif Jlf Xf + iy Xy 


+ mjj H [Xr cos (ijjR - X) + Xl ccc (i^l 




Jzz = ^yy ■ ^xx 
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ROLL EQUATION 

I.... P + P'3) 


XX 


XX 


+ 2,m Y < cos (i '^) 
N N) NR NR 


1 COS '.X 

NL 


i -mI 

NL J 


AERO 


PITCH EQUATION 


i„ q - i.<z 

2 

(P - 

2 

r ) 



^NR ^ 

I 

i I' + 

^ yy 

amj,[-ZR 

sin 

^^NR 

-X) + 

- ^nlI 

I' + 

yy 


sin 

^^NL 

-X) + 

N. 

■'■ ’^AERO 






(iwp 




YAW EQUATION 


I^z r = -J^z pq - (rq “ P) I 


'zz 


xz 


- Yn { Inr sin Unr " ^NL ^in Unl "^)| 
+ N^ero 
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RIGHT NACELLE ACTUATOR PITCHING MOMENT EQUATION 


^NRACT " ^NR 


'yy "S <1 - 


- <1 - 7 


1) r - pr cos 2 -A) + q 

+ (r^ - p2) sin -X) cos -X) 


- (r^ - p2) 




zz ^NR ^NR 


-I' q 

yy 


*AE«0 <^»R - “‘aero '“-HR 


+ i 

- *'"’n [ sin (iNR-^> 


- M 


- (P 


- 

+ rq) cos 


- M 


^NRAERO 

LEFT NACELLE PITCHING MOMENT EQUATION OBTAINED BY CHANGING 
SIGN OF Yjj AND CHANGING SUBSCRIPT PROM R TO L. 

NOTE ; THE ABOVE EQUATION MUST BE CALCULATED FOR WING TORSION 
CALCULATION ONLY. 
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V ' • 


1 


MOTION OF A.C. MASS CENTER 


u = 

^AERO _ 

g 

sin 

e 

- qW + rV 


m 




# 

^AERO 



0 

sin (ji - rU + pW 

V = 

T 

g 

cos 


m 




. 

^AERO 





w = 

+ 

g 

cos 

6 

cos ij) + qU - pV 


m 


EULER ANGLE CALCULATION 

i|) = (r cos <1> + q sin <)))/cos 0 
0 = q cos ij) - r sin <)) 

(ji = p + Ij) sin 0 


I 
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I 

1 

1 

1 

I 

1 

I 

I 

I 

I 

I 

1 

I 

I 
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AIRCRAFT CONDITION CALCULATIONS 


GROUND TRACK 


NORTHWARD VELOCITY 


Xjjorth ” ^ ^ il> + V (sin ((> sin 9 cos 4> 


- cos (ji sin i|i) 


+ W (cos (t> sin 0 cos ij» + sin 41 sin t(») 


EASTWARD VELOCITY 


Yeast “ ^ ® 4' + V (sin (p sin 0 sin + cos 

4 > cos 4') 


+ W (cos 4> sin 9 sin 4< - sin 4' cos 4") 


DOWNWARD VELOCITY 


Zr.rttixT“ “ U sin 9 + V sin 4> cos 9 + w cos 4> cos 9 
DOWN 
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PILOT STATION ACCELERATIONS (BODY AXES) 

^XPA “ + (q + Pr) (Zp^ - Z^q) 

+ (q2 + r2) (X^g - + Yp^ (pq - r) 

m 

” ^ *5 ^CG ” ^CG 

^YPA “ + (P - qr) (Zgg - Zp^) + (r + pq) (Hp^ - 

” ^PA + P^) +2 (pZgg - rX^g) 

^zpA = ^AERO + (q - pr) (X^^ - £p^) + (p^ + q^ : (Z^^ - Zp^) 

m 

+ Yp^ (p + qr) + 2qX^.g - Z^,^ 


PILOT STATION VELOCITIES (BODY AXES) 


UpA 

il 

a 

+ 

‘^ZpA 

1 

> 

VpA 

> 

II 

+ 

^^PA 

P^PA 

WPA 

W 

Wp 

+ 

P^PA 

1 

> 


GUST MODEL 

The gust model will be that represented by NASA-AMES program 
NAPS-80. The output of this program, in the form of gust velo- 
city components Ug, Vg, Wg, pg, qa, rg will be added to the 
aircraft velocity components in clear air as follows: 


U = U' 

+ 


II 

+ Pg 

V = V 

+ 

Vg 

11 

- qg 

II 

+ 

Wg 

U 

II 

U 



1 

I 

I 

1 

1 

1 

1 

I 

I 

I 

1 

I 

1 

I 

I 
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PRELIMINARY CALCULATIONS PREPROCESSOR 


g 

= 

32.174 


m 

= 

[w'+w’t 



= 

/2g 

mf 

= 

[w* + 

W^T + «VT '^C 


= 

[w- + 


1 * 

f 

= 

[ (FS)p 

- (FS)j,1/12 

NT 

= 

[(FS)p 

- (FS) ]/12 

HTCG 

5f* 

— 

[ (FS) 

- (FS) ]/12 

VT 


P 

vrcG 

Z ’ 

_ 

[ (FS) 

- (FS) 1/12 

PA 


P 

PA 

Z * 

c 

= 

l(FS)p 

- (FS)^]/12 

■ .1 


£' = fw'Jl'+w' I' +w’ Z' + w' +w'il']/(m^g) 

f '• f f ht ht vt vt cr pa C C f 

= [ (FS)p - (FS)^j]/12 

‘fuel = 1 

‘iF = I(FS)p - (PS)^^1/12 

* "fuel SuEL * ”iF 
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'f 

= [(WL)p - 

(ML) j, )/12 

, 1 

HT 

- [(WL)p - 


p 1 

'VT 

= [(WL)p - 


‘PA 

= ((WL)p - 

(ML) 1/12 

PA 

, 1 

c 

= [{WL)p - 

(WL)^ ]/12 


h, = 1/(32.174 .f) 4 z'l 


z * 


[(WL)p - (WL)jj]/l2 


FUEL 




- t(WL)p - (WL)^pl/12 


w 


1/(32.174 m ■) fw'z' + W z' + W z' ] 
^ ^ W W FUEL FUEL NF NF' 


X. 


WAG 




“•WAC 


= [(BL)„^^l /12 


WAG 


((WL)p - (WL)^^^] /12 


[ (BL) 1/12 
N 




D238-10002-1 1 


^HT 

= [ 

(FS)p 

- (FS)^^^] 

1/12 

^HT 

= [ 

(WL)p 

- (WL)jjp] 

1/12 

XvT 

= [ 

(FS)p 

- (FS)^pl 

1/12 

Zyi 

= [ 

(WL)p 

- (WL)y^] 

1/12 

A 

= 3. 

14159 

R2 


Y 

^WAC 

= r 

W 

^®1'^WAC^ 1/12 



Xg 2 = = [(FS)p - (FS)g 2 ] 1/12 

Zg2 = 5^g^ = [(WL)p - (WL)g2] V12 
Yg2 = [ (BL)g2] 1/12 

V = -V 

^G 1 ^G 2 

^G 3 = ° 

Yp^ = [(BL)p^] 1/12; POSITIVE FOR PILOT IN RIGHT SEAT 
^fAC ~ [(FS)p - (FS)^^gI 1/12 

ZfAc = nwDp - (WDfAC^ 1/12 


I 

f 

( 
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c/2 

= [(FS)p - 

(FS)^/2^/12 

G3 

= [(WL)p - 

(WL)g3]/12 

G3 

= [(FS)p - 

(FS)g3l/12 

NF 


1/12 

HT 

[(BL)hTCG 

1/12 

W 


1/12 


^FUEL f ^®^^FUELCG^ 


INERTIA CALCULATIONS 

n^. 

= 

£f- 4 

6^. 

= 

hf- 

n^T 

= 

*■£■ *'HT 


= 

hf- z^^ 

nvT 

= 

'VT 

•^VT 

= 

"VT 

'’CR 


*'f" *'PA 

«CR 

= 

^PA 
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)l, - H' 

f c 


« - 
i ' 


6' = h, - z* 

c f c 


l(f) = i(wf.) ^ j(HT) j(VT) ^ j(CR) ^ j(C) ^ 

yy yyo yyo yyo yyo yyo f' 


+ W' /g (n'2 + 6-2) + w* /g (n'2 + 5’2 
f f HT HT HT VT VT VT 


XX 



'"iR " 

6'2) 

CR 

+ W'/g 

C 

(n'2 

c 

+ 

6'2) 

C 

(Wf • ) 

+ I<«> 

+ I 

(VT) ^ 

(CR) 

+ 

j(C) 

xxo 

XXO 


XXO 

xxo 


■^xxo 

W /g 
HT ^ 

( 5'2 + 
HT 

) 

HT 

+ W» /g 

6 ’ 2 

VT 

+ 

W* /g 
CR ^ 


CR 


(f) 


zz 


+ w'/g 5 ' 2 

c c 

!.<?.■> + .. I.<« * W^/g nif 


zzo 


zzo zzo 


zzo 


zzo 


«iii/9 (n;j + + w;yg t w;„/g n;; + w; /g 


VT = VT CR ^ CR 


(f) (Wf 

I = I 

xz xzo 


HT' 


(HT) 

(VT) 

(CR) 

(C) 

+ 1 

+ I 

+ I 

+ I 

xzo 

xzo 

xzo 

xzo 

HT '^HT ^’^VT'^^ 

r, 1 r 1 

Hvt ^vt 

* “iR/? 


+ W’/g n^. 6^, 


, I r I 

’CR ^CR 


+ wyg n’ 6' 


H' 


w w 


w 


I ' 
w 


A', = h - z' 

w w w w 


H' , = n - e' 

w'FUEL w FUEL 


^w'FUEL * ^FUEL 


E-89 
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U 1 = i — £ * 

w’NF w *^NF 


' NF ^NF 

^ -'w/^ 

+ w'FUEL/g (H^tpuEL ^w’FUEL^ ^ ''’nF^^ <”w’NF ^w’NF^ 

j(w) ^ ,(Ww) ^ ^(WFUEL) ^ ,(w-NF) ^ (,.2 ^ y2j 

XX xxo xxo xxo W w w 

+ w'FUEL/g ^FUEL^ ^ NF'^^ ^^w'NF ^ ^NF^ 

(w) ^ (w'w) ^ (w'FUEL) ^ j.(w'NF)_^ (H'?,+ Y^) 

^zz ^zzo zzo zzo w wv 

+ w'FOt!L/g + ^PUEL^ ■*■ '*' ^**w'NF * ^HF^ 

i(w) ^ j(»'w) ^ F wyg h;,„ 4„,„ 


xz 


xzo 


+ W'FUEL/g H;,pupL ^^FUEL '^'nf/^ “w'NF ^w' 


NF 


E-90 
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APPENDIX F 

This appendix contains the nvunerical constants and functions 
required by the equations presented in the preceding pages . 

The data is listed by reference to the page number in Appendix 
E where the numerioal constant or function first appears. 





INPUT DATA 


D238-J, 0002-1 


PAGE NO. 
E-7 


QUANTITY 

VALUE 

UNITS 

1 

'=’STEER 

20. C 

degrees/'inch X 

^•^RUD 

-8.0 

degrees/inc J 


1.0 

I 


20.0 

rad/sec 


1.0 

1 


l.O 

j 

«L 

35.5 

rad/sec 


0.18 

I 


1.0 

1 


0.0 

1 


-3.33 

•i 

deg/inch 

Schedule A 

to be determined 

i 


Schedule A 
" B 

" C 

" D 


to be determined 


Gains and 'Schedules 
to be determined 



PAGE NO. 
E-9 
E-13 


E-15 

E-16 




D238-10002 

QUANTITY 

VALUE 

UNITS 

Gains and Schedules 


to be determined 
Engine Data 

See page A-17 

— 

WDTIND 

0.0 

— 

SHP* 

4120.0 

HP 


1.0 

— 

NIIND 

0.0 

— 

Nt /n? 

I max' 1 

1.0 

— 

NieiND 

1.0 

— 


1.115 

— 

QIND 

1.0 

— 

QMAX/Q* 

1.0 

— 

Ntt /n?t 
I%AX II 

1.214 

— 

N* 

27932 

RPM 

(Nij/Nii ) 

^^MAX REF 

.7547 

— 

^REF 

27.436 

rad/sec 


2.5 

deg/sec/rad/sec 

^2 

2.66 

deg/rad/sec 


0.05 

deg/sec/deg 


150378.0 

slug-f 

K 

-1.0 

— 

^tr 

1.0 

— 


F-3 
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PAGE NO. 

QUANTITY 

VALUE 

UNITS j[ 


Schedule A 
" B 

" C 

" D 

•' E 

II pi 

" G 

See pages A-19 

It 

tl 

tl 

tl 

II 

II 

to A-22 

1 


" H 

II 

1 

E-17 


-12.0 

degrees J 


mf 

1397.5 

slugs V 



0.091288 

ft * 


mw 

560.82 

slugs J 


^'W 

-0.83235 

ft 

'1 


m 

2323.27 

slugs 1, 


a 

7.4 

I 



182.48 

slugs 


X 

0.0 

de. 1 

E-18 


7.0896 

ft 



0.10271 

ft •* 

T 

E-19 

^WAC 

0.54 

ft ^ 


'•^WAC 

17.51 

1 


^WAC 

1.4167 

ft 



36.55 

“ 1 



8.5 

ft 

1 

E-20 


2.0 

degrees 


F-4 




1 ' 

L, 
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f 

I 

1 

1 

r 

f 

u 

PAGE NO. 

QUANTITY 

VALUE 


UNITS 

f 

E-22 

Zht 

-15.583 


ft 

L 


XjjT 

-49.33 


ft 

1 

r* 

t 

i 


ZvT 

-5.25 


ft 

t 

1 — 

( 


XvT 

-40.4167 


ft 

» 

[. 

E-23 

Solutions of 

See Page 

A-23 


{ 

1 

• 


Quartic 





j 


A 

2485.0 


ft" 

i 

4 

I 

i_. 

E-24 

^WRR 

See page 

A-24 

rad 

i 




It 

A-24 

rad 

' 



WLR 




t 


PC 

4.0 


ft 

i' 

i 

t 

t 

1 , 



0.54 


ft 

i 



D 

56.25 


ft 

\ 


^iN 

75.0 


deg 

1— 

E-25 

Cw 

10.23 


ft 

( 




Sw 

747.5 


ft" 

t 

} 

*- > 


^Law 

4.1832 


rad”^ 

! 

f 

i 

»• 

1 

• . 

E-28 

^7 

0.0269 


deg“^ 

*e 



«2 

22.22906 


deg 

1 

• • 


^8 

-2.437137 










f 

¥ 

• 


ag 

0.20607 


deg"“ 

S, 

r 

•* 


aio 

-0.003128 


deg "2 




63 

29.786 


deg 




^11 

0.442188 


— 


* 

] 

4 


/ 


I 

; 1 
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E-29 




1 





\ 

\ 

4 



m 

D238-10002-1 

1 

' 1 
1 

• 



■ ! - 

QUANTITY 

VALUE 

UNITS J 



ai2 

0.0263 

deg"l 1 

^ 1 ; 

^13 

-0.000338 

deg“^ ■ 

' 


az9 

0.000011 

deg"^ 1 

\ 

1 i 


330 

0.00003 

deg ^ 

i 

k j 


^5 

30.0 

deg J 

1 1 

1 ; 

^31 

-0.03105 

1 

1 i 

1 . 

a32 

0.001946 

deg~^ I 

1 

: 

1 

ai4 

-0.00911 

deg"^ ] 

1 

1 ; 

6<?pi 

30.0 

deg - 

1 

1 ■- 

^15 

0.066537 

J 

1 

i ; 

^16 

-0.016342 

deg“^ 1 

I 

\ 

^17 

0.00014 

9 * 

deg”^ 


1 : 

1 

1 

^0 

-0.000088247 

deg“^ j 

1 



0.000008596 

deg”^ 

t 

V 

Fl 

1.003412 

J 

I 


F2 

0.011163 

deg“^ 1 

1 

i 

F3 

0.002168 

deg“^ * 

1 



20.1665 

deg f 

J 

I 



-0.756323 

n 


’ : 

^5 

0.185684 

deg“^ 

1 



-0.002159 

deg"^ 

r 


a 

14.0 

deg 

1 


0 



■ . 


-0.07 

— 

f 

1 



deg 




30.0 


; 


I 


F-6 
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PAGE NO. 

QUANTITY 

VALUE 

UNITS 


32 

11.9 

deg 


33 

-15.0 

deg 



-0.06 

— 



-16.8 

deg 


^6 

0.0888 

— 


Cl 

0.073 

deg“^ 


323 

-2.051491 

— 


324 

0.215913 

deg“^ 


325 

-0.005276 

deg~^ 


Cdo„ 

0.0152 

— 


^26 

0.0036 

— 


^27 

0.054313 

— 


^28 

0.0269 

— 

E-30 

^18 

11.32 

deg 


®20 

-1.910958 

— 


321 

0.211969 

deg”^ 


322 

-0.005391 

deg“^ 

E-31 

319 

10.92 

deg 

E-33 

b 2 

-0.035178 

— 


b3 

-0.000244 

deg”^ 


b4 

0.0 

— 


^5 

-0.010154 

deg~^ 


be 

0.000081 

deg“^ 


F-7 
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r .i ... 

i * , 

i 

PAGE NO. 

QUANTITY 

VALUE 

UNITS ] 

f' »; 

1 ' a '■ 

t ' 

E-33 

C, 

1.2335 


i : 


LmAX 

1 



{a^/a)„ 

See page A-25 

J 

1 

k ; 

1 

E-37 

K 20 

-0.0908 

rad-l ] 

Iji 


K 21 

-0.04 

rad“^ 1 




J 

[ 


bw 

73.1 

ft 

h 1 





1 . 


V 

1.0 

1 

} 




J 



Yac 

17.51 

ft 

1 ' 


K 22 

0.02 

rad“^ j 

i 1 

1 

1 

1 

A ^ 


Kn 

1.0 

1 




J 

i , 

E-38 

fe 

220.0 

ft^ 

i 


U 


1 



Xc /2 

-1.1408 

ft J 

. 5 


Kdi 

0.0 

I 

! 


T 


J 



Kp 2 

0.0 




T 





Kd3 

0.05 

1 



T 


i 

1, 


^ Ml 

0.0 

— 



T 


1 

1 


^M 2 

0.0 

1 

1 * 


T 





Km3 

0.0 

L 

1 ' 


T 


V 


E-39 


0.05 

i 

1 

■ 


T 




^M4 

0.0 

a 

1 


T 

J ■ 


I 

•* 

F-8 

i 

f 


I 
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PAGE NO. 

QUANTITY 

VALUE 

UNITS 

E-42 


See page A-26 

deg 


de/da 

See page A-26 

deg/deg 


^HT 

0 

deg 

E-43 

■^HT 

0.52 

• 


“^’^STALL 

17.15 

deg 


^“ht 

0.0694 

deg“^ 



See page a-25 

— 


‘^DOrt 

0.00734 

— 


ARrt 

5.384 

— 


^HT 

0.8 

— 

E-46 

da/d 8 

-0.557 

— 


Tvt 

0.46 

— 


°*'^^STALL 

26.9 

deg 


CtvT 

®VT 

0.0613 

0.88 

deg~^ 


‘^D0„ 

0.00615 

— 



1.536 


E-50 

^HT 

1.0 

— 


^HT 

0.0 

deg 


®HT 

227.5 

ft^ 


Hvt 

1.0 

— 


Syx 

278.0 

ft" 


F-9 
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PAGE NO. 

QUANTITY 

VALUE 

UNITS 

E-52 

for a^< .5236 rad 
^DON 

0.00268 



^30 

-0.001908 

— 


«31 

for ajj>.5236 rad 

0.061849 

— 


^DON 

-.102244 

— 


^30 

0.271593 

— 


^31 

-0.077763 

— 


K 32 

0.0546 

— 


^MON 

0.0 

— 


^34 

0.02532 

— 


^35 

-0.00206 

— 

E-53 

^36 

-0.0546 

— 


K 37 

0.0 

— 


«36 

-0.0546 

— 


K 37 

0.0 

— 


'“NORN 

0.0 

— 


00 

u: 

-0.001206 

— 


^39 

0.0 

— 


*"NOLN 

0.0 

— 


^0 

0.001206 

— 


^41 

0.0 

— 


F-10 
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1 1. 

PAGE NO. 

QUANTITY 

VALUE 

UNITS 

1, 

t f 

E-54 

^Gl 

-7.41 

ft 

; i 


Xg2 

-7.41 

ft 

' j 


^G3 

33.33 

ft 

; - 


Ygi 

-7.33 

ft 

. i 
\ 


^G2 

7.33 

ft 

\ 


^G3 

0.0 

ft 

1 1.. 

( 


^G1 

13.5 

ft 

i 


^G2 

13.5 

ft 

L 


^G3 

14.0 

ft 

i 

I 



1.375 

ft 

* 


^2 

1.375 

ft 

* 


^3 

1.0 

ft 



^STl 

20000 

Ib/ft 



KgT2 

20000 

Ib/ft 



^ST3 

20000 

Ib/ft 

« 


°ST1 

1200 

Ib/f t/sec 

1 


‘^ST2 

1200 

Ib/f t/sec 

• 


°ST3 

1200 

Ib/f t/sec 

^ > 

E-55 

M 

0 

0.03 

— 

j 



0.005 

— 

i 


U 

s 

0.5 

— 




1 



PAGE 

E-57 


E-58 

E-60 


NO. 


i 
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1 

QUANTITY 

VALUE 

UNITS i 

^DOF 

0.01288 

1 

Ko 

65.14 

rad - 3 * 

K2 

0.4 

J 

Kl 

-0.00553 

rad"^ 

^^DLG 

0.0401 

J 

K3 

1.26 

rad”^ 

K4 

0.0 

rad”^ 1 

K42 

0.0261 

J 

K7 

-0.361 

rad--^ 

KS 

-0.0987 

J 

^MOF 

-0.00455 

1 

K 5 

1.432 

rad -1 J 


-0.494 

rad-^ 1 

^^TALG 

-0.0568 

1 

^NOF 

0.0 

1 

K9 

-0.0051 

rad-l 

KlO 

0 

• 

0 

rad"^ 1 

^FaC 

5.75 

- ] 

^FAC 

1.4167 

ft 

* 

1 

Ti 

0.2434 

rad“^ 

T2 

-0.483 

rad”^ 

T3 

0.5208 

rad"2 

R 

28.125 

ft 


F-12 
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PAGE 

E-62 


E-63 


E-64 


QUANTITY 

VALUE 

UNITS 


0.1 

sec 

"2 

0.1 

sec 


1 — 1 

Q 

0.00425 

deg"^ 

CM 

a 

0.0014483 

deg“^ 


-0.0000734 

deg"l 

^NF, 

4 

0.00002175 

deg”^ 


-0.0006 

deg"^ 


-0.0245 

deg"^ 

%F2 

-0.0017028 

deg“^ 

^NF3 

-0.0010492 

deg“^ 

Enf^ 

0.0000425 

deg“^ 

^NFj 

0.0017892 

deg~^ 


°SFi 

0.0245 

deg“^ 

°SF2 

0.0017028 

deg”^ 

°SF3 

0.0010492 

deg“^ 

°SF4 

-0.0000425 

deg"^ 

°SF5 

-0.001735 

deg“^ 


F-13 




PAGE NO. 


E-65 


E-66 


E-67 




quantity 


VALUE 


D238-10002-1 

UNITS 


DpM, 

.002 

deg"^ 

•1 

DpM, 

-.00072556 

deg"-*- 

DpM, 

.00111967 

deg"^ 

■1 

J 

.0002094 

deg"-*- 

DpMc 

.00036524 

deg"-*- 

J 

DPM. 

-.00007296 

deg"^/ra<3/ sec 

O 

19.1637 

rad/sec 


E 

-.0025 

deg"^ 

E 

0.0004375 

deg"^ 

^PM2 

0.0000729 

deg"^ 

1 

3 

-0.000111245 

deg 

EpM. 

k 

^PM, 

0.00063045 

deg~^ 

3 

1? 

-.00006809 

deg/rad/sec 

EpM^ 





-.0025 

deg"^ 

°YM, 


, -1 

n 

.0004375 

deg 

*^YM2 

n 

.0000792 

deg”^ 

°YM, 


-1 

n 

-0.000111245 

Cleg 

^YM^ 

n 

0.0005 

deg~^ 

n 

-0.00007296 

deg"^/rad/sec 

-"•e 




I 

1 
i 
I 
1 
I 
I 
1 
1 
I 
1 
J 
1 
1 
1 


F-14 



I 
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PAGE NO. 

QUANTITY 

VALUE 

UNITS 



-0.002 

deg”^ 


EyM2 

0.00072556 

deg“^ 


EYM3 

-0.00111967 

deg“^ 


EYM4 

-0.0002094 

deg“^ 


^YM5 

-0.0004702 

deg”^ 



0.00007296 

deg“^/rad/sec 

E-68 


1.0 

— 



1.0 

— 



1.0 

— 



1.0 

— 


'SF 

-1.0 

— 


‘sFl 

-1.0 





1.0 

--- 



1.0 

— 


'yMj 

-1.0 

— 


fYM^ 

-1.0 

— 


^QR 

-1.0 

— 


^QL 

-1.0 

— 

E-69 


0.40 

slug-f t 

E-72 

^W1 

3.707x10"5 

ft/lb 


^W2 

1.023x10"5 

ft/lb 


^W 3 

1.22x10"4 

ft/lb 


F-15 



1 


PAGE NO. 


E-74 
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QUANTITY 

VALUE 

UNITS 

K 

1.09x10"2 

slugs ft/lb 

W4 



■'«s 

4.66x10"^ 

slugs ft/ lb 


1.004x10"^ 

ft/lb 


:.49x10"5 

ft/lo 

oo 

2.37x10*5 

ft/ lb 

KWj 

2.94x10“^ 

slugs ft/lb 

”io 

1.52x10*2 

slugs ft/lb 

?W1 

0.5 

” “ 

<^W1 

16.33 

rad/sec 

^W2 

0.5 

— 

0)^2 

16.33 

rad/sec 

^W3 

0.5 

— 

<^W3 

16.33 

rad/sec 


4.90x10® 

ft Ib/rad 

"'=««=/4 

-0.002 


dCL 



Cl 

4.1832 

rad"^ 

a 



Cl 

-0.0352 

— “ “ 

C2 

-0.010154 

deg*^ 

C3 

0.000081 

deg*2 


F-16 



REFERRED POWER TURBINE SPEED Njj-//0/N 





REFERRED FUEL FLOW wf/6/6/SHP 













OUTPUT GAIN RATIO 
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in 

VD 

O 

iH 

CM 
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